Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,430,392 G→C 5.6% 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 5.6% 410.5 / 7.7 128S189S (TCG→TCCPP_3045ATP‑dependent protease ClpP
Reads supporting (aligned to +/- strand):  ref base G (57/63);  new base C (2/6);  total (59/69)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

GCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGACCGA  >  NC_002947/3430244‑3430541
                                                                                                                                                    |                                                                                                                                                     
gCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCgg                                                                                                                                                      >  1:832363/1‑150 (MQ=255)
     gATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGcc                                                                                                                                                 <  2:1805106/150‑1 (MQ=255)
        gCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTg                                                                                                                                              <  1:570598/150‑1 (MQ=255)
        gCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTg                                                                                                                                              <  1:860941/150‑1 (MQ=255)
         ccGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCGGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCCGAGGCCGTga                                                                                                                                             <  1:1156079/150‑2 (MQ=255)
         ccGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTgg                                                                                                                                             <  1:135542/150‑1 (MQ=255)
         ccGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTgg                                                                                                                                             <  2:1057678/150‑1 (MQ=255)
          cGAGGACCTTCGTCGGGTGGCGGCCGTGCTTGATCAGGCCTTCGATTCAATGATCGCAGCCTATAAGGGCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCGCACGGCCTCCGAGGCCGTgag                                                                                                                                            <  1:515470/150‑3 (MQ=255)
           gAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGgcg                                                                                                                                           >  2:1915783/1‑150 (MQ=255)
           gAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAACGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCCGAGGCCgcgagc                                                                                                                                           <  1:2105164/150‑6 (MQ=255)
            aGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGgcgc                                                                                                                                          >  2:2306219/1‑150 (MQ=255)
               aCCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgg                                                                                                                                       <  1:175710/150‑1 (MQ=255)
                 cTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgggg                                                                                                                                     >  1:2278053/1‑150 (MQ=255)
                 cTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgggg                                                                                                                                     <  1:2379688/150‑1 (MQ=255)
                 cTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgggg                                                                                                                                     >  1:482077/1‑150 (MQ=255)
                 cTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgggg                                                                                                                                     >  2:1766636/1‑150 (MQ=255)
                 cTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCAGCCAGGTGGGTCAGGTCCTGGGAGGGCGGTGCGCTgggg                                                                                                                                     >  2:263747/1‑150 (MQ=255)
                  ttCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTg                                                                                                                                              >  2:157630/1‑140 (MQ=255)
                  ttCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTg                                                                                                                                              <  1:157630/140‑1 (MQ=255)
                        gggTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCg                                                                                                                              >  1:105111/1‑150 (MQ=255)
                          gTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATCAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGAt                                                                                                                            <  2:1612453/150‑1 (MQ=255)
                            ggCGGACGTGCTTGATCAGGCCTCCGAGGCAATGATCGCAGCCTATCCCATCGGTGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGa                                                                                                                          <  1:1353181/150‑1 (MQ=255)
                               ggACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATAGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGt                                                                                                                       <  1:301003/150‑1 (MQ=255)
                                gACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTc                                                                                                                      >  2:625946/1‑150 (MQ=255)
                                     gCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTaa                                                                                                                 >  2:1996223/1‑150 (MQ=255)
                                       ttGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                               >  2:1578644/1‑150 (MQ=255)
                                       ttGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                               >  2:2183994/1‑150 (MQ=255)
                                       ttGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                               <  2:105111/150‑1 (MQ=255)
                                            cAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTc                                                                                                          >  2:2378882/1‑150 (MQ=255)
                                               gCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTAAAg                                                                                                       >  2:1933354/1‑150 (MQ=255)
                                               gCCGTCGATGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAg                                                                                                       >  2:1242827/1‑150 (MQ=255)
                                                  ttCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGtt                                                                                                    >  1:270617/1‑150 (MQ=255)
                                                     gAGGGAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAg                                                                                                 <  1:1242827/150‑1 (MQ=255)
                                                       ggCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGgt                                                                                                             >  2:1472478/1‑136 (MQ=255)
                                                       ggCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGgt                                                                                                             <  1:1472478/136‑1 (MQ=255)
                                                       ggCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGtt                                                                                                    <  1:1348209/145‑1 (MQ=255)
                                                         cAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCgg                                                                                                                    <  1:1186180/127‑1 (MQ=255)
                                                         cAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCgg                                                                                                                    >  2:1186180/1‑127 (MQ=255)
                                                           aTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAGTGGTGTAAAGTTTCTGTCGTGc                                                                                           >  2:1342971/1‑150 (MQ=255)
                                                           aTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGc                                                                                           >  1:1317424/1‑150 (MQ=255)
                                                             gATCGCAGCCTATAAGGCCAAGGCTCCGAACATAGATGACGCGCAGCTCAGACATCTGGTCAATGCCGAAAGTACCCACACGGCCTCGGACGCCGTCGCGCCGGGGCCGGCCGAAGAGGTCGGGAAGGGTGCCAAGGTCGCGGCGTGccg                                                                                         >  2:2023932/1‑149 (MQ=255)
                                                                cGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTgg                                                                                      <  1:625946/150‑1 (MQ=255)
                                                                     ccTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAg                                                                                                 <  2:1326970/134‑1 (MQ=255)
                                                                     ccTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAg                                                                                                 >  1:1326970/1‑134 (MQ=255)
                                                                      cTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAgg                                                                                <  1:2116766/150‑1 (MQ=255)
                                                                      cTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAgg                                                                                <  1:263747/150‑1 (MQ=255)
                                                                            ggCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAg                                                                                                                         <  1:141300/103‑1 (MQ=255)
                                                                            ggCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAg                                                                                                                         >  2:141300/1‑103 (MQ=255)
                                                                             gCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCgcggcg                                                                         >  2:1587976/1‑150 (MQ=255)
                                                                              ccAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGa                                                                        <  1:1725395/150‑1 (MQ=255)
                                                                               cAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAAGCGTGCCTTGGGCAAGGCGCGGCGag                                                                       >  1:362500/1‑149 (MQ=255)
                                                                               cAAGGCGCCGAACAGCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGAt                                                                       >  2:1677046/1‑150 (MQ=255)
                                                                                  ggCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTGCATGCCGAAAGTTTGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGgc                                                                    <  1:38059/150‑1 (MQ=255)
                                                                                    cTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGgcgc                                                                  <  2:470829/150‑1 (MQ=255)
                                                                                     tCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGACGGcgcg                                                                 >  1:2291568/1‑150 (MQ=255)
                                                                                       cGAACATCGATGACGCGGAGCTGAGACGTCTGGTCATTGGCGAACGTTGGCGCCCGGGCTGGTGGTCCGGTGGGCGGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGtt                                                               >  1:657350/1‑150 (MQ=255)
                                                                                            aTCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTc                                                                                                                      >  1:829559/1‑90 (MQ=255)
                                                                                            aTCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTc                                                                                                                      <  2:829559/90‑1 (MQ=255)
                                                                                            aTCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTg                                                                                            >  1:2235058/1‑116 (MQ=255)
                                                                                                     gcgGAGCTCAGACGCCTGGTCAATGCCGCAAGGTGGCTCCCGGCCTCCGCGGCCGTGGCGCTGGGGCTCGCCGTTTAGGGCGGGAATGGTGCCTAGGTTCAGGCGTGCCTTGGGCAGAGCGCGGCCGTGGGGCGGTACCGGAATACCCCg                                                 >  2:1907743/1‑150 (MQ=255)
                                                                                                       ggAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTaa                                                                                                                 <  1:1439496/84‑1 (MQ=255)
                                                                                                        gAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCaa                                              <  1:1578644/150‑1 (MQ=255)
                                                                                                           cTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGcc                                           <  1:1677046/150‑1 (MQ=255)
                                                                                                             gagaCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCt                                         <  2:832363/150‑1 (MQ=255)
                                                                                                                aCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGTAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                                      <  2:1363744/150‑1 (MQ=255)
                                                                                                                aCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                                      <  1:1888911/150‑1 (MQ=255)
                                                                                                                aCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                                      <  2:1317424/150‑1 (MQ=255)
                                                                                                                aCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                                      <  1:691125/150‑1 (MQ=255)
                                                                                                                 cGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTg                                     <  1:933431/150‑1 (MQ=255)
                                                                                                                  gTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCt                                                                                                                                         >  1:1847719/1‑49 (MQ=255)
                                                                                                                  gTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCt                                                                                                                                         <  2:1847719/49‑1 (MQ=255)
                                                                                                                  gTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAg                                                                                                                         <  2:1799031/65‑1 (MQ=255)
                                                                                                                  gTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAg                                                                                                                         >  1:1799031/1‑65 (MQ=255)
                                                                                                                  gTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGtt                                                                                                    >  2:1366897/1‑86 (MQ=255)
                                                                                                                  gTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATAcc                                                    >  2:1725211/1‑134 (MQ=255)
                                                                                                                  gTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATAcc                                                    <  1:1725211/134‑1 (MQ=255)
                                                                                                                  gTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCaa                                              <  2:417022/140‑1 (MQ=255)
                                                                                                                  gTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCaa                                              >  1:417022/1‑140 (MQ=255)
                                                                                                                     tGGTCAATGCCGAAAGTTGGCTCACGGTCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATc                                 <  2:563489/150‑1 (MQ=255)
                                                                                                                     tGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATg                                                                      >  2:1343457/1‑113 (MQ=255)
                                                                                                                      ggTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCa                                >  2:69178/1‑150 (MQ=255)
                                                                                                                          aaTGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGc                                                                                           <  1:1316080/87‑1 (MQ=255)
                                                                                                                          aaTGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGc                                                >  1:678620/1‑130 (MQ=255)
                                                                                                                          aaTGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGc                                                <  2:678620/130‑1 (MQ=255)
                                                                                                                          aaTGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTc                            >  2:1724571/1‑150 (MQ=255)
                                                                                                                             gCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAg                         <  1:1587976/150‑1 (MQ=255)
                                                                                                                                gAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATg                                                                      >  2:897071/1‑102 (MQ=255)
                                                                                                                                gAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGcc                      >  2:1155108/1‑150 (MQ=255)
                                                                                                                                        gCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGGCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGa              <  2:270617/150‑1 (MQ=255)
                                                                                                                                          tCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGGGTGCCTTGGGAAAGGCGCGGCGATCGCACCTTACCCAAATACCGTCCAAGCTCTATGTTATGTGGTCCAGTCCAGTAAGGATg            >  2:266473/1‑150 (MQ=255)
                                                                                                                                          tCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGt                                                                                             <  1:585814/69‑1 (MQ=255)
                                                                                                                                          tCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGt                                                                                             >  2:585814/1‑69 (MQ=255)
                                                                                                                                            aCGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATggg          >  2:301798/1‑150 (MQ=255)
                                                                                                                                            aCGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATggg          >  1:2175571/1‑150 (MQ=255)
                                                                                                                                                ccTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                                      <  2:2252400/118‑1 (MQ=255)
                                                                                                                                                ccTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                                      >  1:2252400/1‑118 (MQ=255)
                                                                                                                                                 cTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGAc     <  2:1607750/150‑1 (MQ=255)
                                                                                                                                                  tCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGAt                                  <  2:672919/120‑1 (MQ=255)
                                                                                                                                                  tCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGAt                                  >  1:672919/1‑120 (MQ=255)
                                                                                                                                                    ggAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGActta  <  2:2278053/150‑4 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGACCGA  >  NC_002947/3430244‑3430541

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