Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,464,601 G→A 5.8% T136M (ACG→ATG)  mlaA ← ABC transporter maintaining OM lipid asymmetry, OM lipoprotein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,464,6010GA 5.8% 219.5 / 2.8 69T136M (ACG→ATG) mlaAABC transporter maintaining OM lipid asymmetry, OM lipoprotein component
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base A (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

ATTTACCCACAGACATCGGCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCGCCCTGCAA  >  NC_000913/2464455‑2464743
                                                                                                                                                  |                                                                                                                                              
aTTTACCCACAGACATCGGCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGta                                                                                                                                              >  8:73948/1‑149 (MQ=255)
 tttACCCACAGACATCGGCCAGGTCAGCCAGGAAAGCACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGTCCCATAACCCACGCCATAATGACCAAGCGtac                                                                                                                                             >  8:34372/1‑149 (MQ=255)
      ccACAGACATCGGCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTAcc                                                                                                                                         >  2:286295/1‑148 (MQ=255)
       cacaGACATCGGCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCg                                                                                                                                        >  1:297071/1‑148 (MQ=255)
         caGACATCGGCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCCCGCAGCGTGAAGCTACCGAAGACCGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAg                                                                                                                                     >  2:267197/1‑149 (MQ=255)
              aTCGGCCAGGTCAGCCAGGAAAGAACCGGGTAAAACCCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACACACGCCATAATGCCCAAGCATACTACCGAAGCGGt                                                                                                                                 >  3:182599/1‑148 (MQ=255)
              aTCGGCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGt                                                                                                                                 >  4:182371/1‑148 (MQ=255)
              aTCGGCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGt                                                                                                                                 >  4:324154/1‑148 (MQ=255)
              aTCGGCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCAAGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGGCATAATGACAAAGCATACTACCGAAGCGGt                                                                                                                                 >  3:178091/1‑148 (MQ=255)
                 ggCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAg                                                                                                                              <  2:329176/148‑1 (MQ=255)
                  gCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGt                                                                                                                            >  7:37827/1‑149 (MQ=255)
                   ccAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGtt                                                                                                                           <  1:286295/149‑1 (MQ=255)
                   ccAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGtt                                                                                                                           <  8:37827/149‑1 (MQ=255)
                   ccAGGTAAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCCCCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGtt                                                                                                                           <  5:264098/149‑1 (MQ=255)
                     aGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCa                                                                                                                         <  7:35039/149‑1 (MQ=255)
                       gTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGACCGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGt                                                                                                                       >  5:83034/1‑149 (MQ=255)
                           gCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGt                                                                                                                   <  4:125947/149‑1 (MQ=255)
                              aGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGc                                                                                                                <  3:55783/149‑1 (MQ=255)
                                gAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAg                                                                                                              >  2:7753/1‑149 (MQ=255)
                                 aaaGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCCCGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGt                                                                                                             >  2:238455/1‑149 (MQ=255)
                                      aCCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCgg                                                                                                        <  1:344072/149‑1 (MQ=255)
                                      aCCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGGGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCgg                                                                                                        <  7:303551/149‑1 (MQ=255)
                                       ccGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCggg                                                                                                       <  8:160364/149‑1 (MQ=255)
                                       ccGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCggg                                                                                                       <  8:185632/149‑1 (MQ=255)
                                       ccGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCggg                                                                                                       >  8:87897/1‑149 (MQ=255)
                                          ggTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTc                                                                                                    <  2:156089/149‑1 (MQ=255)
                                          ggTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTc                                                                                                    <  1:281793/149‑1 (MQ=255)
                                               aaaCCATCCGCCATATCTCCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTTCCGaa                                                                                                                                      <  6:265918/110‑1 (MQ=255)
                                                 aCCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATcc                                                                                             <  7:73948/149‑1 (MQ=255)
                                                 aCCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATcc                                                                                             >  2:177794/1‑149 (MQ=255)
                                                  ccATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGAACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATccc                                                                                            >  6:291304/1‑149 (MQ=255)
                                                     tCCGCCATATCCCCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGc                                                                                         <  4:230534/149‑1 (MQ=255)
                                                     tCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGc                                                                                         <  3:182371/149‑1 (MQ=255)
                                                       cGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCaa                                                                                       <  7:290614/149‑1 (MQ=255)
                                                        gCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCaa                                                                                       >  7:27964/1‑148 (MQ=255)
                                                           aaatCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGACCGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGACGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGGTCGCCATCCCTGCAACATc                                                                                   <  6:339455/147‑1 (MQ=255)
                                                            tatCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATc                                                                                   >  3:172720/1‑148 (MQ=255)
                                                            tatCACAACCGTCATCCCGCAGCGTGAAGCTACAGTAGAACGGTAACGGAACGTAAGGCCCCTTCGCCCCGCCAGAATGACCACGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCACTTTCGGGGTCGCGCTCCCGGCACCATCa                                                                                  >  4:245876/1‑149 (MQ=255)
                                                              tCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCaa                                                                                 >  5:139480/1‑148 (MQ=255)
                                                                     cGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGCCCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCg                                                                         >  2:66254/1‑149 (MQ=255)
                                                                      gTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGc                                                                        <  1:66254/149‑1 (MQ=255)
                                                                              gCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATccc                                                                 >  1:290289/1‑148 (MQ=255)
                                                                                      aaGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGgt                                                        <  2:297071/149‑1 (MQ=255)
                                                                                         cTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTCCCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGtt                                                     >  5:314924/1‑149 (MQ=255)
                                                                                            ccGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAg                                                  <  8:27964/149‑1 (MQ=255)
                                                                                                     ggTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCg                                         >  8:329215/1‑149 (MQ=255)
                                                                                                     ggTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCg                                         <  7:87897/149‑1 (MQ=255)
                                                                                                     ggTAACTGAACGTAAGGCCCATAACCCACGCCATAATGAACAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCg                                         <  2:135628/149‑1 (MQ=255)
                                                                                                      gTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCgg                                        <  3:324154/149‑1 (MQ=255)
                                                                                                         aCTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTa                                     <  5:154546/149‑1 (MQ=255)
                                                                                                           tGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTg                                                                                                                                <  4:184761/56‑1 (MQ=255)
                                                                                                           tGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTa                                     >  5:124333/1‑147 (MQ=255)
                                                                                                           tGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTa                                     >  1:221930/1‑147 (MQ=255)
                                                                                                             aaCGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTaaa                                   <  1:238455/147‑1 (MQ=255)
                                                                                                             aaCGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGAAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCAAATCCCCAAATTGGTGTTCAGGAAAAAGCGGGTAAAGt                                 <  1:267197/149‑1 (MQ=255)
                                                                                                                gTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTATCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCACATCCCCAAAATGGAGTTCAGGAAAAAGAGGGTAAAGTCGa                              >  8:238124/1‑149 (MQ=255)
                                                                                                                 tAAGGCCCATAACCAACGCCAAAATTACCAATCATACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGAc                             <  4:178091/149‑1 (MQ=255)
                                                                                                                 tAAGGCCCATAACACACGACATAATGACCAAGCATACTACCGAAGAGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTGCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGAc                             <  4:182599/149‑1 (MQ=255)
                                                                                                                   aGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGAcc                            <  2:5698/148‑1 (MQ=255)
                                                                                                                           aaCCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGa                   <  6:83034/149‑1 (MQ=255)
                                                                                                                              ccACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGAt                  >  1:56350/1‑147 (MQ=255)
                                                                                                                               cACGCCATTATGACCAAGCGTACTACCGAAGCGGTGATGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAg               <  6:314924/149‑1 (MQ=255)
                                                                                                                                  gCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGCCCATCCCCTGATAAGGGt            >  1:169146/1‑149 (MQ=255)
                                                                                                                                      tAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAACCCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCTCCTGATAAGGGTCGcc        >  3:321070/1‑149 (MQ=255)
                                                                                                                                      tAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCGcc        >  6:287779/1‑149 (MQ=255)
                                                                                                                                      tAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCGcc        >  4:154258/1‑149 (MQ=255)
                                                                                                                                      tAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCGcc        >  4:105716/1‑149 (MQ=255)
                                                                                                                                         tGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCGCCCTg     <  7:329215/149‑1 (MQ=255)
                                                                                                                                            ccAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCGCCCTGCaa  <  1:177794/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
ATTTACCCACAGACATCGGCCAGGTCAGCCAGGAAAGAACCGGGTAAAAACCATCCGCCATATCACCACCGTCATCACGCAGCGTGAAGCTACCGTAGAACGGTAACTGAACGTAAGGCCCATAACCCACGCCATAATGACCAAGCGTACTACCGAAGCGGTGAGGTTCAGTCCGTTGCAGTTTCGGGTTCGCCATCCCTGCAACATCAATAAAACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCGCCCTGCAA  >  NC_000913/2464455‑2464743

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑