Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,962,586 G→C 5.2% V69V (GTG→GTC dapF → diaminopimelate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,962,5860GC 5.2% 303.6 / 6.8 97V69V (GTG→GTCdapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base G (52/39);  new base C (2/3);  total (55/42)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

TGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGCGGCATCATCGTGCT  >  NC_002947/5962439‑5962735
                                                                                                                                                    |                                                                                                                                                     
tgggctgggcAACGACTTCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGa                                                                                                                                                     >  1:474889/1‑150 (MQ=255)
  ggctgggcAACGACTTCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACt                                                                                                                                                   >  1:713848/1‑150 (MQ=255)
    ctgggcAACGACTTCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTGTACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTc                                                                                                                                                 >  1:1056622/1‑150 (MQ=255)
     tgggcAACGACTTCATGGTCCTGGACCTGGTCAGACAGCACTCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTcc                                                                                                                                                <  1:1534104/150‑1 (MQ=255)
       ggcAACGACTTCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGCGGCGCCTAACCACCCGGAAGTGCACTTACGc                                                                                                                                              <  2:72313/150‑1 (MQ=255)
        gcAACGACTTCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTCCACTTCCggg                                                                                                                                             <  1:1219045/150‑3 (MQ=255)
            cGACTTCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGCTCGTCAACATGCAGAACAACCCGGAAGTGGACTTCCGCTAccc                                                                                                                                         >  2:16923/1‑149 (MQ=255)
             gACTTCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGc                                                                                                                                        >  2:1175107/1‑150 (MQ=255)
                ttCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGa                                                                                                                                                     >  1:706380/1‑134 (MQ=255)
                ttCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGa                                                                                                                                                     <  2:706380/134‑1 (MQ=255)
                ttCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATc                                                                                                                                     >  2:1702646/1‑150 (MQ=255)
                ttCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATc                                                                                                                                     <  2:564104/150‑1 (MQ=255)
                ttCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTCCACTTCCggg                                                                                                                                             >  2:1791434/1‑140 (MQ=255)
                ttCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTCCACTTCCggg                                                                                                                                             <  1:1791434/142‑3 (MQ=255)
                ttCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCAACCAAGCGGCAGATTCGACCAGTTGCTGATCGTCGACGCGACGAACAACCCGGAAGTTGACTTCCGCTACCGCATc                                                                                                                                     >  1:329443/1‑150 (MQ=255)
                      gTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAAc                                                                                                                               >  2:1796888/1‑150 (MQ=255)
                      gTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGAACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAAc                                                                                                                               >  2:1263437/1‑150 (MQ=255)
                       tCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGGATGGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACg                                                                                                                              >  2:548551/1‑150 (MQ=255)
                       tCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACg                                                                                                                              >  1:1139017/1‑150 (MQ=255)
                         cTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGcc                                                                                                                            >  2:861646/1‑150 (MQ=255)
                         cTGGACCTGGTCAGCCAGCACGCGCAAATCCCGAAGCAAAACCCATTGTTATAAGATGACGACCATACCGTCTTCGACATCGCTAAGTTTCCTATCGTCGATGCGCCGAACAACCCGGAAGTGGCCTTCCGATACCGCATCATCAACGcc                                                                                                                            >  1:1382661/1‑150 (MQ=255)
                          tGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCACTGGATGTGCGGCCGTACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCg                                                                                                                           >  2:168270/1‑150 (MQ=255)
                            gcccTTGTCAGCCAGACCGCGCACCTCCAGTCCTATCACGCCAAGCTATGTGGTGACTGCCATACCGGTGGCGGCTTCCACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGAc                                                                                                                         <  2:663906/148‑1 (MQ=255)
                             aCCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGGAGTGGACTTCCGCTACCGCATCTTCAACGCCGACg                                                                                                                        >  2:1159219/1‑150 (MQ=255)
                                 ggTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTc                                                                                                                    >  1:1521944/1‑150 (MQ=255)
                                 ggTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTc                                                                                                                    >  1:1241268/1‑150 (MQ=255)
                                  gTCNGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACGTCCGCTACCGCATCTTCAACGCCGACGGCTcc                                                                                                                   >  2:1193036/1‑150 (MQ=255)
                                    cagccagcACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGa                                                                                                                 <  2:478252/150‑1 (MQ=255)
                                      gccagcACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAgg                                                                                                               <  2:1007744/150‑1 (MQ=255)
                                       ccagcACGCGCACATCCAGCCCAAGCACGCCACGCTAGGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGt                                                                                                              >  1:1095917/1‑150 (MQ=255)
                                       ccagcACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTCCCGCATCTTCAACGCCGACGGCTCCGAggg                                                                                                              >  2:423733/1‑149 (MQ=255)
                                       ccagcACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGt                                                                                                              >  2:237970/1‑150 (MQ=255)
                                        cagcACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTc                                                                                                             <  1:1097632/150‑1 (MQ=255)
                                        cagcACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTc                                                                                                             <  2:1581557/150‑1 (MQ=255)
                                          gcACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGa                                                                                                           >  2:520093/1‑150 (MQ=255)
                                            acacgcACGTCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGc                                                                                                         >  2:1324523/4‑150 (MQ=255)
                                            aCGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATc                                                                                                                                     >  2:758600/1‑122 (MQ=255)
                                            aCGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATc                                                                                                                                     <  1:758600/122‑1 (MQ=255)
                                            aCGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGc                                                                                                         <  2:1758606/150‑1 (MQ=255)
                                              gcgcACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCaa                                                                                                       >  1:647769/1‑150 (MQ=255)
                                              gcgcACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCaa                                                                                                       >  1:671443/1‑150 (MQ=255)
                                                           ccAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTg                                                                                          >  2:1255872/1‑150 (MQ=255)
                                                           ccAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCCCTACCGCATCTTCAAc                                                                                                                               >  1:1083635/1‑113 (MQ=255)
                                                           ccAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCCCTACCGCATCTTCAAc                                                                                                                               <  2:1083635/113‑1 (MQ=255)
                                                            cAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAGCGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTgc                                                                                         >  1:781698/1‑150 (MQ=255)
                                                              aGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTAGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTgcgc                                                                                       <  1:16923/150‑1 (MQ=255)
                                                                cACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTgcgcgc                                                                                     <  1:1796888/150‑1 (MQ=255)
                                                                cACGCCAAGCAATGGGGTGACCGCAATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTgcgcgc                                                                                     >  2:1288235/1‑150 (MQ=255)
                                                                 aCGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCgct                                                                                    >  1:947735/1‑150 (MQ=255)
                                                                      aaGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCttt                                                                               <  2:647769/150‑1 (MQ=255)
                                                                          aaTGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGccc                                                                           >  2:1668694/1‑150 (MQ=255)
                                                                            tGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGc                                                                             >  1:1652522/1‑146 (MQ=255)
                                                                            tGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGc                                                                             <  2:1652522/146‑1 (MQ=255)
                                                                                     ccGCCATACCGGTGTCGGCTTCGACCAGTTGATGATCGTCGCGGCGCCTAACAACTCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCt                                                                 <  2:781698/150‑1 (MQ=255)
                                                                                      cGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGTTCCGAGGGCGAGCAATGCGGCAACGGTGCGCGCTGGTTTGCCCGCTTCGTGCTg                                                                <  2:465671/150‑1 (MQ=255)
                                                                                      cGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTg                                                                <  1:1263437/150‑1 (MQ=255)
                                                                                      cGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTg                                                                >  1:161660/1‑150 (MQ=255)
                                                                                        ccATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGGAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGa                                                              >  2:911582/1‑150 (MQ=255)
                                                                                          aTACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACa                                                            >  2:760950/1‑150 (MQ=255)
                                                                                          aTACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACa                                                            >  2:1457237/1‑150 (MQ=255)
                                                                                          aTACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTACGCCACCGCATCTTCAACGACGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACa                                                            >  2:1901169/1‑150 (MQ=255)
                                                                                           tACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGTGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACaa                                                           <  2:1336952/150‑1 (MQ=255)
                                                                                           tACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACaa                                                           <  1:1599486/150‑1 (MQ=255)
                                                                                            aCCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAg                                                          <  1:1175107/150‑1 (MQ=255)
                                                                                             ccGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCTGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAgc                                                         <  1:1796213/150‑1 (MQ=255)
                                                                                             ccGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAgc                                                         <  2:1750120/150‑1 (MQ=255)
                                                                                                   tCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGa                                                   <  1:861646/150‑1 (MQ=255)
                                                                                                     ggCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGGcc                                                 >  1:1615668/1‑150 (MQ=255)
                                                                                                      gCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCTGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCCATGCGGCAACGGTGCGCGCTTCTCTGCCCGCTTCGTGCTGGACAAGCGCCTGACcg                                                <  2:83359/150‑1 (MQ=255)
                                                                                                      gCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGATGTGGAATTCCGCTACCACATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACcg                                                >  1:976971/1‑150 (MQ=255)
                                                                                                      gCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGATGGACAAGCGCCTGACcg                                                >  2:1340789/1‑150 (MQ=255)
                                                                                                       cTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCACTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACcgc                                               <  2:1032517/150‑1 (MQ=255)
                                                                                                        ttCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACcgcg                                              >  1:1812449/1‑150 (MQ=255)
                                                                                                             ccAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGaaaaa                                         >  2:560158/1‑150 (MQ=255)
                                                                                                                  tGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCa                                    <  2:671443/150‑1 (MQ=255)
                                                                                                                       aTCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCcgc                               >  2:1534832/1‑150 (MQ=255)
                                                                                                                         cgtcgAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGt                             <  2:1521944/150‑1 (MQ=255)
                                                                                                                           tcgAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACTCCTACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAgc                                       <  2:536525/138‑1 (MQ=255)
                                                                                                                           tcgAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACTCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAgc                                       >  1:536525/1‑138 (MQ=255)
                                                                                                                           tcgAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTgg                           >  1:636444/1‑150 (MQ=255)
                                                                                                                             gAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGaa                         >  1:538283/1‑150 (MQ=255)
                                                                                                                             gAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGaa                         >  1:751883/1‑150 (MQ=255)
                                                                                                                               ggCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAAc                       <  1:584826/150‑1 (MQ=255)
                                                                                                                                  gcCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCaa                    >  2:683561/1‑150 (MQ=255)
                                                                                                                                     gAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAagag                 <  2:1241268/150‑1 (MQ=255)
                                                                                                                                      aacaacCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGt                             >  2:827712/1‑137 (MQ=255)
                                                                                                                                      aacaacCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGt                             <  1:827712/137‑1 (MQ=255)
                                                                                                                                      aacaacCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGc                <  1:309828/150‑1 (MQ=255)
                                                                                                                                       acaacCCGGAAGGGGACTTCCGCTACCGCATCTTCGACGCCGACGGCTCCGATGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGGCAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGATCg               <  2:858152/150‑1 (MQ=255)
                                                                                                                                         aacCCGGAAGTGGACTTCCGCTACCGCATCTTCAATGGCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGCGGc             <  1:1442855/150‑1 (MQ=255)
                                                                                                                                         aacCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGCGGc             <  1:548551/150‑1 (MQ=255)
                                                                                                                                           cccGGAAGTGGACTTCCGCTACCACATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGCGGcat           >  1:378295/1‑150 (MQ=255)
                                                                                                                                           agcGGAAGTCCACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGt                             >  1:1600069/3‑132 (MQ=255)
                                                                                                                                           agcGGAAGTCCACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGt                             <  2:1600069/130‑1 (MQ=255)
                                                                                                                                              ggAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGCGGcatcat        >  2:734208/1‑150 (MQ=255)
                                                                                                                                               gAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGCGGCatcatc       <  2:538283/150‑1 (MQ=255)
                                                                                                                                                 aGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGCGGCATCATCGt     <  1:683561/150‑1 (MQ=255)
                                                                                                                                                    ggACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGCGGCATCATCGTGCt  <  1:560158/150‑1 (MQ=255)
                                                                                                                                                    ggACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGCGGCATCATCGTGCt  <  1:237970/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTCAGCCAGCACGCGCACATCCAGCCCAAGCACGCCAAGCAATGGGGTGACCGCCATACCGGTGTCGGCTTCGACCAGTTGCTGATCGTCGAGGCGCCGAACAACCCGGAAGTGGACTTCCGCTACCGCATCTTCAACGCCGACGGCTCCGAGGTCGAGCAATGCGGCAACGGTGCGCGCTGCTTTGCCCGCTTCGTGCTGGACAAGCGCCTGACCGCGAAAAAGCGCATCCGCGTGGAAACCAAGAGCGGCATCATCGTGCT  >  NC_002947/5962439‑5962735

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: