Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 1,271,763 Δ1 bp intergenic (+20/+69) PA14_RS05965 → / ← PA14_RS05970 YidB family protein/aldo/keto reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084631,271,7610G.99.8% 2261.7 / ‑1.1 462intergenic (+18/+71)PA14_RS05965/PA14_RS05970YidB family protein/aldo/keto reductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (253/206);  total (256/206)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 497 total aligned reads displayed.

GGAGTCAGCCAGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTGGCGGCGATCAGGTC  >  NC_008463/1271612‑1271906
                                                                                                                                                     |                                                                                                                                                 
ggAGTCAGCCAGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGACCTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCg                                                                                                                                                   <  2:217019/150‑1 (MQ=255)
ggAGTCAGCCAGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCg                                                                                                                                                   >  2:3967127/1‑150 (MQ=255)
ggAGTCAGCCAGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCg                                                                                                                                                   >  2:9531739/1‑150 (MQ=255)
 gAGTCAGCCAGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCgg                                                                                                                                                  >  1:8481601/1‑150 (MQ=255)
   gTCAGCCAGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCggag                                                                                                                                                <  1:8894423/150‑3 (MQ=255)
     cAGCCAGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCggagcg                                                                                                                                              <  1:674877/150‑5 (MQ=255)
     cAGCCAGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCggagcg                                                                                                                                              <  1:4860516/150‑5 (MQ=255)
       gCCAGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCGAAGGTTTTCCCGTCGATgcggagcgga                                                                                                                                           >  1:5996049/1‑150 (MQ=255)
          aGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAg                                                                                                                                        <  2:684004/150‑1 (MQ=255)
            ggCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGcc                                                                                                                                      >  1:2159498/1‑150 (MQ=255)
              cGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCcgc                                                                                                                                    <  2:863850/150‑1 (MQ=255)
                aGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATCGCCGCGa                                                                                                                                  >  1:1373608/1‑150 (MQ=255)
                aGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGa                                                                                                                                  >  1:38134/1‑150 (MQ=255)
                 ggcggcCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGAt                                                                                                                                 >  2:9441596/1‑150 (MQ=255)
                 ggcggcCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGAt                                                                                                                                 <  1:7256804/150‑1 (MQ=255)
                       cGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGtt                                                                                                                           <  2:8047590/150‑1 (MQ=255)
                        gAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTa                                                                                                                          <  2:10636163/150‑1 (MQ=255)
                        gAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTa                                                                                                                          >  2:5869405/1‑150 (MQ=255)
                         aaCACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTAt                                                                                                                         >  1:6622086/1‑150 (MQ=255)
                         aaCACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTAt                                                                                                                         >  1:9776294/1‑150 (MQ=255)
                          acacCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATc                                                                                                                        <  1:2960455/150‑1 (MQ=255)
                              cTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCt                                                                                                                    >  2:9785599/1‑150 (MQ=255)
                                 agcagcCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCa                                                                                                                 <  1:10687732/150‑1 (MQ=255)
                                   cagcCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTGCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACt                                                                                                               >  1:2077681/1‑150 (MQ=255)
                                    agcCTGCTGCCGGAACTGGTCCACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGAGATCGGTTCGCTGCTGGCGCGCTTGAGCTATGGTGTTCCCGTCGTTGCGGAGCGGTTAGCCGCGTTACGGTGTTCCGGTCCAcgc                                                                                                              >  2:3127531/1‑148 (MQ=255)
                                    agcCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTc                                                                                                              >  1:271026/1‑150 (MQ=255)
                                    agcCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTc                                                                                                              <  2:5704689/150‑1 (MQ=255)
                                     gcCTGCTGCCGGACCTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGGTATCCGCTCCACTCa                                                                                                             >  2:6049602/1‑150 (MQ=255)
                                      ccTGCTGCCGGCACTGGTCAACACGCTGCCGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCaa                                                                                                            <  1:10738117/150‑1 (MQ=255)
                                      ccTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCaa                                                                                                            >  2:410459/1‑150 (MQ=255)
                                       ctgctgCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAAc                                                                                                           <  1:7445287/150‑1 (MQ=255)
                                        tgctgcCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACg                                                                                                          >  1:5573732/1‑150 (MQ=255)
                                        tgctgcCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACg                                                                                                          <  1:7935577/150‑1 (MQ=255)
                                           tgcCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACgggg                                                                                                       <  1:1551867/150‑1 (MQ=255)
                                            gcCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGc                                                                                                      >  2:6610408/1‑150 (MQ=255)
                                             ccGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCg                                                                                                     <  1:1221384/150‑1 (MQ=255)
                                              cGGAACTGGTCAACACGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGa                                                                                                    <  2:1019384/150‑1 (MQ=255)
                                                 aaCTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGgc                                                                                                 >  2:2047501/1‑150 (MQ=255)
                                                    tGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGgccgc                                                                                              >  2:6942866/1‑150 (MQ=255)
                                                       tCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCt                                                                                           <  1:5070149/150‑1 (MQ=255)
                                                         aacaaGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGTTACGGGTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCg                                                                                         >  2:3706387/1‑150 (MQ=255)
                                                         aacaaGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCg                                                                                         <  1:7790276/150‑1 (MQ=255)
                                                         aacaaGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCg                                                                                         <  1:585801/150‑1 (MQ=255)
                                                          acaaGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGa                                                                                        <  2:5559610/150‑1 (MQ=255)
                                                           caaGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGAc                                                                                       <  2:4825419/150‑1 (MQ=255)
                                                             aGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACgg                                                                                     >  1:909110/1‑150 (MQ=255)
                                                             aGCGGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTGTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACgg                                                                                     >  2:5155716/1‑150 (MQ=255)
                                                               cTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCa                                                                                   <  1:6125339/150‑1 (MQ=255)
                                                                tGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAg                                                                                  <  1:6526051/150‑1 (MQ=255)
                                                                 gACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAgg                                                                                 >  1:3136032/1‑150 (MQ=255)
                                                                 gACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAgg                                                                                 >  2:2506599/1‑150 (MQ=255)
                                                                  aCGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAggt                                                                                <  1:10068212/150‑1 (MQ=255)
                                                                      cGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAgggcggt                                                                            >  1:24722/1‑150 (MQ=255)
                                                                      cGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAAAGGGGCGAGGCCGCCCTCGACGGCAggtcggt                                                                            >  2:3797093/1‑150 (MQ=255)
                                                                        gACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCa                                                                          <  1:11031891/150‑1 (MQ=255)
                                                                         aCGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCaa                                                                         <  1:4418515/150‑1 (MQ=255)
                                                                          cGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAAt                                                                        <  1:3468312/150‑1 (MQ=255)
                                                                           ggCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATa                                                                       <  1:925133/150‑1 (MQ=255)
                                                                                 gcgcCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAggcgct                                                                 <  1:10168261/150‑1 (MQ=255)
                                                                                     cgcAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAggcgctggcg                                                             >  2:1743320/1‑150 (MQ=255)
                                                                                     cgcAGGGCGACCTGGACATCGGTTAGCTGCAGGCGCGCTTCAGCAAAGGTGTTCCCGTCGATGCGGAGCGGATAGCCTCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGTTCGGTCAATAggcgctggcg                                                             >  1:5391810/1‑150 (MQ=255)
                                                                                       cAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAggcgctggcgct                                                           >  1:86601/1‑150 (MQ=255)
                                                                                        aGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGgcgctggcgctg                                                          >  1:8308900/1‑150 (MQ=255)
                                                                                         gggCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGt                                                         >  1:17517/1‑150 (MQ=255)
                                                                                              aCCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATCGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGGCAATAGGCGCTGGCGCTGTCCAgg                                                    >  2:11010293/1‑150 (MQ=255)
                                                                                              aCCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAgg                                                    <  2:9118874/150‑1 (MQ=255)
                                                                                                cTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGGTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATTGGCGCTGGCGCTGTCCAGGCg                                                  >  2:7336767/1‑150 (MQ=255)
                                                                                                 tGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCgg                                                 >  1:3916316/1‑150 (MQ=255)
                                                                                                 tGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGGCCAGGCgg                                                 >  2:1675848/1‑150 (MQ=255)
                                                                                                   gACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTg                                               <  2:8868692/150‑1 (MQ=255)
                                                                                                   gACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTg                                               >  2:5656325/1‑150 (MQ=255)
                                                                                                      aTCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGAt                                            <  2:1181531/150‑1 (MQ=255)
                                                                                                       tCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATc                                           <  2:3496670/150‑1 (MQ=255)
                                                                                                         ggTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGc                                         <  1:5319295/150‑1 (MQ=255)
                                                                                                          gTTCGCTGCTGGCGCGCTTCAGGTAAGCGTATCCCGTCGATGCGGAGACGATAGCCGCGCTACGGGTATCCGCTCCACACCACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCGGGCGGTGGATCggg                                        >  2:8220195/1‑148 (MQ=255)
                                                                                                          gTTCGCTGCTGGCGCGCTTCAGCTCAGGGTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCGGGCGCTGTCCAGGCGGTGGATCGcc                                        >  2:2644789/1‑150 (MQ=255)
                                                                                                          gTTCGCGGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGAGAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGcc                                        <  2:10883387/150‑1 (MQ=255)
                                                                                                            tCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTg                                      <  2:7642357/150‑1 (MQ=255)
                                                                                                            tCGCTGCTGGCGCGCGTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTg                                      <  1:9901994/150‑1 (MQ=255)
                                                                                                               ctgctgGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTtc                                   >  1:10444926/1‑150 (MQ=255)
                                                                                                               ctgctgGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTtc                                   >  2:10303392/1‑150 (MQ=255)
                                                                                                                tgctgGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTtcg                                  <  2:4726142/150‑1 (MQ=255)
                                                                                                                  ctgGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTtcgtc                                <  1:10851801/150‑1 (MQ=255)
                                                                                                                     gcgcgcTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCACTAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAg                             >  2:9292456/1‑150 (MQ=255)
                                                                                                                     gcgcgcTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAg                             <  2:2867759/150‑1 (MQ=255)
                                                                                                                         gcTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTc                         >  2:1394362/1‑150 (MQ=255)
                                                                                                                         gcTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTc                         >  2:8466090/1‑150 (MQ=255)
                                                                                                                          cTTCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCa                        >  1:413589/1‑150 (MQ=255)
                                                                                                                            tCAGCTAAGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAgg                      <  2:7771949/150‑1 (MQ=255)
                                                                                                                                  aaGGTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGATTGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAggtgggtg                >  2:4603471/1‑150 (MQ=255)
                                                                                                                                    ggTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTggc              >  2:6343180/1‑150 (MQ=255)
                                                                                                                                    ggTTTTCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTggc              >  2:4437169/1‑150 (MQ=255)
                                                                                                                                      ttttCCCGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTggcgg            >  1:3614980/1‑150 (MQ=255)
                                                                                                                                           ccGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTGGCGGCGATc       >  2:4190800/1‑150 (MQ=255)
                                                                                                                                           ccGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTGGCGGCGATc       <  1:2540990/150‑1 (MQ=255)
                                                                                                                                            cGTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTGGCGGCGATCa      <  2:10062168/150‑1 (MQ=255)
                                                                                                                                             gTCGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCTGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTGGCGGCGATCAg     >  2:2977308/1‑150 (MQ=255)
                                                                                                                                               cGATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTGGCGGCGATCAGgt   <  1:4272082/150‑1 (MQ=255)
                                                                                                                                                gATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTGGCGGCGATCAGgtc  >  1:4563217/1‑150 (MQ=255)
                                                                                                                                                gATGCGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGGCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTGCAGGTGGGTGGCGGCGATCAGgtc  >  2:4350146/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGAGTCAGCCAGGGCGAGGCGGCCGAACACCTGAGCAGCCTGCTGCCGGAACTGGTCAACAAGCTGACGCCGGACGGCCAGGCGCCGCAGGGCGACCTGGACATCGGTTCGCTGCTGGCGCGCTTCAGCTAAGGTTTTCCCGTCGATGCGGGAGCGGATAGCCGCGATACGGTTATCCGCTCCACTCAACGGGGCGAGGCCGCCCTCGACGGCAGGTCGGTCAATAGGCGCTGGCGCTGTCCAGGCGGTGGATCGCCTGTTCGTCGAGCTTCAGGTGGGTGGCGGCGATCAGGTC  >  NC_008463/1271612‑1271906

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