Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,866,730 T→G A733A (GCA→GCC PA14_RS26890 ← ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,866,7300TG99.5% 2080.4 / NA 550A733A (GCA→GCCPA14_RS26890ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  major base G (338/207);  minor base C (3/0);  total (341/209)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 550 total aligned reads displayed.

CGGATGAGCGGATTATTGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCTGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGACAGGCGATGG  >  NC_008463/5866582‑5866877
                                                                                                                                                    |                                                                                                                                                   
  gATGAGCGGATTATTGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGgcgc                                                                                                                                                 >  1:5363475/1‑151 (MQ=255)
  gATGAGCGGATTATTGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGgcgc                                                                                                                                                 <  2:7140425/151‑1 (MQ=255)
   aTGAGCGGATTATTGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGcc                                                                                                                                                <  1:11043900/151‑1 (MQ=255)
     gAGCGGATTATTGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGt                                                                                                                                              >  1:10065945/1‑151 (MQ=255)
      aGCGGATTATTGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTc                                                                                                                                             >  1:7587092/1‑151 (MQ=255)
      aGCGGATTATTGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGGTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTc                                                                                                                                             <  2:2529833/151‑1 (MQ=255)
          gATTATTGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCgcgg                                                                                                                                         >  2:545431/1‑151 (MQ=255)
          gATTATTGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCgcgg                                                                                                                                         >  2:6023878/1‑151 (MQ=255)
           attattGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCgcggc                                                                                                                                        >  2:6621328/1‑151 (MQ=255)
                 ggcggcGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGc                                                                                                                                  <  1:78853/151‑1 (MQ=255)
                     gcggATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTAcc                                                                                                                              >  2:9767842/1‑151 (MQ=255)
                     gcggATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTAcc                                                                                                                              >  2:4658395/1‑151 (MQ=255)
                      cggATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTAccc                                                                                                                             >  1:9977442/1‑151 (MQ=255)
                       ggATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCa                                                                                                                            <  2:2150936/151‑1 (MQ=255)
                        gATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAg                                                                                                                           <  1:8129544/151‑1 (MQ=255)
                                aCGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTgcgc                                                                                                                   >  1:8905251/1‑151 (MQ=255)
                                aCGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTgcgc                                                                                                                   >  1:2421363/1‑151 (MQ=255)
                                 cGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCGGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGcc                                                                                                                  <  1:8176242/151‑1 (MQ=255)
                                   ccGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGcccc                                                                                                                <  1:2233310/151‑1 (MQ=255)
                                    cGTTGAACGCGCGGATCGAAGATGTAAAATGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCa                                                                                                               >  2:592981/1‑151 (MQ=255)
                                    cGTAGAACGCGCGGATCGAAGATGAAAAAAGCTGTAAGGAAGTTGATCGCCGCGACGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGACGCGGCGCAGGCTACCCAGGCGAGCGCCCCa                                                                                                               >  2:3404294/1‑151 (MQ=255)
                                      ttGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAgg                                                                                                             <  1:10681622/151‑1 (MQ=255)
                                      ttGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAgg                                                                                                             <  2:276633/151‑1 (MQ=255)
                                        gAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCt                                                                                                           <  1:875982/151‑1 (MQ=255)
                                        gAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCt                                                                                                           <  1:7057589/151‑1 (MQ=255)
                                                 gATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCg                                                                                                  >  2:1651668/1‑151 (MQ=255)
                                                  aTCGAAGATGTAAAAAGCTGTAAGGACGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGa                                                                                                 >  2:5808681/1‑151 (MQ=255)
                                                  aTCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGa                                                                                                 >  2:8057534/1‑151 (MQ=255)
                                                       aGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCg                                                                                            >  2:1290394/1‑151 (MQ=255)
                                                       aGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCg                                                                                            >  2:7841905/1‑151 (MQ=255)
                                                        gATGTAAAAAGCTGTAAGGAAGTTGGTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCgg                                                                                           <  2:6340732/151‑1 (MQ=255)
                                                         aTGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGc                                                                                          >  2:318021/1‑151 (MQ=255)
                                                         aTGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGc                                                                                          >  1:7390342/1‑151 (MQ=255)
                                                          tGTAAAAAGCTGTAAGGACGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGcc                                                                                         >  1:853821/1‑151 (MQ=255)
                                                          tGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGcc                                                                                         >  2:8888260/1‑151 (MQ=255)
                                                          tGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGcc                                                                                         <  1:4999725/151‑1 (MQ=255)
                                                             aaaaaGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCcgcg                                                                                      >  1:4520392/1‑151 (MQ=255)
                                                               aaaGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGcc                                                                                    >  1:6458757/1‑151 (MQ=255)
                                                                aaGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGccc                                                                                   >  1:6880861/1‑151 (MQ=255)
                                                                    tGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAggg                                                                               >  1:7863497/1‑151 (MQ=255)
                                                                       aaGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCt                                                                            >  1:6122429/1‑151 (MQ=255)
                                                                        aGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTc                                                                           <  2:7997343/151‑1 (MQ=255)
                                                                        aGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTc                                                                           <  2:10084745/151‑1 (MQ=255)
                                                                          gAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTcgc                                                                         <  1:9541867/151‑1 (MQ=255)
                                                                          gAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTcgc                                                                         >  1:3938545/1‑151 (MQ=255)
                                                                          gAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTcgc                                                                         <  1:7722793/151‑1 (MQ=255)
                                                                           aaGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGCGCGTCCGAGCCGTCGCGGCTCCCGCTCCCCCCCAGTGACCGACGGGCGGCCGGGTCGGCCCCGGCCGCGCACTGGGACTcgcg                                                                        >  2:3516282/1‑151 (MQ=255)
                                                                             gttgttCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCt                                                                      >  2:3959866/1‑151 (MQ=255)
                                                                              ttgttCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTc                                                                     >  1:8985271/1‑151 (MQ=255)
                                                                                gttCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCa                                                                   <  1:3591024/151‑1 (MQ=255)
                                                                                gttCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCa                                                                   >  2:4264511/1‑151 (MQ=255)
                                                                                gttCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTGCGCCGAGGACGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCGACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCa                                                                   <  2:8387002/151‑1 (MQ=255)
                                                                                gttCGCCGCGCCGCTCAACCGGCCGCGGAGTCGGCGTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTccc                                                                   >  2:5329254/1‑150 (MQ=255)
                                                                                  tCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGc                                                                 >  1:1999316/1‑151 (MQ=255)
                                                                                    gccgcgccgcTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTc                                                               >  2:18646/1‑151 (MQ=255)
                                                                                    gccgcgccgcTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTc                                                               >  2:15780/1‑151 (MQ=255)
                                                                                    gccgcgccgcTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTc                                                               <  1:9771944/151‑1 (MQ=255)
                                                                                     ccgcgccgcTCAACCGGCCGATGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGCGCTTCGGCGCCGTCGCGGCGAAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCg                                                              >  2:9024053/1‑151 (MQ=255)
                                                                                     ccgcgccgcTCAACCGGCCACTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCg                                                              >  1:1417729/1‑151 (MQ=255)
                                                                                          ccgcTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCCGGGCCTCGCGCTCCAGCTCGCCCAg                                                         >  2:9464469/1‑151 (MQ=255)
                                                                                           cgcTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAgg                                                        >  1:5066028/1‑151 (MQ=255)
                                                                                            gcTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGc                                                       >  2:421979/1‑151 (MQ=255)
                                                                                               cAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCgg                                                    >  1:5255301/1‑151 (MQ=255)
                                                                                                aaCCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCTGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGgc                                                   >  2:6118087/1‑151 (MQ=255)
                                                                                                aaCCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCTGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCTGGCCGgc                                                   <  2:7316539/151‑1 (MQ=255)
                                                                                                   cGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGcc                                                >  2:4826640/1‑151 (MQ=255)
                                                                                                   cGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGcc                                                <  1:4205625/151‑1 (MQ=255)
                                                                                                    ggCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCg                                               >  2:4466273/1‑151 (MQ=255)
                                                                                                      ccGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGTCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAg                                             <  2:7575321/151‑1 (MQ=255)
                                                                                                       cGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAgc                                            >  1:6991565/1‑151 (MQ=255)
                                                                                                        gCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCACGCGCTCCAGCTCGCCCAGGCCGGCGCCGAgcg                                           <  1:4078486/151‑1 (MQ=255)
                                                                                                           gAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGcca                                        >  1:3079851/1‑151 (MQ=255)
                                                                                                            aGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGccac                                       <  1:2653215/151‑1 (MQ=255)
                                                                                                            aGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGccac                                       <  2:2300523/151‑1 (MQ=255)
                                                                                                             gTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGccacc                                      <  1:2340269/151‑1 (MQ=255)
                                                                                                                 gCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGc                                  >  1:9410144/1‑151 (MQ=255)
                                                                                                                  cTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGcc                                 <  1:6645802/151‑1 (MQ=255)
                                                                                                                    tCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCa                               <  1:3461981/151‑1 (MQ=255)
                                                                                                                     cGCCGAGGAAGCCGGCGCGGCGCAGGGCGTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGAGCCCCAGGCTCGCGGGTCGAGCACGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGGCGAGCGCCACCAGGCCCag                              >  2:10394274/1‑150 (MQ=255)
                                                                                                                      gCCGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATg                             <  1:5865715/151‑1 (MQ=255)
                                                                                                                       ccGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATgg                            <  1:2075546/151‑1 (MQ=255)
                                                                                                                        cGAGGAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGa                           <  2:1421183/151‑1 (MQ=255)
                                                                                                                            gAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTgg                       >  1:11106773/1‑151 (MQ=255)
                                                                                                                            gAAGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCACGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTgt                       <  2:9230078/151‑2 (MQ=255)
                                                                                                                              aGCCGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCa                     <  1:7447373/151‑1 (MQ=255)
                                                                                                                                ccGGCGCGGCGCAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGc                   <  2:6808331/151‑1 (MQ=255)
                                                                                                                                  ggcgcggcgcAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCtt                 <  1:1698989/151‑1 (MQ=255)
                                                                                                                                    cgcggcgcAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTcc               <  1:554156/151‑1 (MQ=255)
                                                                                                                                     gcggcgcAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCg              <  1:10819734/151‑1 (MQ=255)
                                                                                                                                      cggcgcAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCgg             <  1:10339000/151‑1 (MQ=255)
                                                                                                                                       ggcgcAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGa            <  1:1049410/151‑1 (MQ=255)
                                                                                                                                       ggcgcAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGa            >  2:10919245/1‑151 (MQ=255)
                                                                                                                                         cgcAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGACa          >  2:2456344/1‑151 (MQ=255)
                                                                                                                                           cAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGACAgg        <  1:1280193/151‑1 (MQ=255)
                                                                                                                                           cAGGGCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGACAgg        <  2:3031048/151‑1 (MQ=255)
                                                                                                                                             gggCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGACAggcg      <  2:5592369/151‑1 (MQ=255)
                                                                                                                                              ggCTTCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGACAggcga     >  2:11072266/1‑151 (MQ=255)
                                                                                                                                                 ttCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGACAggcgatgg  >  2:10669539/1‑151 (MQ=255)
                                                                                                                                                 ttCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGACAggcgatgg  <  2:8241399/151‑1 (MQ=255)
                                                                                                                                                 ttCGGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGACAggcgatgg  <  1:10572786/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGATGAGCGGATTATTGGCGGCGGATGTAAGACGCCGTTGAACGCGCGGATCGAAGATGTAAAAAGCTGTAAGGAAGTTGTTCGCCGCGCCGCTCAACCGGCCGCTGAGTCGGCTTCGCCGAGGAAGCCGGCGCGGCGCAGGGCTTCTGCGCCGTCGCGGCGCAGGCTACCCAGGCGTGCGCCCCAGGCTCGCGGGTCGACCCCGGCCGCGCCCAGGGCCTCGCGCTCCAGCTCGCCCAGGCCGGCGCCGAGCGCCACCAGGCCCATGGAGTGGCAGCTTCCGGACAGGCGATGG  >  NC_008463/5866582‑5866877

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: