Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,248 G→A 7.1% A429A (GCC→GCT dnaK ← chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,2480GA 7.1% 596.3 / 24.1 197A429A (GCC→GCTdnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base G (75/108);  new base A (4/10);  total (79/118)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 200 total aligned reads displayed.

TCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGC  >  NC_002947/5378099‑5378388
                                                                                                                                                     |                                                                                                                                            
tCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGtcgtcg                                                                                                                                              <  2:5240399/150‑1 (MQ=255)
      tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                                                                                                        <  1:3534648/150‑1 (MQ=255)
      tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                                                                                                        >  2:4015637/1‑150 (MQ=255)
        gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                                                                                                                      >  2:1192294/1‑150 (MQ=255)
             tCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAaca                                                                                                                                 >  1:4489858/1‑150 (MQ=255)
             tCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAaca                                                                                                                                 >  2:1312848/1‑150 (MQ=255)
              cACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAacac                                                                                                                                >  1:4694574/1‑150 (MQ=255)
               aCTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACAcc                                                                                                                               >  1:687812/1‑150 (MQ=255)
                 ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                                                                                                                             >  1:3686955/1‑150 (MQ=255)
                 ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                                                                                                                             >  1:3357684/1‑150 (MQ=255)
                      tttGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACGTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCGCGGGTGGTCGTCGGCAGTCGAGAACACCTGCGAct                                                                                                                        >  2:4335744/1‑150 (MQ=255)
                        tGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtc                                                                                                                      <  2:514602/150‑1 (MQ=255)
                          cGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtctt                                                                                                                    <  2:1611958/150‑1 (MQ=255)
                           ggCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTg                                                                                                                   >  1:398183/1‑150 (MQ=255)
                           ggCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTg                                                                                                                   <  2:2245278/150‑1 (MQ=255)
                               cgccgcGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCg                                                                                                               <  2:4770822/150‑1 (MQ=255)
                                gccgcGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCgg                                                                                                              >  2:518669/1‑150 (MQ=255)
                                  cgcgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGa                                                                                                            >  2:1013088/1‑150 (MQ=255)
                                    cgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtg                                                                                                          <  1:2412622/150‑1 (MQ=255)
                                       ggTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtg                                                                                                       <  2:1458800/150‑1 (MQ=255)
                                       ggTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtg                                                                                                       >  2:3763025/1‑150 (MQ=255)
                                       ggTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACCGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtg                                                                                                       <  1:559191/150‑1 (MQ=255)
                                          gCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGgtggtg                                                                                                    <  1:130684/150‑1 (MQ=255)
                                            tGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtgg                                                                                                         >  2:2092826/1‑143 (MQ=255)
                                            tGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGtt                                                                                                  >  2:2702662/1‑150 (MQ=255)
                                             ggCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGttc                                                                                                 <  1:5030180/150‑1 (MQ=255)
                                             ggCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGttc                                                                                                 <  2:687812/150‑1 (MQ=255)
                                              gCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCgg                                                                                                                                             <  2:2907571/105‑1 (MQ=255)
                                              gCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGttgc                                                                                                <  1:1080889/150‑3 (MQ=255)
                                                 gAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTtcttct                                                                                             <  1:4766858/150‑1 (MQ=255)
                                                  aaTGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTtctc                                                                                            <  1:1227248/150‑1 (MQ=255)
                                                  aaTGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGttgcgttc                                                                                            <  1:2767954/150‑7 (MQ=255)
                                                       cAGCCAGGTCGTACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGACACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCa                                                                                       <  2:1591362/150‑1 (MQ=255)
                                                        aacgaGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAg                                                                                      >  1:5295366/5‑150 (MQ=255)
                                                        aacgaGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAg                                                                                      >  2:5370768/5‑150 (MQ=255)
                                                            aGGTCGAACTTGCCCAGCGACTTGTTCTGCGCATCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGcgcg                                                                                  <  2:3373106/150‑1 (MQ=255)
                                                            aGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGcgcg                                                                                  <  2:3037164/150‑1 (MQ=255)
                                                              gTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGt                                                                                <  1:1022832/150‑1 (MQ=255)
                                                              gTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGt                                                                                <  2:3935071/150‑1 (MQ=255)
                                                               tCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGtc                                                                               <  1:2199384/150‑1 (MQ=255)
                                                                 gAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGtcat                                                                             >  1:5343834/1‑150 (MQ=255)
                                                                       gCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCAcgcc                                                                       >  2:4652727/1‑150 (MQ=255)
                                                                         ccAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCAcgccgc                                                                     <  2:1062923/150‑1 (MQ=255)
                                                                          cAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACgccgcc                                                                    >  1:231207/1‑150 (MQ=255)
                                                                          cAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACgccgcc                                                                    >  1:3534776/1‑150 (MQ=255)
                                                                           aGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAGGGTGGTGTTCTGCTCGATCAGCGCGGTCATCACGCCGccc                                                                   >  2:4843129/1‑150 (MQ=255)
                                                                            gCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtg                                                                                                       <  2:2232947/113‑1 (MQ=255)
                                                                            gCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTtcttct                                                                                             <  2:1385741/123‑1 (MQ=255)
                                                                             cGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAt                                                                 >  1:3839807/1‑150 (MQ=255)
                                                                               aCTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATgg                                                               >  1:3611950/1‑150 (MQ=255)
                                                                               aCTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATgg                                                               >  1:1783789/1‑150 (MQ=255)
                                                                                     tCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGa                                                         <  1:3947745/150‑1 (MQ=255)
                                                                                      cTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGAt                                                        <  2:955358/150‑1 (MQ=255)
                                                                                       tGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATa                                                       >  1:4455218/1‑150 (MQ=255)
                                                                                        gcgcAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATAc                                                      >  1:5200629/1‑150 (MQ=255)
                                                                                            aGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAg                                                  <  2:636455/150‑1 (MQ=255)
                                                                                             gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                 <  2:4267956/150‑1 (MQ=255)
                                                                                             gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                 <  1:221667/150‑1 (MQ=255)
                                                                                             gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                 <  1:4078634/150‑1 (MQ=255)
                                                                                             gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                 <  2:1085434/150‑1 (MQ=255)
                                                                                             gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                 <  1:1326279/150‑1 (MQ=255)
                                                                                                tGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTc                                              >  1:1689438/1‑150 (MQ=255)
                                                                                                tGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTc                                              >  1:2596485/1‑150 (MQ=255)
                                                                                                   ttACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCATTCGAGAACACCTGCGACTTCTTGGTCGGGATggtggt                                                                                                     >  1:2974917/1‑92 (MQ=255)
                                                                                                   ttACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATggtggt                                                                                                     <  2:2974917/92‑1 (MQ=255)
                                                                                                   ttACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGc                                           <  2:2843011/150‑1 (MQ=255)
                                                                                                     aCGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGTCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCgg                                         >  2:3200358/1‑150 (MQ=255)
                                                                                                       gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                       >  1:2071799/1‑150 (MQ=255)
                                                                                                       gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                       <  2:442630/150‑1 (MQ=255)
                                                                                                           aCCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGAc                                   <  1:1192294/150‑1 (MQ=255)
                                                                                                           aCCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGAc                                   <  2:4088311/150‑1 (MQ=255)
                                                                                                            ccTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACg                                  >  2:34370/1‑150 (MQ=255)
                                                                                                            ccTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACg                                  >  2:2546242/1‑150 (MQ=255)
                                                                                                              ttGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTc                                <  2:4489759/150‑1 (MQ=255)
                                                                                                              ttGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTc                                >  2:3542578/1‑150 (MQ=255)
                                                                                                              ttGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTc                                >  1:5269734/1‑150 (MQ=255)
                                                                                                                gcagcaCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGACAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCca                              <  2:2767857/150‑1 (MQ=255)
                                                                                                                 cagcaCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcag                             <  2:4869438/150‑1 (MQ=255)
                                                                                                                    caCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcag                             >  1:3923805/1‑147 (MQ=255)
                                                                                                                    caCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcag                             <  2:3923805/147‑1 (MQ=255)
                                                                                                                     aCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagc                         >  1:1574073/1‑150 (MQ=255)
                                                                                                                       gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagcag                       <  2:98162/150‑1 (MQ=255)
                                                                                                                       gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagcag                       <  1:4431275/150‑1 (MQ=255)
                                                                                                                       gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagcag                       <  2:3326904/150‑1 (MQ=255)
                                                                                                                       gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagcag                       <  1:3675398/150‑1 (MQ=255)
                                                                                                                        tGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAgg                      <  2:1783789/150‑1 (MQ=255)
                                                                                                                          aaTGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGAc                    <  2:3248575/150‑1 (MQ=255)
                                                                                                                            gggTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGAGACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGt                  <  1:4941927/149‑1 (MQ=255)
                                                                                                                             cgTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTc                 <  1:4848448/149‑1 (MQ=255)
                                                                                                                              gTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCTGCTGCTGGACGTca                <  2:5448859/150‑2 (MQ=255)
                                                                                                                                cACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCttt              <  1:4462899/150‑1 (MQ=255)
                                                                                                                                 aCGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTc             >  1:4561315/1‑150 (MQ=255)
                                                                                                                                   ggccgactGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTCCTTGAc           <  1:4512200/144‑1 (MQ=255)
                                                                                                                                   ggCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCAc           <  2:3839807/150‑1 (MQ=255)
                                                                                                                                   ggCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCAc           <  2:100934/150‑1 (MQ=255)
                                                                                                                                            gTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc  >  1:3606300/1‑150 (MQ=255)
                                                                                                                                            gTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc  <  2:4694574/150‑1 (MQ=255)
                                                                                                                                            gTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc  <  1:3090941/150‑1 (MQ=255)
                                                                                                                                            gTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc  <  2:5269734/150‑1 (MQ=255)
                                                                                                                                            gTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGc  <  1:2829555/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
TCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGC  >  NC_002947/5378099‑5378388

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