Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,430,392 G→C 6.2% S189S (TCG→TCC)  PP_3045 → ATP‑dependent protease ClpP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,430,3920GC 6.2% 258.1 / 7.1 81S189S (TCG→TCC) PP_3045ATP‑dependent protease ClpP
Reads supporting (aligned to +/- strand):  ref base G (34/42);  new base C (3/2);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

GGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGACC  >  NC_002947/3430246‑3430539
                                                                                                                                                  |                                                                                                                                                   
ggTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAg                                                                                                                                                  <  2:966735/150‑1 (MQ=255)
     tGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGt                                                                                                                                             <  2:1195964/150‑1 (MQ=255)
       ccGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTgg                                                                                                                                           <  2:1061861/150‑1 (MQ=255)
             aCCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgg                                                                                                                                     <  1:446431/150‑1 (MQ=255)
             aCCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgg                                                                                                                                     <  1:677716/150‑1 (MQ=255)
              ccTTCGTCGGGTGGCGGACGTGCTGGATCAGGCATTCGAGGCAATGATCGCAGCCCATAAGGCCACGGCTCCGAACAGCGATGACGCGGAGCTGAGACGTCTGGTCAAGGCCGAAAGTTGGCTCACGGACTCGGAGGCCGTGGCACTggg                                                                                                                                    >  1:144338/1‑150 (MQ=255)
                    tcgGGTGTCGGCCGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCCCGGCCTCGGAGGCCGTGGCGCTGGGGCTGGc                                                                                                                              >  2:134094/1‑150 (MQ=255)
                              gACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGa                                                                                                                                                   >  2:236426/1‑119 (MQ=255)
                              gACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGa                                                                                                                                                   <  1:236426/119‑1 (MQ=255)
                               aCGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCg                                                                                                                   <  2:880648/150‑1 (MQ=255)
                                 gTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGt                                                                                                                 <  1:458838/150‑1 (MQ=255)
                                   gCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCCGAGGCCGt                                                                                                                                             >  1:1103532/1‑120 (MQ=255)
                                   gCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCCGAGGCCGt                                                                                                                                             <  2:1103532/120‑1 (MQ=255)
                                     ttGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                             <  2:144338/150‑1 (MQ=255)
                                       gATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGgt                                                                                                           <  2:1129255/150‑1 (MQ=255)
                                       gATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGgt                                                                                                           <  1:1072284/150‑1 (MQ=255)
                                       gATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGgt                                                                                                           >  1:912015/1‑150 (MQ=255)
                                           aGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                             <  1:1339046/144‑1 (MQ=255)
                                           aGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                             >  2:1339046/1‑144 (MQ=255)
                                               cTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGt                                                                                                   <  2:548249/150‑1 (MQ=255)
                                                ttCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTTGCCGATGAGGTCGGTAATGGTGTCAAGGtt                                                                                                  >  1:83840/1‑150 (MQ=255)
                                                ttCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGtt                                                                                                  >  1:997509/1‑150 (MQ=255)
                                                ttCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGtt                                                                                                  >  1:933033/1‑150 (MQ=255)
                                                 tCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGt                                                                                                                                             >  2:1206625/1‑106 (MQ=255)
                                                 tCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGt                                                                                                                                             <  1:1206625/106‑1 (MQ=255)
                                                    aGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAgg                                                                                              >  1:572589/1‑150 (MQ=255)
                                                         aTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGc                                                                                         >  2:336775/1‑150 (MQ=255)
                                                         aTGATCGCAGCCTATAAGGCCAAGGCTCCGAACAGCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCGCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGc                                                                                         >  1:1182689/1‑150 (MQ=255)
                                                          tGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGcc                                                                                        <  1:336775/150‑1 (MQ=255)
                                                             tCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTg                                                                                                                                            >  2:1253636/1‑95 (MQ=255)
                                                             tCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTg                                                                                                                                            <  1:1253636/95‑1 (MQ=255)
                                                             tCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTg                                                                                     >  2:1325478/1‑150 (MQ=255)
                                                                  gCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCaa                                                                                <  2:758776/150‑1 (MQ=255)
                                                                   ccTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAg                                                                               <  2:572589/150‑1 (MQ=255)
                                                                     tataAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                             <  2:670688/118‑1 (MQ=255)
                                                                     tataAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                             >  1:670688/1‑118 (MQ=255)
                                                                     tataAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGgc                                                                             >  2:772486/1‑150 (MQ=255)
                                                                         aGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCgcgg                                                                         <  2:997509/150‑1 (MQ=255)
                                                                           gCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCgcggcg                                                                       >  2:1402737/1‑150 (MQ=255)
                                                                             cAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGAt                                                                     <  2:1099702/150‑1 (MQ=255)
                                                                              aaGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGgc                                                                                                                                          <  2:1408691/80‑1 (MQ=255)
                                                                              aaGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGgc                                                                                                                                          >  1:1408691/1‑80 (MQ=255)
                                                                              aaGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATg                                                                    >  1:396346/1‑150 (MQ=255)
                                                                                ggCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGgc                                                                  <  2:1137389/150‑1 (MQ=255)
                                                                                  cTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGgcgc                                                                <  2:912015/150‑1 (MQ=255)
                                                                                   tCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGAt                                                                     >  2:66584/1‑144 (MQ=255)
                                                                                   tCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGAt                                                                     <  1:66584/144‑1 (MQ=255)
                                                                                    ccGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGt                                                              <  1:569394/150‑1 (MQ=255)
                                                                                      gAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAg                                                                                                     <  1:441134/109‑1 (MQ=255)
                                                                                      gAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAg                                                                                                     >  2:441134/1‑109 (MQ=255)
                                                                                                   gcgGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGgcgc                                                                >  1:89447/1‑133 (MQ=255)
                                                                                                   gcgGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGgcgc                                                                <  2:89447/133‑1 (MQ=255)
                                                                                                      gAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCaa                                            <  1:134094/150‑1 (MQ=255)
                                                                                                           gagaCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCt                                       >  2:260999/1‑150 (MQ=255)
                                                                                                                ttCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGAAAATATGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTgg                                  <  1:865658/149‑1 (MQ=255)
                                                                                                                gTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGCCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTgg                                  >  1:373719/1‑150 (MQ=255)
                                                                                                                 tCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCACGGCGCGGCGATGGCGCGTTACCCCAATTCCCCGCAAGTCCTATTGGa                                 <  2:414761/150‑1 (MQ=255)
                                                                                                                 tCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                                    <  1:383807/147‑1 (MQ=255)
                                                                                                                 tCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                                    >  2:383807/1‑147 (MQ=255)
                                                                                                                  cTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGATGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCTTTGCGAGAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGAt                                <  1:1106299/150‑1 (MQ=255)
                                                                                                                   tGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTGAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATc                               >  1:1157043/1‑150 (MQ=255)
                                                                                                                   tGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATc                               >  2:323088/1‑150 (MQ=255)
                                                                                                                   tGGTCAATGCCGAAAGTTGGCTCACGGCCACGGAGCCCGTGGCTCGGGGAGTGGCCTGGGAGGTCGGGTTTGGTCTCAAGAAACCGCGCCGCTGTTGTTTAGGCGCGGCGATGGCGCGTTACCGCAATACGCCGCAAGGCCTGTTTGATc                               >  1:678645/1‑150 (MQ=255)
                                                                                                                     gTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAg                             >  2:1176490/1‑150 (MQ=255)
                                                                                                                     gTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGACGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAg                             >  1:481566/1‑150 (MQ=255)
                                                                                                                      tCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGc                            <  1:203726/150‑1 (MQ=255)
                                                                                                                      tCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGc                            <  1:772486/150‑1 (MQ=255)
                                                                                                                                     ggCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGCCAGTAATGGTGTCAAGGTTCAGGCTTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAggg             >  1:603510/1‑150 (MQ=255)
                                                                                                                                       cTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGAt           >  2:573037/1‑150 (MQ=255)
                                                                                                                                        tCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATg          <  2:550735/150‑1 (MQ=255)
                                                                                                                                        tCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATg          >  1:550735/1‑150 (MQ=255)
                                                                                                                                        ccACGGCCTCCGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTAcc                                                          <  2:81648/101‑1 (MQ=255)
                                                                                                                                        ccACGGCCTCCGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTAcc                                                          >  1:81648/2‑102 (MQ=255)
                                                                                                                                         cACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATgg         <  1:323088/150‑1 (MQ=255)
                                                                                                                                         cACGGCCTCCGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATgg         >  1:1404700/1‑150 (MQ=255)
                                                                                                                                          aCGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATggg        <  1:1402737/150‑1 (MQ=255)
                                                                                                                                            ggCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGt                                                                                           >  2:1302844/1‑65 (MQ=255)
                                                                                                                                            ggCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGt                                                                                           <  1:1302844/65‑1 (MQ=255)
                                                                                                                                             gCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTg     <  1:1176490/150‑1 (MQ=255)
                                                                                                                                               cTCGGAGGCCGTGGCGCTGGGGCTGGCCGCTGCCGTGCTCAAGGTTTTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGAc   <  2:1182689/150‑1 (MQ=255)
                                                                                                                                                tCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGAcc  <  2:396346/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGACC  >  NC_002947/3430246‑3430539

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