Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 127,121 C→G 8.6% intergenic (‑31/+39) thrB ← / ← PP_0122 homoserine kinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947127,1210CG 8.6% 899.3 / 53.3 276intergenic (‑31/+39)thrB/PP_0122homoserine kinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (137/115);  new base G (19/5);  total (156/120)
Fisher's exact test for biased strand distribution p-value = 2.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

AAGAAATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGTTTT  >  NC_002947/126972‑127269
                                                                                                                                                     |                                                                                                                                                    
aagaaATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGTAGACTGACATGATGAAAAATTGCCCATACGGGCACTTccgc                                                                                                                                                      >  2:5127521/1‑150 (MQ=255)
aagaaATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCAGTTccgc                                                                                                                                                      >  2:5127718/1‑150 (MQ=255)
  gaaATTGCTGTTCTCGGTACCGGCGGCAATGCCCCGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTccgccg                                                                                                                                                    >  2:5910204/1‑150 (MQ=255)
    aaTTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCggg                                                                                                                                                  <  2:5168284/150‑1 (MQ=255)
    aaTTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCggg                                                                                                                                                  <  1:6148330/150‑1 (MQ=255)
    aaTTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCgga                                                                                                                                                  <  2:3224142/150‑2 (MQ=255)
       tGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAg                                                                                                                                               <  1:4024401/150‑1 (MQ=255)
          tGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGc                                                                                                                                            >  1:1778790/1‑150 (MQ=255)
                 ggTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTg                                                                                                                                     <  1:1853426/150‑1 (MQ=255)
                   tACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAt                                                                                                                                   <  2:5249684/150‑1 (MQ=255)
                     ccGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAtgt                                                                                                                                 <  1:4593935/150‑1 (MQ=255)
                          ggcAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACg                                                                                                                            >  1:8612971/1‑150 (MQ=255)
                           gcAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACgg                                                                                                                           >  1:8775306/1‑150 (MQ=255)
                             caggCCCGGGAAGTCGAGCAGACGGCCCAGCTCGTCCGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTg                                                                                                                         <  1:7107213/147‑1 (MQ=255)
                             aaTGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTg                                                                                                                         >  2:7782600/1‑150 (MQ=255)
                                    tGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCa                                                                                                                  <  2:1676599/150‑1 (MQ=255)
                                    aggAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCa                                                                                                                  <  2:1227452/149‑1 (MQ=255)
                                          tCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTAc                                                                                                            >  1:6795399/1‑150 (MQ=255)
                                          tCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTAc                                                                                                            <  2:7582388/150‑1 (MQ=255)
                                          tCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTAc                                                                                                            <  2:8621584/150‑1 (MQ=255)
                                                  aCGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCTCGGGGGTGAAGAGTGTCATGATGATAAATTGCCCATCCGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTGGaa                                                                                                    <  2:2333768/150‑1 (MQ=255)
                                                  aCGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGaa                                                                                                    <  2:7752921/150‑1 (MQ=255)
                                                    ggCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGa                                                                                                  >  1:5469910/1‑150 (MQ=255)
                                                      cccAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATc                                                                                                <  2:3106473/150‑1 (MQ=255)
                                                      cccAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATc                                                                                                <  2:4863359/150‑1 (MQ=255)
                                                       ccAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCt                                                                                               >  1:7779288/1‑150 (MQ=255)
                                                         aGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTc                                                                                             <  2:7062610/150‑1 (MQ=255)
                                                           cTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCa                                                                                           >  1:9089070/1‑150 (MQ=255)
                                                            tCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGCAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTGCCACTCGCAGTTCTTCCAg                                                                                          >  2:7012227/1‑150 (MQ=255)
                                                            tCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAg                                                                                          <  1:156110/150‑1 (MQ=255)
                                                                   gcgcCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGaa                                                                                   <  2:6797228/150‑1 (MQ=255)
                                                                     gcCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATc                                                                                 <  1:7293534/150‑1 (MQ=255)
                                                                       cAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAg                                                                               >  2:1041069/1‑150 (MQ=255)
                                                                        aGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGc                                                                              >  1:2700268/1‑150 (MQ=255)
                                                                        aGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGc                                                                              >  2:854639/1‑150 (MQ=255)
                                                                         gAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCa                                                                             >  2:874690/1‑150 (MQ=255)
                                                                          aaaGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCAt                                                                            >  2:5828952/1‑150 (MQ=255)
                                                                          aaaGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCAt                                                                            >  1:3507977/1‑150 (MQ=255)
                                                                          aaaGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCAt                                                                            >  1:7271583/1‑150 (MQ=255)
                                                                               tttCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCg                                                                       >  2:2322134/1‑150 (MQ=255)
                                                                               tttCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCg                                                                       <  1:1240696/150‑1 (MQ=255)
                                                                                 tCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGc                                                                     <  1:4266857/150‑1 (MQ=255)
                                                                                 tCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGc                                                                     <  1:3787561/150‑1 (MQ=255)
                                                                                  ccggCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCt                                                                    >  1:66490/4‑150 (MQ=255)
                                                                                  ccAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCt                                                                    >  2:2840758/1‑150 (MQ=255)
                                                                                   cAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTg                                                                   >  2:7994083/1‑150 (MQ=255)
                                                                                   cAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTg                                                                   <  1:5797856/150‑1 (MQ=255)
                                                                                   cAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACCGGCACTTCCCGGCGGAAGTGcccg                                                                                                                                         <  1:5378574/80‑3 (MQ=37)
                                                                                         aGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGa                                                             >  1:8636983/1‑150 (MQ=255)
                                                                                               ggTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAAt                                                       >  2:2609452/1‑150 (MQ=255)
                                                                                               ggTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAAt                                                       >  1:4825655/1‑150 (MQ=255)
                                                                                                gTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATg                                                      <  1:435902/150‑1 (MQ=255)
                                                                                                     gggggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAtt                                                 <  1:6420355/150‑1 (MQ=255)
                                                                                                     gggggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAtt                                                 <  2:1482195/150‑1 (MQ=255)
                                                                                                       gggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTgg                                               >  2:2193767/1‑150 (MQ=255)
                                                                                                        ggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGc                                              >  2:819187/1‑150 (MQ=255)
                                                                                                        ggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGc                                              >  2:5487564/1‑150 (MQ=255)
                                                                                                          tGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCAt                                            >  2:8250359/1‑150 (MQ=255)
                                                                                                          tGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCAt                                            <  1:1190658/150‑1 (MQ=255)
                                                                                                             aGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAg                                         <  2:1058278/150‑1 (MQ=255)
                                                                                                              gACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGt                                        <  2:5702922/150‑1 (MQ=255)
                                                                                                              gACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGt                                        <  2:6253285/150‑1 (MQ=255)
                                                                                                                 tGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGcc                                     >  1:7506842/1‑150 (MQ=255)
                                                                                                                  gACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCa                                    >  2:4802026/1‑150 (MQ=255)
                                                                                                                   aCATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCAt                                   >  1:7997194/1‑150 (MQ=255)
                                                                                                                   aCATGATGAAAAATTGCCCATACGGGCACTACCGCCGGTAAGTGCCGGTTTGATGTGAACGGTGAGCGTAAGCTTAACACACGAAGATCTTCCAGGACGGAATCATCATGTCCGGCTGGTCCGATCGAATGAAATTTGCATAAGAGCCAt                                   >  1:3382232/1‑150 (MQ=255)
                                                                                                                     atgatgAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCa                                 <  1:6997305/150‑1 (MQ=255)
                                                                                                                      cgatgaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAg                                <  1:5999642/149‑1 (MQ=255)
                                                                                                                       gatgaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAgc                               <  1:2156753/150‑1 (MQ=255)
                                                                                                                        atgaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAgcg                              >  1:8995273/1‑150 (MQ=255)
                                                                                                                        atgaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAgcg                              >  2:1364145/1‑150 (MQ=255)
                                                                                                                          gaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGcgcg                            >  2:2419606/1‑150 (MQ=255)
                                                                                                                           aaaaaTTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGt                           <  1:1689915/150‑1 (MQ=255)
                                                                                                                            aaaaTTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTa                          <  2:3486132/150‑1 (MQ=255)
                                                                                                                            aaaaTTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTa                          <  1:2325945/150‑1 (MQ=255)
                                                                                                                               aTTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCa                       <  1:1431482/150‑1 (MQ=255)
                                                                                                                                   cccATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGaaa                   >  1:2799396/1‑150 (MQ=255)
                                                                                                                                    ccATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGc                                                                              <  2:4712277/90‑1 (MQ=255)
                                                                                                                                    ccATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGc                                                                              >  1:4712277/1‑90 (MQ=255)
                                                                                                                                    ccATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGaaaa                  <  2:4311776/150‑1 (MQ=255)
                                                                                                                                      aTACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATa                <  1:7660924/150‑1 (MQ=255)
                                                                                                                                      aTACGGGCACTTCCGCCGGGAAGTGCCGGGATGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATa                <  2:9119405/150‑1 (MQ=255)
                                                                                                                                         ccggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg             >  2:4088013/3‑150 (MQ=255)
                                                                                                                                         cGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg             >  2:7413895/1‑150 (MQ=255)
                                                                                                                                         acggCACTTCCCGCCGGAAGACCCGGGTTGATGAGAACGGTGAGCGACAGCTTACCAATCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGATGGTAGGATCGAATCAAATTTGCATAAGTGCCATAAGCGCGTACCAGAAAATAAgg             >  2:6647242/3‑150 (MQ=255)
                                                                                                                                          gggCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGc            >  2:6456481/1‑150 (MQ=255)
                                                                                                                                          gggCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGc            >  1:5639468/1‑150 (MQ=255)
                                                                                                                                          cggCACTTCCCGGCGGAAGTGCCCGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGc            >  1:6335796/2‑150 (MQ=255)
                                                                                                                                           ggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCt           >  1:3922906/1‑150 (MQ=255)
                                                                                                                                           ggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCt           >  2:2526160/1‑150 (MQ=255)
                                                                                                                                           ggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCt           >  2:3320562/1‑150 (MQ=255)
                                                                                                                                            gCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCtt          >  1:2472514/1‑150 (MQ=255)
                                                                                                                                            gCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCtt          >  2:2699942/1‑150 (MQ=255)
                                                                                                                                             cACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTg         <  1:8197629/150‑1 (MQ=255)
                                                                                                                                              aCTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGc        <  2:4642212/150‑1 (MQ=255)
                                                                                                                                                  ccgccgGGAAGTGCCGGGTTGATGTGATCGGTGCGCGTCAGCTTACCTCTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGtt    <  2:1778790/150‑1 (MQ=255)
                                                                                                                                                  ccgccgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGtt    >  1:3313886/1‑150 (MQ=255)
                                                                                                                                                  ccgccgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGtt    >  1:7489178/1‑150 (MQ=255)
                                                                                                                                                   cgccgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGttt   <  1:4914853/150‑1 (MQ=255)
                                                                                                                                                    gccgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGtttt  <  1:5910204/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGAAATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGTTTT  >  NC_002947/126972‑127269

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: