Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,808,132 G→T 5.5% K47N (AAG→AAT) ‡ kdsA‑I → 2‑dehydro‑3‑deoxyphosphooctonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,808,1320GT 5.5% 625.2 / 18.3 182K47N (AAG→AAT) ‡kdsA‑I2‑dehydro‑3‑deoxyphosphooctonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (76/96);  new base T (6/4);  total (82/100)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 184 total aligned reads displayed.

GCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCA  >  NC_002947/1807983‑1808281
                                                                                                                                                     |                                                                                                                                                     
gCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAg                                                                                                                                                       >  2:3321545/1‑150 (MQ=255)
 ccTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGc                                                                                                                                                      >  1:3911401/1‑150 (MQ=255)
  cTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCt                                                                                                                                                     <  2:1809712/150‑1 (MQ=255)
   tGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTc                                                                                                                                                    <  1:4100827/150‑1 (MQ=255)
   tGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTc                                                                                                                                                    >  1:640751/1‑150 (MQ=255)
   tGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTc                                                                                                                                                    <  1:2330317/150‑1 (MQ=255)
    gATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCg                                                                                                                                                   <  2:1078687/150‑1 (MQ=255)
    gATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCg                                                                                                                                                   <  2:1262855/150‑1 (MQ=255)
    gATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCg                                                                                                                                                   <  1:2080590/150‑1 (MQ=255)
    gATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCg                                                                                                                                                   >  1:1951857/1‑150 (MQ=255)
     aTCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCgg                                                                                                                                                  >  2:855153/1‑150 (MQ=255)
     aTCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCgg                                                                                                                                                  >  1:3773299/1‑150 (MQ=255)
     aTCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCgg                                                                                                                                                  <  2:1442107/150‑1 (MQ=255)
       cGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTa                                                                                                                                                <  1:543778/150‑1 (MQ=255)
         aTGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATc                                                                                                                                              <  1:337442/150‑1 (MQ=255)
           gACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGATACCGAGCTTCTCGGTAAccc                                                                                                                                            <  2:1581992/150‑1 (MQ=255)
              ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGCAGCGCGGTCGCCCGtg                                                                                                                                         >  1:4181158/1‑149 (MQ=255)
              ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTa                                                                                                                                         >  1:4181191/1‑150 (MQ=255)
              ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTa                                                                                                                                         >  1:3269464/1‑150 (MQ=255)
              ccAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGATACCGAGCTTCTCGGTAACCCgca                                                                                                                                         <  1:1673321/150‑3 (MQ=255)
               cAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTAc                                                                                                                                        >  1:2496386/1‑150 (MQ=255)
               cAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTAc                                                                                                                                        <  1:3359679/150‑1 (MQ=255)
               cAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTAc                                                                                                                                        <  1:806617/150‑1 (MQ=255)
                   agaTCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACgtgt                                                                                                                                    >  1:257703/1‑150 (MQ=255)
                       cATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCATTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCaa                                                                                                                                <  1:2724315/150‑1 (MQ=255)
                        aTTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGCTTCTCGGTAACCCGCACgtactcttc                                                                                                                               <  2:258641/150‑8 (MQ=38)
                          tCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGc                                                                                                                             >  2:1390008/1‑150 (MQ=255)
                              gTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGt                                                                                                                         >  2:3048469/1‑150 (MQ=255)
                                cGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGttt                                                                                                                       >  2:3442514/1‑150 (MQ=255)
                                cGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGttt                                                                                                                       <  2:1125464/150‑1 (MQ=255)
                                 ggTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTc                                                                                                                      <  2:4224117/150‑1 (MQ=255)
                                 ggTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTc                                                                                                                      <  2:4292606/150‑1 (MQ=255)
                                 ggTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTc                                                                                                                      <  2:3167452/150‑1 (MQ=255)
                                 ggTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTc                                                                                                                      <  1:2142731/150‑1 (MQ=255)
                                       aTCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAg                                                                                                                >  1:3193557/1‑150 (MQ=255)
                                       aTCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAg                                                                                                                >  1:4398141/1‑150 (MQ=255)
                                          gagaTCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGcc                                                                                                             >  1:1080795/1‑150 (MQ=255)
                                             aTCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAAc                                                                                                          <  2:3927933/150‑1 (MQ=255)
                                                gCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGt                                                                                                       >  2:3100546/1‑150 (MQ=255)
                                                 ccAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGtt                                                                                                      >  1:2596561/1‑150 (MQ=255)
                                                  cAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGTATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTtc                                                                                                     <  2:3355578/150‑1 (MQ=255)
                                                  cAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTtc                                                                                                     >  2:3396658/1‑150 (MQ=255)
                                                  cAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTtc                                                                                                     >  1:739887/1‑150 (MQ=255)
                                                      gACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTtcgtcg                                                                                                 >  2:784709/1‑150 (MQ=255)
                                                         aaGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTa                                                                                              >  1:1589948/1‑150 (MQ=255)
                                                         aaGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTa                                                                                              <  2:4237653/150‑1 (MQ=255)
                                                           gCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGACCAACCGTTCGTCGGTAAc                                                                                            >  2:1732100/1‑150 (MQ=255)
                                                              gTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTc                                                                                         <  1:1255105/150‑1 (MQ=255)
                                                               ttCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTcc                                                                                        <  1:3048643/150‑1 (MQ=255)
                                                                  gTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTAc                                                                                     <  1:2728751/150‑1 (MQ=255)
                                                                      tGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTg                                                                                 >  1:3678198/1‑150 (MQ=255)
                                                                       gTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTgg                                                                                <  2:231167/150‑1 (MQ=255)
                                                                        ttCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGc                                                                               >  2:1228405/1‑150 (MQ=255)
                                                                         tCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGcc                                                                              <  1:4345586/150‑1 (MQ=255)
                                                                           ggcggcATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCt                                                                            >  1:266341/1‑150 (MQ=255)
                                                                              ggcATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGc                                                                         <  2:2713567/150‑1 (MQ=255)
                                                                                cATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCAt                                                                       <  1:1947703/150‑1 (MQ=255)
                                                                                   gAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGa                                                                    >  2:459175/1‑150 (MQ=255)
                                                                                   gAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGa                                                                    >  2:4050818/1‑150 (MQ=255)
                                                                                    aaCGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGaa                                                                   <  1:1117492/150‑1 (MQ=255)
                                                                                        ttCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGc                                                                                                                             <  2:4498357/88‑1 (MQ=255)
                                                                                        ttCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGc                                                                                                                             >  1:4498357/1‑88 (MQ=255)
                                                                                        ttCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAaggac                                                               <  2:4273904/150‑4 (MQ=255)
                                                                                         tCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAggg                                                              <  1:3633386/150‑1 (MQ=255)
                                                                                         tCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAggg                                                              <  2:1953049/150‑1 (MQ=255)
                                                                                                      gACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTc                                                 <  2:3776888/150‑1 (MQ=255)
                                                                                                      gACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTc                                                 <  1:1098580/150‑1 (MQ=255)
                                                                                                      gACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCAGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTc                                                 <  1:2807703/150‑1 (MQ=255)
                                                                                                         cTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                                              <  2:2272046/150‑1 (MQ=255)
                                                                                                         cTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                                              <  2:1700104/150‑1 (MQ=255)
                                                                                                         cTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGaa                                              <  1:2205670/150‑1 (MQ=255)
                                                                                                          tGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAag                                             >  1:930443/1‑150 (MQ=255)
                                                                                                          tGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAag                                             >  1:1443552/1‑150 (MQ=255)
                                                                                                          tGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTGCGAag                                             >  1:2242244/1‑150 (MQ=255)
                                                                                                               aTGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATc                                        <  2:640751/150‑1 (MQ=255)
                                                                                                                tGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCa                                       <  2:1635803/150‑1 (MQ=255)
                                                                                                                 gAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCaa                                      <  1:1390008/150‑1 (MQ=255)
                                                                                                                   aGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAgc                                    <  2:2036244/150‑1 (MQ=255)
                                                                                                                        tGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCAcc                               >  2:3045176/1‑150 (MQ=255)
                                                                                                                          cGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCtt                             >  1:1184128/1‑150 (MQ=255)
                                                                                                                             agagTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCaa                          <  2:4029722/150‑1 (MQ=255)
                                                                                                                               agTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACg                        >  1:751092/1‑150 (MQ=255)
                                                                                                                                     tacGGGATACCGAGCTTCTCGGTAACCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCgg                  >  2:4180160/3‑150 (MQ=255)
                                                                                                                                      gCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGt                 >  1:136500/1‑150 (MQ=255)
                                                                                                                                      gCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGt                 >  1:2136360/1‑150 (MQ=255)
                                                                                                                                      gCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGt                 >  1:355892/1‑150 (MQ=255)
                                                                                                                                         ggatACCGAGCTTCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGAt              >  2:2103295/4‑150 (MQ=255)
                                                                                                                                         ggatACCGAGCTTCTCGGTAACCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACg                                                                                           >  1:1531539/4‑73 (MQ=25)
                                                                                                                                         ggatACCGAGCTTCTCGGTAACCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGAt              >  2:3276073/4‑150 (MQ=255)
                                                                                                                                         ggTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGAt              >  1:3721480/1‑150 (MQ=255)
                                                                                                                                          gatACCGAGCTTCTCGGTAACCCGCACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATc             >  2:2129465/3‑150 (MQ=255)
                                                                                                                                          gTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATc             >  1:4003311/1‑150 (MQ=255)
                                                                                                                                           ttACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCa            <  1:3321545/150‑1 (MQ=255)
                                                                                                                                               cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGa        >  2:1359267/1‑150 (MQ=255)
                                                                                                                                               cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGa        >  2:1245892/1‑150 (MQ=255)
                                                                                                                                               cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGa        >  1:1866200/1‑150 (MQ=255)
                                                                                                                                               cGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGa        >  2:292696/1‑150 (MQ=255)
                                                                                                                                                  gaAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGAtgt     <  1:2096944/150‑1 (MQ=255)
                                                                                                                                                  gaAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGAtgt     <  1:370310/150‑1 (MQ=255)
                                                                                                                                                     gCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCa  >  2:3854513/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCTGATCGATGACCCAGAAGATCATTCGCGTCGGTAACATCGAGATCGCCAACGACAAGCCGTTCGTCCTGTTCGGCGGCATGAACGTTCTGGAATCCCGTGACCTGGCGATGAAGGTCTGCGAAGAGTACGTGCGGGTTACCGAGAAGCTCGGTATCCCGTACGTGTTCAAGGCCAGTTTCGACAAGGCCAACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCCGGTGATCACCGATGTGCA  >  NC_002947/1807983‑1808281

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: