Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 1,082,923 Δ1 bp pseudogene (230/2518 nt) hrpB → ATP‑dependent helicase HrpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084631,082,9210G.99.3% 2044.8 / 5.5 419pseudogene (228/2518 nt)hrpBATP‑dependent helicase HrpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (207/208);  minor base T (3/0);  total (210/209)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 447 total aligned reads displayed.

GCCCCGCCGGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGGTCGGCGAGACCGTCGGCTACGGGATCCGCCTGGAAAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTCGACG  >  NC_008463/1082772‑1083067
                                                                                                                                                     |                                                                                                                                                  
  cccGCCGGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAggt                                                                                                                                                  <  2:3508706/150‑2 (MQ=255)
    cgccgGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAggtcg                                                                                                                                                <  1:5297113/150‑4 (MQ=255)
     gccgGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAggtcgg                                                                                                                                               <  2:7817821/150‑5 (MQ=255)
      ccgGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGc                                                                                                                                             <  1:10014346/150‑1 (MQ=255)
        gggTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCga                                                                                                                                           >  2:4438942/1‑150 (MQ=255)
        gggTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCga                                                                                                                                           >  2:2679833/1‑150 (MQ=255)
            gCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCGGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGAcc                                                                                                                                       <  1:2225802/150‑1 (MQ=255)
              cGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGt                                                                                                                                     <  1:9181088/150‑1 (MQ=255)
                    gACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTa                                                                                                                               >  2:5908206/1‑150 (MQ=255)
                       caccCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTAccg                                                                                                                            <  2:9125922/150‑3 (MQ=255)
                        accCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACcgg                                                                                                                           <  2:10583504/150‑1 (MQ=255)
                         cccGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGa                                                                                                                          >  1:4293769/1‑150 (MQ=255)
                           cGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATc                                                                                                                        >  2:10498262/1‑150 (MQ=255)
                            gggTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATcc                                                                                                                       <  1:8201976/150‑1 (MQ=255)
                            gggTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATcc                                                                                                                       <  2:9258472/150‑1 (MQ=255)
                              gTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGc                                                                                                                     >  2:9055142/1‑150 (MQ=255)
                                    cTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGag                                                                                                               >  1:10607400/1‑149 (MQ=255)
                                      ggCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGagag                                                                                                             <  1:2742048/150‑1 (MQ=255)
                                       gCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGc                                                                                                            <  2:706501/150‑1 (MQ=255)
                                        cccTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCa                                                                                                           >  2:4218309/1‑150 (MQ=255)
                                        cccTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCa                                                                                                           <  1:6602648/150‑1 (MQ=255)
                                        cccTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCa                                                                                                           <  2:3543718/150‑1 (MQ=255)
                                         ccTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAg                                                                                                          >  1:2895830/1‑150 (MQ=255)
                                         ccTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAg                                                                                                          <  2:1558672/150‑1 (MQ=255)
                                         ccTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAg                                                                                                          >  2:5600842/1‑150 (MQ=255)
                                          cTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAgg                                                                                                         <  1:419700/150‑1 (MQ=255)
                                            gCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGt                                                                                                       <  2:2433521/150‑1 (MQ=255)
                                            gCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGt                                                                                                       <  1:2371654/150‑1 (MQ=255)
                                            gCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGt                                                                                                       >  1:3841336/1‑150 (MQ=255)
                                            gCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGt                                                                                                       <  1:6172578/150‑1 (MQ=255)
                                            gCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGt                                                                                                       >  2:9471302/1‑150 (MQ=255)
                                              tCGACGAGCCCTGGCTCGCCGGCCACCGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCg                                                                                                     <  1:9695415/150‑1 (MQ=255)
                                               cgacgaGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCgg                                                                                                    <  1:8379857/150‑1 (MQ=255)
                                                   gaGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCg                                                                                                >  2:5731989/1‑150 (MQ=255)
                                                    aGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGa                                                                                               >  1:4096772/1‑150 (MQ=255)
                                                     gCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGaa                                                                                              >  2:7555478/1‑150 (MQ=255)
                                                     gCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGaa                                                                                              >  1:6115790/1‑150 (MQ=255)
                                                     gCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGaa                                                                                              <  1:7701754/150‑1 (MQ=255)
                                                      cccTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAg                                                                                             >  1:5741411/1‑150 (MQ=255)
                                                          ggCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGAccc                                                                                         >  2:2901670/1‑150 (MQ=255)
                                                             tCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCa                                                                                      <  2:6946569/150‑1 (MQ=255)
                                                             tCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCa                                                                                      >  1:9001819/1‑150 (MQ=255)
                                                             tCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCa                                                                                      >  1:4936660/1‑150 (MQ=255)
                                                              cGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCAt                                                                                     >  2:3167896/1‑150 (MQ=255)
                                                               gCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATc                                                                                    >  2:6894859/1‑150 (MQ=255)
                                                               gCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATc                                                                                    >  2:2045372/1‑150 (MQ=255)
                                                               gCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATc                                                                                    >  1:4705094/1‑150 (MQ=255)
                                                                 cGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGa                                                                                  <  1:124527/150‑1 (MQ=255)
                                                                   gCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAg                                                                                >  2:10187778/1‑150 (MQ=255)
                                                                     cAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTg                                                                              >  1:9925853/1‑150 (MQ=255)
                                                                      aaCGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTgg                                                                             >  1:7309479/1‑150 (MQ=255)
                                                                      aaCGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTgg                                                                             >  2:8579764/1‑150 (MQ=255)
                                                                        cGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTc                                                                           >  2:397050/1‑150 (MQ=255)
                                                                         gCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCa                                                                          <  1:8617624/150‑1 (MQ=255)
                                                                            tCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCg                                                                       >  2:192902/1‑150 (MQ=255)
                                                                            tCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCg                                                                       >  2:552304/1‑150 (MQ=255)
                                                                              cTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAg                                                                     >  2:8351078/1‑150 (MQ=255)
                                                                               tCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGTCCCGCATCGAAGGGGTCACCGAgg                                                                    >  1:5061268/1‑150 (MQ=255)
                                                                                  tGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTa                                                                 <  1:3358516/150‑1 (MQ=255)
                                                                                      cGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCt                                                             <  2:6564254/150‑1 (MQ=255)
                                                                                       gAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCtt                                                            <  1:4113152/150‑1 (MQ=255)
                                                                                        aaCCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTg                                                           >  1:1585597/1‑150 (MQ=255)
                                                                                            gcgACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCg                                                       <  1:6554390/150‑1 (MQ=255)
                                                                                             cgACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGt                                                      >  1:5385635/1‑150 (MQ=255)
                                                                                             cgACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGt                                                      <  1:206972/150‑1 (MQ=255)
                                                                                              gACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTc                                                     >  2:4579613/1‑150 (MQ=255)
                                                                                              gACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTc                                                     <  1:10143918/150‑1 (MQ=255)
                                                                                               aCGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCg                                                    >  1:3142272/1‑150 (MQ=255)
                                                                                                cGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGc                                                   >  1:641504/1‑150 (MQ=255)
                                                                                                 gCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGcc                                                  <  1:3631888/150‑1 (MQ=255)
                                                                                                   cTGGCGGCGCGGGCCGCGGCCGAGCGGCGGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTg                                                <  1:9659177/150‑1 (MQ=255)
                                                                                                      gcggcgCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCaa                                             >  1:8852049/1‑150 (MQ=255)
                                                                                                        ggcgCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAga                                           <  1:9426864/150‑1 (MQ=255)
                                                                                                         gcgcGGGCCGCGGCCGCGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAgac                                          <  2:5063738/150‑1 (MQ=255)
                                                                                                         gcgcGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAATACCCGCAACGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAgac                                          >  2:8886424/1‑150 (MQ=255)
                                                                                                         gcgcGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAgac                                          >  1:4527589/1‑150 (MQ=255)
                                                                                                           gcgGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAgacga                                        >  2:1224555/1‑150 (MQ=255)
                                                                                                             gggCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGAcc                                      <  1:6795075/150‑1 (MQ=255)
                                                                                                              ggCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGAccc                                     >  1:7605661/1‑150 (MQ=255)
                                                                                                                ccGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCgg                                   >  2:8229362/1‑150 (MQ=255)
                                                                                                                  gcgGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGgcg                                 >  1:7148114/1‑150 (MQ=255)
                                                                                                                    ggCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCt                               <  2:8509561/150‑1 (MQ=255)
                                                                                                                       cGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGa                            <  2:4955777/150‑1 (MQ=255)
                                                                                                                       cGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGa                            >  1:792790/1‑150 (MQ=255)
                                                                                                                        gAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGAc                           >  2:8662238/1‑150 (MQ=255)
                                                                                                                           cGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGgt                        <  1:594499/150‑1 (MQ=255)
                                                                                                                            ggCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGgtg                       <  2:9624830/150‑1 (MQ=255)
                                                                                                                            ggCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGgtg                       <  2:1654497/150‑1 (MQ=255)
                                                                                                                            ggCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTAGCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGgtg                       <  1:1148854/150‑1 (MQ=255)
                                                                                                                             gCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGgtgt                      <  1:1855615/150‑1 (MQ=255)
                                                                                                                               tGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCg                    >  2:6743063/1‑150 (MQ=255)
                                                                                                                                 gccgccGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGc                  <  2:8100094/150‑1 (MQ=255)
                                                                                                                                  ccgccgAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGcc                 >  1:3702653/1‑150 (MQ=255)
                                                                                                                                   cgccgAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCt                <  2:2354355/150‑1 (MQ=255)
                                                                                                                                         gCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGAt          >  1:5875402/1‑150 (MQ=255)
                                                                                                                                         gCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGAt          <  2:644032/150‑1 (MQ=255)
                                                                                                                                           tGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCt        <  1:5611248/150‑1 (MQ=255)
                                                                                                                                             ggCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTc      >  1:9990303/1‑150 (MQ=255)
                                                                                                                                                 gaaaaGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTcgac   <  2:6297198/149‑1 (MQ=255)
                                                                                                                                                 gaaaaGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTcgac   <  1:1524820/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCCCCGCCGGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGGTCGGCGAGACCGTCGGCTACGGGATCCGCCTGGAAAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTCGACG  >  NC_008463/1082772‑1083067

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: