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.1% 2254.3 / 13.2 471pseudogene (228/2518 nt)hrpBATP‑dependent helicase HrpB
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base . (228/234);  minor base T (4/0);  total (233/238)
Fisher's exact test for biased strand distribution p-value = 6.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 488 total aligned reads displayed.

CCCCGCCGGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGGTCGGCGAGACCGTCGGCTACGGGATCCGCCTGGAAAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTCGACG  >  NC_008463/1082773‑1083067
                                                                                                                                                    |                                                                                                                                                  
  ccgccgGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAggtc                                                                                                                                                 <  1:5253124/150‑3 (MQ=255)
  ccgccgGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAggtc                                                                                                                                                 >  1:3504477/1‑148 (MQ=255)
   cgccgGGTGCCGGCAAGACCACCCGGGGGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAggtcg                                                                                                                                                <  2:5707576/150‑4 (MQ=255)
   cgccgGGCGCCGGCACGACCCCCCGGCTGCCGCTGGCCCTGCTCGACGAGCCCCGGATCGCCGGCCAACGCATCCTCCACCTCGAACCGCGCCGCCAGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAggtcg                                                                                                                                                <  2:8813838/150‑4 (MQ=255)
    gccgGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAggtcgg                                                                                                                                               <  1:2478556/150‑5 (MQ=255)
    gccgGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAggtcgg                                                                                                                                               >  1:4315931/1‑146 (MQ=255)
        ggTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCgag                                                                                                                                          <  1:8345580/150‑1 (MQ=255)
         gTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCgaga                                                                                                                                         >  2:9070460/1‑150 (MQ=255)
            ccGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCg                                                                                                                                      <  2:10643468/150‑1 (MQ=255)
            ccGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCg                                                                                                                                      >  2:3849466/1‑150 (MQ=255)
             cGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGt                                                                                                                                     <  2:4163819/150‑1 (MQ=255)
               gCAAGACCCCCCGGGTGCCGCTGGCCCTGCTCGCCGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCg                                                                                                                                   <  1:7722902/150‑1 (MQ=255)
               gCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCg                                                                                                                                   <  2:9270367/150‑1 (MQ=255)
                  aGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCt                                                                                                                                <  1:1828613/150‑1 (MQ=255)
                  aGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCt                                                                                                                                >  2:8165461/1‑150 (MQ=255)
                  aGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGACAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCt                                                                                                                                >  1:6511710/1‑150 (MQ=255)
                        cccGGGTGCCGCTGGCCCTGCTCGACGAGCCCTCGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGa                                                                                                                          <  2:5644607/150‑1 (MQ=255)
                          cGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATc                                                                                                                        <  2:47664/150‑1 (MQ=255)
                           gggTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATcc                                                                                                                       >  1:3620483/1‑150 (MQ=255)
                                  gCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGa                                                                                                                >  2:2880993/1‑150 (MQ=255)
                                  gCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGa                                                                                                                >  1:665949/1‑150 (MQ=255)
                                  gCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATACGCCTGGa                                                                                                                >  2:10495643/1‑150 (MQ=255)
                                   cTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGag                                                                                                               >  1:10106603/1‑149 (MQ=255)
                                   cTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGag                                                                                                               >  1:10382996/1‑149 (MQ=255)
                                     ggCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGagag                                                                                                             >  1:9604706/1‑150 (MQ=255)
                                     ggCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGagag                                                                                                             <  1:2910044/150‑1 (MQ=255)
                                          tGCTCGACGAGCCCTGGCTCGCCGGCCACCGCATCCTCATGCTCGACCCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAggg                                                                                                        <  2:1675677/150‑1 (MQ=255)
                                          tGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGACCCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAggg                                                                                                        <  1:4004212/150‑1 (MQ=255)
                                          tGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAggg                                                                                                        >  1:3825033/1‑150 (MQ=255)
                                           gCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGt                                                                                                       <  2:7539870/150‑1 (MQ=255)
                                           gCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGt                                                                                                       <  1:8813693/150‑1 (MQ=255)
                                           gCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGt                                                                                                       >  2:8579608/1‑150 (MQ=255)
                                             tCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGACCCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCg                                                                                                     <  1:6820945/150‑1 (MQ=255)
                                               gacgaGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCGGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGc                                                                                                   <  1:2062081/150‑1 (MQ=255)
                                                   aGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGa                                                                                               >  1:5753406/1‑150 (MQ=255)
                                                     cccTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAg                                                                                             <  2:4379977/150‑1 (MQ=255)
                                                      ccTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGa                                                                                            <  1:8303561/150‑1 (MQ=255)
                                                          gCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCg                                                                                        >  1:9066362/1‑150 (MQ=255)
                                                          gCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCg                                                                                        <  2:9775652/150‑1 (MQ=255)
                                                           cTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGc                                                                                       <  1:7893329/150‑1 (MQ=255)
                                                            tCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCa                                                                                      <  1:4416580/150‑1 (MQ=255)
                                                                cGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGACGACCCGCATCGa                                                                                  >  1:5162276/1‑150 (MQ=255)
                                                                 ggCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGaa                                                                                 <  1:5879537/150‑1 (MQ=255)
                                                                   ccAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGt                                                                               >  1:10317336/1‑150 (MQ=255)
                                                                    cAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTg                                                                              >  1:8078835/1‑150 (MQ=255)
                                                                    cAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCGGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTg                                                                              <  1:2600058/150‑1 (MQ=255)
                                                                      aCGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGt                                                                            <  1:1430573/150‑1 (MQ=255)
                                                                      aCGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGt                                                                            <  2:1428603/150‑1 (MQ=255)
                                                                       cGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTc                                                                           >  1:10684388/1‑150 (MQ=255)
                                                                           tCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCg                                                                       >  2:3533957/1‑150 (MQ=255)
                                                                               cATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAggg                                                                   <  2:5499873/150‑1 (MQ=255)
                                                                                 tGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTa                                                                 <  2:2247535/150‑1 (MQ=255)
                                                                                   cTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATc                                                               >  2:776960/1‑150 (MQ=255)
                                                                                     cGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGACGTGGTCACCGAGGGTATCCt                                                             >  1:10223804/1‑150 (MQ=255)
                                                                                            cgACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGt                                                      >  1:6002136/1‑150 (MQ=255)
                                                                                             gACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGACCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCCCCGAGGGTATCCTTGCCCGTc                                                     >  1:7523857/1‑150 (MQ=255)
                                                                                               cGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGc                                                   >  1:847117/1‑150 (MQ=255)
                                                                                               cGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGc                                                   >  2:4908009/1‑150 (MQ=255)
                                                                                                 ccTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCt                                                 >  1:2718876/1‑150 (MQ=255)
                                                                                                 ccTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCt                                                 >  1:6333282/1‑150 (MQ=255)
                                                                                                     gcggcgCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCaa                                             >  1:9501567/1‑150 (MQ=255)
                                                                                                       ggcgCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAga                                           <  2:6297415/150‑1 (MQ=255)
                                                                                                       ggagcgGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGTTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGTGTATCCTTGCCCGTCGCCTGCAAga                                           <  1:2210978/147‑1 (MQ=255)
                                                                                                        gcgcGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCCCCGAGGGTATCCTTGCCCGTCGCCTGCAAgac                                          >  2:3172871/1‑150 (MQ=255)
                                                                                                        gcgcGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAgac                                          >  2:8043682/1‑150 (MQ=255)
                                                                                                         cgcgGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAgacg                                         >  2:5053277/1‑150 (MQ=255)
                                                                                                          gcgGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAACAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAgacga                                        <  2:6724025/150‑1 (MQ=255)
                                                                                                           cgGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAgacgac                                       <  2:8248238/150‑1 (MQ=255)
                                                                                                            gggCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGAcc                                      >  1:11138376/1‑150 (MQ=255)
                                                                                                            gggCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGAcc                                      <  2:2530319/150‑1 (MQ=255)
                                                                                                             ggCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGAccc                                     >  1:9238840/1‑150 (MQ=255)
                                                                                                             ggCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGAccc                                     <  2:6206236/150‑1 (MQ=255)
                                                                                                                cgcgGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGTCCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGgc                                  >  2:6002636/1‑150 (MQ=255)
                                                                                                                 gcgGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGgcg                                 >  1:9140056/1‑150 (MQ=255)
                                                                                                                      cGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGa                            <  2:1218558/150‑1 (MQ=255)
                                                                                                                      cGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGa                            >  2:184838/1‑150 (MQ=255)
                                                                                                                         gCGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACgg                         <  1:7117026/150‑1 (MQ=255)
                                                                                                                          cGGCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGgt                        <  2:5911696/150‑1 (MQ=255)
                                                                                                                            gCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGgtgt                      <  2:9625009/150‑1 (MQ=255)
                                                                                                                            gCTGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGgtgt                      <  2:7243664/150‑1 (MQ=255)
                                                                                                                              tGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCg                    <  1:16091/150‑1 (MQ=255)
                                                                                                                              tGGCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCg                    >  1:8639929/1‑150 (MQ=255)
                                                                                                                               ggCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCgg                   <  1:3174440/150‑1 (MQ=255)
                                                                                                                               ggCCGCCGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCgg                   >  2:6836284/1‑150 (MQ=255)
                                                                                                                                gccgccGAGCTGGGCGAACAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGc                  <  1:8019651/150‑1 (MQ=255)
                                                                                                                                gccgccGAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGc                  >  1:2818797/1‑150 (MQ=255)
                                                                                                                                 ccgccgAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGcc                 <  2:6912699/150‑1 (MQ=255)
                                                                                                                                 ccgccgAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGcc                 <  2:10922931/150‑1 (MQ=255)
                                                                                                                                  cgccgAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCt                >  1:9833158/1‑150 (MQ=255)
                                                                                                                                    ccgAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTgg              >  1:5406888/1‑150 (MQ=255)
                                                                                                                                      gAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTg            >  1:552455/1‑150 (MQ=255)
                                                                                                                                      gAGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTg            >  2:3237020/1‑150 (MQ=255)
                                                                                                                                       aGCTGGGCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGa           >  2:3642614/1‑150 (MQ=255)
                                                                                                                                           gggCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCtt       <  1:4753021/150‑1 (MQ=255)
                                                                                                                                           gggCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCtt       <  1:6210999/150‑1 (MQ=255)
                                                                                                                                           gggCGAAAAGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCtt       <  1:7392952/150‑1 (MQ=255)
                                                                                                                                               cgaaaaGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTcga    <  2:2169079/148‑1 (MQ=255)
                                                                                                                                               cgaaaaGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTcga    >  2:2796661/3‑150 (MQ=255)
                                                                                                                                               cgaaaaGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTcga    >  2:10352561/3‑150 (MQ=255)
                                                                                                                                                 aaaaGGTCGGCGAGACCGTCGGCTACCGGATCCGCCTGGAGAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGGCCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTcgacg  >  1:4589384/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCCCGCCGGGTGCCGGCAAGACCACCCGGGTGCCGCTGGCCCTGCTCGACGAGCCCTGGCTCGCCGGCCAACGCATCCTCATGCTCGAACCGCGACGCCTGGCGGCGCGGGCCGCGGCCGAGCGGCTGGCCGCCGAGCTGGGCGAAAAGGGTCGGCGAGACCGTCGGCTACGGGATCCGCCTGGAAAGCAGGGTCGGCCCGAAGACCCGCATCGAAGTGGTCACCGAGGGTATCCTTGCCCGTCGCCTGCAAGACGACCCGGCGCTGGACGGTGTCGGCCTGGTGATCTTCGACG  >  NC_008463/1082773‑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: