Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 264,177 G→T 5.4% A310S (GCC→TCC)  gabD‑I → succinate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947264,1770GT 5.4% 500.5 / 12.8 149A310S (GCC→TCC) gabD‑Isuccinate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (59/82);  new base T (5/3);  total (64/85)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 155 total aligned reads displayed.

CCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAA  >  NC_002947/264028‑264325
                                                                                                                                                     |                                                                                                                                                    
ccccGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTgg                                                                                                                                                      <  2:878868/150‑1 (MQ=255)
ccccGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTgg                                                                                                                                                      <  2:4195669/150‑1 (MQ=255)
 cccGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGc                                                                                                                                                     >  1:1420307/1‑150 (MQ=255)
 cccGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGc                                                                                                                                                     >  1:4943141/1‑150 (MQ=255)
  ccGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGcc                                                                                                                                                    >  2:2712089/1‑150 (MQ=255)
  ccGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGcc                                                                                                                                                    >  2:1866792/1‑150 (MQ=255)
        aTCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCa                                                                                                                                              >  1:1294339/1‑150 (MQ=255)
        aTCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCa                                                                                                                                              >  1:1294355/1‑150 (MQ=255)
        aTCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCa                                                                                                                                              <  2:2517855/150‑1 (MQ=255)
          cGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGt                                                                                                                                            <  1:2412087/150‑1 (MQ=255)
          cGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGt                                                                                                                                            >  1:4239681/1‑150 (MQ=255)
             gtTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGc                                                                                                                                         <  1:2033860/150‑1 (MQ=255)
               tCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCa                                                                                                                                       >  2:4529660/1‑150 (MQ=255)
                  acgacgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGc                                                                                                                                    <  1:2029679/150‑1 (MQ=255)
                  acgacgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGc                                                                                                                                    <  1:3690564/150‑1 (MQ=255)
                  acgacgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGc                                                                                                                                    <  1:3690630/150‑1 (MQ=255)
                   cgacgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCt                                                                                                                                   >  2:767349/1‑150 (MQ=255)
                       gCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAg                                                                                                                               >  1:2453393/1‑150 (MQ=255)
                         cGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGAt                                                                                                                             <  1:2519111/150‑1 (MQ=255)
                         cGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGAt                                                                                                                             <  1:749825/150‑1 (MQ=255)
                          gACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATc                                                                                                                            >  1:1627621/1‑150 (MQ=255)
                          gACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATc                                                                                                                            >  1:1627380/1‑150 (MQ=255)
                           aCCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCg                                                                                                                           <  1:1033305/150‑1 (MQ=255)
                           aCCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCg                                                                                                                           <  1:4457006/150‑1 (MQ=255)
                           aCCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCg                                                                                                                           <  1:5079819/150‑1 (MQ=255)
                             cTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGt                                                                                                                         >  2:3573236/1‑150 (MQ=255)
                             cTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGt                                                                                                                         >  2:56306/1‑150 (MQ=255)
                                gACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAAc                                                                                                                      >  1:3948916/1‑150 (MQ=255)
                                  cAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGACAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACgg                                                                                                                    >  2:2825377/1‑150 (MQ=255)
                                   aaGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACggt                                                                                                                   >  1:3511453/1‑149 (MQ=255)
                                     ggCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCt                                                                                                                 <  1:4130546/150‑1 (MQ=255)
                                     ggCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCt                                                                                                                 >  2:1111484/1‑150 (MQ=255)
                                     ggCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCt                                                                                                                 <  2:1793626/150‑1 (MQ=255)
                                      gCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTg                                                                                                                <  2:871804/150‑1 (MQ=255)
                                      gCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACAGTATCTACGAGCAGGACGGAGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTg                                                                                                                >  2:1823443/1‑150 (MQ=255)
                                        ggTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGga                                                                                                              >  1:3474875/1‑150 (MQ=255)
                                         gTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGgaa                                                                                                             <  1:436248/150‑1 (MQ=255)
                                          tCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGgaag                                                                                                            >  2:1489979/1‑150 (MQ=255)
                                           cGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGgaaga                                                                                                           <  1:2896935/150‑1 (MQ=255)
                                            gAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGgaagaa                                                                                                          <  2:4569253/150‑1 (MQ=255)
                                            gAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGgaagaa                                                                                                          <  1:3201447/150‑1 (MQ=255)
                                              gggCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAgg                                                                                                        >  2:2422001/1‑150 (MQ=255)
                                               ggCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCACGCTGAAGATCGGTAACGGCCTGGAAGAAGGc                                                                                                       <  1:2467682/150‑1 (MQ=255)
                                               ggCGCGATCATATCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGc                                                                                                       >  2:2023795/1‑150 (MQ=255)
                                                gcgcGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGcc                                                                                                      >  2:3708538/1‑149 (MQ=255)
                                                    gATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGcaccac                                                                                                  <  2:4752472/150‑1 (MQ=255)
                                                    gATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGcaccac                                                                                                  <  2:4752501/150‑1 (MQ=255)
                                                      tcatcTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGcaccacca                                                                                                >  1:240584/1‑150 (MQ=255)
                                                            ccAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGcc                                                                                          <  2:123899/150‑1 (MQ=255)
                                                            ccAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGcc                                                                                          >  2:939710/1‑150 (MQ=255)
                                                            ccAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGcc                                                                                          <  1:4025179/150‑1 (MQ=255)
                                                             cAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGccc                                                                                         >  2:1231353/1‑150 (MQ=255)
                                                                ggaCCGTACCCCCGGCCAGACCTGCGTCTGCGCCACCCGGCTCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGACGGCACCACCACTGGCCCGCt                                                                                      <  2:3854256/148‑1 (MQ=255)
                                                                 tACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTg                                                                                     >  2:854582/1‑150 (MQ=255)
                                                                    cGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATc                                                                                  >  2:2349112/1‑150 (MQ=255)
                                                                     gcaacaacGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCg                                                                                 >  1:2793585/3‑150 (MQ=255)
                                                                        acaacGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGATg                                                                              <  1:1079911/150‑1 (MQ=255)
                                                                         caccGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGATga                                                                             <  2:4844797/147‑2 (MQ=255)
                                                                          aacGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGATgag                                                                            <  1:1545036/150‑3 (MQ=255)
                                                                             ggCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAg                                                                         <  2:1242078/150‑1 (MQ=255)
                                                                             ggCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAg                                                                         <  1:3508096/150‑1 (MQ=255)
                                                                              gCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAgg                                                                        >  1:301816/1‑150 (MQ=255)
                                                                               ccAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGATgaggcggc                                                                       >  1:2053030/1‑143 (MQ=255)
                                                                                cAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGTACGGCGTATACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGTAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCt                                                                      <  1:337368/150‑1 (MQ=255)
                                                                                cAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCt                                                                      >  1:1944774/1‑150 (MQ=255)
                                                                                   aCCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTc                                                                   <  2:3557138/150‑1 (MQ=255)
                                                                                        cGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCaa                                                              <  2:1957408/150‑1 (MQ=255)
                                                                                         gTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAg                                                             <  1:897791/150‑1 (MQ=255)
                                                                                            tGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTc                                                          >  2:5006577/1‑150 (MQ=255)
                                                                                             gcgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTcc                                                         <  2:4008619/150‑1 (MQ=255)
                                                                                              cgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCa                                                        <  2:1234806/150‑1 (MQ=255)
                                                                                              cgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCa                                                        <  2:1487734/150‑1 (MQ=255)
                                                                                              cgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCa                                                        <  2:1631468/150‑1 (MQ=255)
                                                                                                ccAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAgg                                                      >  2:4933152/1‑150 (MQ=255)
                                                                                                    ccGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAaca                                                  >  1:4932311/1‑150 (MQ=255)
                                                                                                         tCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCg                                             >  2:4164722/1‑150 (MQ=255)
                                                                                                          cTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGa                                            <  1:4102834/150‑1 (MQ=255)
                                                                                                          cTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGa                                            <  1:283016/150‑1 (MQ=255)
                                                                                                           tACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAg                                           <  1:4917792/150‑1 (MQ=255)
                                                                                                            aCGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAgg                                          <  1:3740733/150‑1 (MQ=255)
                                                                                                             cGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGa                                         <  1:4240025/150‑1 (MQ=255)
                                                                                                              gTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGAc                                        >  1:5037878/1‑150 (MQ=255)
                                                                                                                           tCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAg                           >  2:5117732/1‑150 (MQ=255)
                                                                                                                            cTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAgg                          <  1:4529660/150‑1 (MQ=255)
                                                                                                                             tACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGgc                         >  1:4361887/1‑150 (MQ=255)
                                                                                                                             tACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGgc                         <  1:3566310/150‑1 (MQ=255)
                                                                                                                               cgacgCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGgcgc                       <  2:3553991/150‑1 (MQ=255)
                                                                                                                                gacgCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGcc                      >  2:3694771/1‑150 (MQ=255)
                                                                                                                                gacgCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGcc                      >  1:112801/1‑150 (MQ=255)
                                                                                                                                  cgcgTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGAt                                                                               >  2:4852241/1‑91 (MQ=255)
                                                                                                                                  cgcgTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGAt                                                                               <  1:4852241/91‑1 (MQ=255)
                                                                                                                                   gcgTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCaaa                   <  1:1853090/150‑1 (MQ=255)
                                                                                                                                   gcgTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCaaa                   <  2:4113681/150‑1 (MQ=255)
                                                                                                                                      ttCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTg                <  1:3140683/150‑1 (MQ=255)
                                                                                                                                      ttCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCACCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTg                <  2:2734884/150‑1 (MQ=255)
                                                                                                                                      ttCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCACCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTg                <  2:2275999/150‑1 (MQ=255)
                                                                                                                                         gCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTg             >  2:2552215/1‑150 (MQ=255)
                                                                                                                                           cGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTc           >  2:4994041/1‑150 (MQ=255)
                                                                                                                                               aaGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGt       <  1:1263880/150‑1 (MQ=255)
                                                                                                                                                    ggCCGCTGCAGTTGCCAAGCTGAAGATCGGTATCGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCATAGGCGCCAAAGTGCTGTCCGGTGGCaa  <  2:3483581/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAA  >  NC_002947/264028‑264325

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: