Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,497,184 T→A 6.7% M110K (ATG→AAG)  opgD → OPG biosynthetic periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,497,1840TA 6.7% 159.9 / 2.4 60M110K (ATG→AAG) opgDOPG biosynthetic periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base A (2/2);  total (31/29)
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.63e-01

GTGGCGTGGTGCGCCTCGTCCGTTACCTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGC  >  NC_000913/1497038‑1497331
                                                                                                                                                  |                                                                                                                                                   
gtggcgtggTGCGCCTCGTCCGTTACCTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATatgggaatgg                                                                                                                                                   <  7:16262/149‑1 (MQ=255)
   gcgtggTGCGCCTCGTCCGTTACCTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCT‑AGTTCTTCCATATGGGAATGGGAtt                                                                                                                                               <  1:287924/149‑1 (MQ=255)
      tggtgCGCCTCGTCCGTTACCTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTcc                                                                                                                                             <  7:124288/149‑1 (MQ=255)
       ggtgCGCCTCGTCCGTTACCTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCg                                                                                                                                            >  1:187637/1‑149 (MQ=255)
            gcCTCGTCCGTTACCTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTcgcc                                                                                                                                       <  8:210312/149‑1 (MQ=255)
               tCGTCCGTTACCTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCcgcg                                                                                                                                    >  5:147273/1‑149 (MQ=255)
                  tCCGTTACCTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTc                                                                                                                                 <  2:160943/149‑1 (MQ=255)
                      ttACCTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTAt                                                                                                                             <  1:138137/149‑1 (MQ=255)
                          cTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGttt                                                                                                                         >  3:239859/1‑149 (MQ=255)
                          cTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCG‑ATGtttt                                                                                                                        >  1:279212/1‑149 (MQ=255)
                              cacGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTg                                                                                                                     <  5:34093/149‑1 (MQ=255)
                              cacGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTg                                                                                                                     <  5:350230/149‑1 (MQ=255)
                              cacGCTGGAGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTg                                                                                                                     <  4:60465/149‑1 (MQ=255)
                                 gCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAg                                                                                                                  <  2:187637/149‑1 (MQ=255)
                                     gCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATcc                                                                                                              <  4:229039/149‑1 (MQ=255)
                                      cGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCTGCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATc                                                                                                               >  8:173602/1‑147 (MQ=255)
                                          aaTGACGCCGCAGGCTTATAACAGTATTCAATACGATGCCGAAAAATTGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCaa                                                                                                         >  2:278552/1‑149 (MQ=255)
                                          aaTGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCaa                                                                                                         >  4:162897/1‑149 (MQ=255)
                                          aaTGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCaa                                                                                                         >  4:115672/1‑149 (MQ=255)
                                                gccgcAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATc                                                                                                   >  2:282639/1‑149 (MQ=255)
                                                   gcAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTgg                                                                                                <  5:174209/149‑1 (MQ=255)
                                                       gCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGcgcg                                                                                            >  4:49723/1‑149 (MQ=255)
                                                         ttATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGt                                                                                           >  3:366849/1‑148 (MQ=255)
                                                          tataACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAAAGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTg                                                                                          >  8:163862/1‑148 (MQ=255)
                                                          tataACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAAAGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTg                                                                                          >  8:169348/1‑148 (MQ=255)
                                                            taACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCG‑ATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGaaa                                                                                       <  2:279212/148‑1 (MQ=255)
                                                            taACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGATGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGACGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACAAATCTGGCGCGTGaaa                                                                                       >  1:293836/1‑149 (MQ=255)
                                                                aGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAAAGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCa                                                                                   <  7:169348/149‑1 (MQ=255)
                                                                aGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAAAGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCa                                                                                   <  7:163862/149‑1 (MQ=255)
                                                                 gTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGc                                                                                                           <  2:81776/124‑1 (MQ=255)
                                                                 gTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGc                                                                                                           >  1:81776/1‑124 (MQ=255)
                                                                  tATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGGTCGTATGTTTTCTGTAGCTCCAGCAACACATCTGGCGCGTGAAATTCACt                                                                                 >  4:164149/1‑149 (MQ=255)
                                                                  tATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTAAACTGGACGCTCAGTTCTTCCATATGGCAATGGTATTCCGACCCCCCGTTCGTATGTTTTCTGAAGAACCCGCAACACATGTGGCGCGTGAAAATCACt                                                                                 >  1:298347/1‑149 (MQ=255)
                                                                     tCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTc                                                                              >  5:39065/1‑149 (MQ=255)
                                                                      cAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTc                                                                              <  2:101196/148‑1 (MQ=255)
                                                                        aTACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGcc                                                                           >  3:173999/1‑149 (MQ=255)
                                                                         tACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGccc                                                                          >  4:76636/1‑149 (MQ=255)
                                                                             acgCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAAtt                                                                                     >  2:228465/1‑134 (MQ=255)
                                                                             acgCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAAtt                                                                                     <  1:228465/134‑1 (MQ=255)
                                                                                ccGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAgttg                                                                   >  8:32411/1‑149 (MQ=255)
                                                                                 cGAAAAATCGCTCTGGCATATCTTTGCTAACAGTCAACTTGACGCTCAGTTCTTCCATATGTGAATGGGATTACGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGgcgc                                                                                             <  6:94660/122‑1 (MQ=255)
                                                                                 cGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCATCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGgcgc                                                                                             >  5:94660/1‑122 (MQ=255)
                                                                                      aaTCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCaaa                                                             <  3:52889/149‑1 (MQ=255)
                                                                                           ctctGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACaa                                                        <  4:239859/149‑1 (MQ=255)
                                                                                                gCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGa                                                     <  3:76636/147‑1 (MQ=255)
                                                                                                gCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATg                                                   >  1:236547/1‑149 (MQ=255)
                                                                                                   tAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAg                                                <  5:44956/149‑1 (MQ=255)
                                                                                                   tAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAg                                                <  7:6925/149‑1 (MQ=255)
                                                                                                   tAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAg                                                <  4:82800/149‑1 (MQ=255)
                                                                                                        ttGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGtt                                           <  1:282639/149‑1 (MQ=255)
                                                                                                            gaACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCACATACAACGATGCAGGTGTTGATa                                       >  2:99731/1‑149 (MQ=255)
                                                                                                                     aCTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAAt                              <  2:236547/149‑1 (MQ=255)
                                                                                                                      cTGGTCGCTATGTTGTTCGATATGGGAATGTGATTACGTCGCCTCGTTCGTATGTTTTCTGTAGTTCCAGCAACACATCTGCCGCGTGAAATTTACTTTCGCCCGGAGTTGTTCAAATACATCGATGCAGGTGTTTATACataaaa                                <  2:293836/146‑3 (MQ=255)
                                                                                                                      cTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAAtt                             <  3:62718/149‑1 (MQ=255)
                                                                                                                              tCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTc                                                                              <  6:130757/92‑1 (MQ=255)
                                                                                                                              tCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTc                                                                              >  5:130757/1‑92 (MQ=255)
                                                                                                                              tCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGGTTTCTGTAGATCCAGCAACCCATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAACCAATTAGAAggg                      >  6:411/1‑148 (MQ=255)
                                                                                                                              tCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCCTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACATTCGCCCGGAGTTGTTCAAATACACCGATGCACGTGTTGATACAAAACAATTAGAAggg                      >  6:291478/1‑148 (MQ=255)
                                                                                                                              tCAGTTATTCCATATGGGAATGTGCTTCCGGCGCCGCGTTCGTATGTTTTCTGGAGATACAGCAAAACATCTGGCGCGTGAAATTCACTTTCGCCAGGAGTTGGTAAAATACAACGATGCAGGTGTTGATACAAAACAATTAGCAggg                      >  5:369349/1‑148 (MQ=255)
                                                                                                                                     ttcCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGa              >  1:68819/1‑149 (MQ=255)
                                                                                                                                      tcCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGAt             >  1:367529/1‑149 (MQ=255)
                                                                                                                                            tgggaatgggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGCTCCAGCAACACATCTGGCGCGTGAAATTCACTTTAGCCCGGAGTTGTTCAAATACAACGATGCAGGAGTTGATACAAAACAATTAGAACGGCAAAGCGATCTCGGc       >  2:124231/1‑149 (MQ=255)
                                                                                                                                             gggaatgggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCt      >  2:20721/1‑149 (MQ=255)
                                                                                                                                                 atgggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAAGGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGc  >  5:162287/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTGGCGTGGTGCGCCTCGTCCGTTACCTGACACGCTGGCGACAATGACGCCGCAGGCTTATAACAGTATTCAATACGACGCCGAAAAATCGCTCTGGCATAACGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGC  >  NC_000913/1497038‑1497331

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: ‑