Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,571,422 A→T 6.3% intergenic (+16/+22) PP_3150 → / ← sad‑II hypothetical protein/NAD+‑dependent succinate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,571,4220AT 6.3% 510.4 / 15.6 153intergenic (+16/+22)PP_3150/sad‑IIhypothetical protein/NAD+‑dependent succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (51/92);  new base T (6/4);  total (57/96)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

CACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAGTAGC  >  NC_002947/3571273‑3571571
                                                                                                                                                     |                                                                                                                                                     
     gTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAagca                                                                                                                                                  <  1:1144617/150‑1 (MQ=255)
     gTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAagca                                                                                                                                                  >  1:646151/1‑150 (MQ=255)
     gTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAagca                                                                                                                                                  <  2:2079252/150‑1 (MQ=255)
     gTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAagca                                                                                                                                                  <  2:306700/150‑1 (MQ=255)
      tGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAagcag                                                                                                                                                 >  1:2332891/1‑150 (MQ=255)
       ggCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAagcagc                                                                                                                                                <  1:4498336/150‑1 (MQ=255)
       ggCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAagcagc                                                                                                                                                <  1:3048097/150‑1 (MQ=255)
        gCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGcc                                                                                                                                               >  1:2659101/1‑150 (MQ=255)
          gggTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGcccc                                                                                                                                             >  2:2750285/1‑150 (MQ=255)
          gggTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGcccc                                                                                                                                             >  2:902939/1‑150 (MQ=255)
            gTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCTGCTTTGCGGCCCCAg                                                                                                                                           <  2:4012358/150‑1 (MQ=255)
            gTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCTGCTTTGCGGCCCCAg                                                                                                                                           <  2:4012273/150‑1 (MQ=255)
             tAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAgg                                                                                                                                          <  1:392079/150‑1 (MQ=255)
             tAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAgg                                                                                                                                          <  2:1983981/150‑1 (MQ=255)
             tAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCCGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAgg                                                                                                                                          <  1:2733127/150‑1 (MQ=255)
             tAACCGCAATGAACGCGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAgg                                                                                                                                          <  2:472828/150‑1 (MQ=255)
                ccGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTc                                                                                                                                       <  2:380464/150‑1 (MQ=255)
                ccGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTc                                                                                                                                       <  2:787563/150‑1 (MQ=255)
                 cGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCTGCTTTGCGGCCCCAgctgct                                                                                                                                      <  1:3894837/150‑6 (MQ=255)
                  gCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGTGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCAt                                                                                                                                     <  2:2354697/150‑1 (MQ=255)
                  gCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCAt                                                                                                                                     <  2:2980528/150‑1 (MQ=255)
                    aaTGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGa                                                                                                                                   <  2:2076274/150‑1 (MQ=255)
                    aaTGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGa                                                                                                                                   <  1:3302737/150‑1 (MQ=255)
                      tGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATCAAt                                                                                                                                 <  2:765077/150‑1 (MQ=255)
                       gAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATc                                                                                                                                <  1:5063521/150‑1 (MQ=255)
                        aaCGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCa                                                                                                                               >  1:3148164/1‑150 (MQ=255)
                            ggCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGTCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTa                                                                                                                           <  2:4374982/150‑1 (MQ=255)
                                cgTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGt                                                                                                                       >  2:2254079/1‑150 (MQ=255)
                                 gTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGACTCAGTTACGGTc                                                                                                                      >  1:4676995/1‑150 (MQ=255)
                                   aTGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCa                                                                                                                    <  2:2685227/150‑1 (MQ=255)
                                    tGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAg                                                                                                                   <  1:2152060/150‑1 (MQ=255)
                                    tGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAg                                                                                                                   >  1:4939844/1‑150 (MQ=255)
                                       tGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGccc                                                                                                                >  2:5045434/1‑149 (MQ=255)
                                         gCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCcaca                                                                                                              >  2:2370797/1‑150 (MQ=255)
                                          cGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCacac                                                                                                             <  2:2209522/150‑1 (MQ=255)
                                           ggggTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACg                                                                                                            >  2:672261/1‑150 (MQ=255)
                                               ttGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCt                                                                                                        <  1:4851777/150‑1 (MQ=255)
                                                tGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTg                                                                                                       <  2:2659101/150‑1 (MQ=255)
                                                 gTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTgg                                                                                                      <  1:2006760/150‑1 (MQ=255)
                                                 gTATACCGCCTGTTCGCGGGCGCGGCATGAAGTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTgg                                                                                                      <  2:982739/150‑1 (MQ=255)
                                                    taCCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCa                                                                                                   <  1:5099252/150‑1 (MQ=255)
                                                     aCCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGCATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCAt                                                                                                  <  2:2135186/150‑1 (MQ=255)
                                                      ccGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCAtt                                                                                                 >  2:4298999/1‑150 (MQ=255)
                                                      ccGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCAtt                                                                                                 >  2:2900800/1‑150 (MQ=255)
                                                      ccGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCAtt                                                                                                 >  1:892249/1‑150 (MQ=255)
                                                      ccGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGGTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGACTCAGTTACGGTCCAGCCACACGGTCTGGGCAtt                                                                                                 >  1:3826250/1‑150 (MQ=255)
                                                       cGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTg                                                                                                <  1:1777018/150‑1 (MQ=255)
                                                        gCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGt                                                                                              <  2:651498/150‑1 (MQ=255)
                                                         gcTGTGCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGt                                                                                              <  2:3426800/149‑1 (MQ=255)
                                                         ccTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGTGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACTGTCCAGCCACACGGTCTGGGCATTGGt                                                                                              <  1:3038243/150‑1 (MQ=255)
                                                         ccTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGt                                                                                              <  1:2605398/150‑1 (MQ=255)
                                                            gTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGaa                                                                                           <  2:1556382/150‑1 (MQ=255)
                                                            gTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGaa                                                                                           <  1:3582978/150‑1 (MQ=255)
                                                            gTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGaa                                                                                           >  2:697289/1‑150 (MQ=255)
                                                            gTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGaa                                                                                           <  2:1196936/150‑1 (MQ=255)
                                                                 cgGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTcgc                                                                                      <  1:4958330/150‑1 (MQ=255)
                                                                 cgGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTcgc                                                                                      <  1:4958515/150‑1 (MQ=255)
                                                                      gcgGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGcaca                                                                                 <  1:2138623/150‑1 (MQ=255)
                                                                      gcgGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGcaca                                                                                 <  1:669768/150‑1 (MQ=255)
                                                                      gcgGCATGAATAGATCGTGCCGGGGGGGATGCTGGACTGGGTCAACGACCACGTACCGGATTGAAGCAGCTCGCGCCGCAAAGCAGCCCCAGGGACATGACTCAGGTACGGTCCAGCCACACGGACTGGGCATTGGTGAACTCGCGcaca                                                                                 >  2:1715823/1‑150 (MQ=255)
                                                                           aTGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCaaa                                                                            <  1:664312/150‑1 (MQ=255)
                                                                           aTGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCaaa                                                                            <  1:3145662/150‑1 (MQ=255)
                                                                                       tGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTc                                                                >  2:1302340/1‑150 (MQ=255)
                                                                                       tGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTc                                                                <  1:1657530/150‑1 (MQ=255)
                                                                                         ccGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCAc                                                              <  2:3716556/150‑1 (MQ=255)
                                                                                          cGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACg                                                             >  1:3520491/1‑150 (MQ=255)
                                                                                                 gATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGaaa                                                      <  1:1459230/150‑1 (MQ=255)
                                                                                                  aTGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAAc                                                     <  1:2984896/150‑1 (MQ=255)
                                                                                                    gCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCa                                                   >  2:3170139/1‑150 (MQ=255)
                                                                                                    gCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCa                                                   >  2:2425529/1‑150 (MQ=255)
                                                                                                       ggACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCActc                                                <  1:1041053/150‑1 (MQ=255)
                                                                                                       ggACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCActc                                                <  2:665885/150‑1 (MQ=255)
                                                                                                            gggTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTcttctt                                           <  2:3161754/150‑1 (MQ=255)
                                                                                                              gTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCa                                         <  2:1942499/150‑1 (MQ=255)
                                                                                                               tCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCAc                                        <  2:2332891/150‑1 (MQ=255)
                                                                                                               tCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCAc                                        <  1:2750285/150‑1 (MQ=255)
                                                                                                                     aCCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGccacc                                  <  1:4041238/150‑1 (MQ=255)
                                                                                                                         ggTTCGGGATTGAAGCATCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAAc                              <  1:23864/150‑1 (MQ=255)
                                                                                                                           ttCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGc                            <  1:2296481/150‑1 (MQ=255)
                                                                                                                                    gAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCgggg                   <  1:3946232/150‑1 (MQ=255)
                                                                                                                                    gAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCgggg                   <  1:3946286/150‑1 (MQ=255)
                                                                                                                                     aaGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCggggg                  >  1:3348958/1‑150 (MQ=255)
                                                                                                                                     aaGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCggggg                  >  2:4479460/1‑150 (MQ=255)
                                                                                                                                     aaGCAGCTGCGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCggggg                  >  2:585344/1‑150 (MQ=255)
                                                                                                                                       gcagcTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTgg                                                                                                      >  2:1780178/1‑64 (MQ=255)
                                                                                                                                       gcagcTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTgg                                                                                                      <  1:1780178/64‑1 (MQ=255)
                                                                                                                                       gcagcTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTc                >  2:2521513/1‑150 (MQ=255)
                                                                                                                                       gcagcTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTc                <  1:297981/150‑1 (MQ=255)
                                                                                                                                       gcagcTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTc                >  1:3456856/1‑150 (MQ=255)
                                                                                                                                       gcagcTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTc                <  1:4737641/150‑1 (MQ=255)
                                                                                                                                        acccTGGGGCTGCTTTGCGGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTcc               >  1:4314418/4‑150 (MQ=255)
                                                                                                                                          ccTGGGGCTGCTTTGCGGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGa             >  1:2386228/2‑150 (MQ=255)
                                                                                                                                           cTGGGGCTGCTTTGCGGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAGCTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGaa            >  1:3744594/1‑150 (MQ=255)
                                                                                                                                           cTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGaa            <  2:2428406/150‑1 (MQ=255)
                                                                                                                                            tGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAg           <  1:902939/150‑1 (MQ=255)
                                                                                                                                                gCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCa       <  1:1237270/150‑1 (MQ=255)
                                                                                                                                                 ccGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAg      >  2:4955570/1‑150 (MQ=255)
                                                                                                                                                  cGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCGTTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAGt     >  2:637384/1‑150 (MQ=255)
                                                                                                                                                  cGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAGt     <  1:2936437/150‑1 (MQ=255)
                                                                                                                                                   gCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAGTa    >  1:1609764/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAGCAGCCCCAGGGTCATGAATCAGTTACGGTCCAGCCACACGGTCTGGGCATTGGTGAACTCGCGCACACCAAAGTGCGACAGCTCACGGCCGAAACCACTCTTCTTCACGCCACCAAACGCCACGCGGGGGTCCGAAGCGCAGTAGC  >  NC_002947/3571273‑3571571

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