Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,957,549 2 bp→AT intergenic (‑371/‑300) PA14_RS22700 ← / → PA14_RS22705 16S ribosomal RNA/exonuclease SbcCD subunit D C‑terminal domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,957,5490GA99.6% 409.1 / NA 255intergenic (‑371/‑301)PA14_RS22700/PA14_RS2270516S ribosomal RNA/exonuclease SbcCD subunit D C‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (98/156);  minor base C (0/1);  total (98/157)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084634,957,5500CT99.2% 404.2 / NA 255intergenic (‑372/‑300)PA14_RS22700/PA14_RS2270516S ribosomal RNA/exonuclease SbcCD subunit D C‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (97/156);  minor base A (1/1);  total (98/157)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

TTGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGATGAAGCTTCGTTTCGAGCCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACTG  >  NC_008463/4957454‑4957694
                                                                                               ||                                                                                                                                                
  gcgcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTcc                                                                                          >  1:7831784/1‑151 (MQ=9)
  gcgcCGACCTTCAAACCATCGTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTcc                                                                                          <  2:560619/151‑1 (MQ=9)
   cgcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCt                                                                                         <  1:8777724/151‑1 (MQ=9)
   cgcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCt                                                                                         <  1:4692374/151‑1 (MQ=9)
    gcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTg                                                                                        >  2:1626524/1‑151 (MQ=9)
    gcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTg                                                                                        >  1:11033184/1‑151 (MQ=9)
       gACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCACGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGga                                                                                     >  2:4017155/1‑151 (MQ=9)
       gACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGga                                                                                     >  2:3887943/1‑151 (MQ=9)
       gACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGga                                                                                     <  2:7406280/151‑1 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    >  2:6577572/1‑151 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    >  2:1921413/1‑151 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    <  1:9443047/151‑1 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    >  2:9009275/1‑151 (MQ=9)
         ccTTCAAACCCTCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgaga                                                                                   <  1:7753877/151‑1 (MQ=9)
         ccTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgaga                                                                                   <  1:5465212/151‑1 (MQ=9)
           ttCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagaga                                                                                 <  1:6472620/151‑1 (MQ=9)
           ttCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagaga                                                                                 >  2:263111/1‑151 (MQ=9)
           ttCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagaga                                                                                 >  1:3247923/1‑151 (MQ=9)
           ttCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagaga                                                                                 <  2:11298660/151‑1 (MQ=9)
            tCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGaa                                                                                >  2:11540247/1‑151 (MQ=9)
            tCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGaa                                                                                <  1:11196817/151‑1 (MQ=9)
            tCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGaa                                                                                >  2:8185206/1‑151 (MQ=9)
            tCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAAACCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGaa                                                                                >  2:8586447/1‑151 (MQ=9)
              aatccttCTTCCGATCAGCGGGAGGCGGATTCTACAGCGTTCAAATTCGCTGTCTAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGTCCGAACCGCTGCGGCGTAAAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGTATc                                                                              <  2:6997041/145‑1 (MQ=11)
               aaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATcc                                                                             <  1:9045910/151‑1 (MQ=9)
                aCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCg                                                                            <  2:10242027/151‑1 (MQ=9)
                 ccATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGt                                                                           <  2:2444086/151‑1 (MQ=9)
                  cATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTa                                                                          <  2:8007286/151‑1 (MQ=9)
                   aTCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTaa                                                                         <  2:3264761/151‑1 (MQ=9)
                   aTCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTaa                                                                         <  1:2701063/151‑1 (MQ=9)
                   aTCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTaa                                                                         <  1:4914295/151‑1 (MQ=9)
                    tCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTaaa                                                                        <  1:3343774/151‑1 (MQ=9)
                    tCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTaaa                                                                        <  2:7741432/151‑1 (MQ=9)
                     cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAt                                                                       <  1:958362/151‑1 (MQ=9)
                     cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAt                                                                       >  1:5145632/1‑151 (MQ=9)
                     cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAt                                                                       <  1:4950137/151‑1 (MQ=9)
                     cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAt                                                                       <  1:1057977/151‑1 (MQ=9)
                     cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAt                                                                       <  1:3583629/151‑1 (MQ=9)
                     cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAt                                                                       <  1:3577877/151‑1 (MQ=9)
                     cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAt                                                                       <  1:1739429/151‑1 (MQ=9)
                        ccGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCt                                                                    <  1:1025379/151‑1 (MQ=9)
                         cGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCtt                                                                   <  2:9527547/151‑1 (MQ=9)
                          gATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTg                                                                  >  1:2277815/1‑151 (MQ=9)
                           aTCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGa                                                                 >  2:1079342/1‑151 (MQ=9)
                            tCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAt                                                                <  1:4064139/151‑1 (MQ=9)
                             cAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtt                                                               <  2:6118283/151‑1 (MQ=9)
                                cGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTc                                                            <  1:1049194/151‑1 (MQ=9)
                                 gggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCa                                                           >  1:11418401/1‑151 (MQ=9)
                                  ggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCaa                                                          <  1:11552261/151‑1 (MQ=9)
                                   gAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAg                                                         >  2:514328/1‑151 (MQ=9)
                                   gAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAg                                                         >  2:3476828/1‑151 (MQ=9)
                                    aGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGACCAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTGGATTTTCAAGc                                                        >  1:3212825/1‑151 (MQ=255)
                                    aGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGc                                                        <  1:7461119/151‑1 (MQ=9)
                                     ggCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCa                                                       >  1:5964490/1‑151 (MQ=9)
                                      gcgcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAt                                                      <  1:1519233/151‑1 (MQ=9)
                                       cgcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGTCCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAtt                                                     <  1:11517339/151‑1 (MQ=9)
                                       cgcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAtt                                                     <  1:9057481/151‑1 (MQ=9)
                                        gcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAttt                                                    <  2:2371311/151‑1 (MQ=9)
                                         cATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTa                                                   >  2:10666566/1‑151 (MQ=9)
                                         cATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTa                                                   >  1:4275382/1‑151 (MQ=9)
                                            tCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActc                                                <  1:6166045/151‑1 (MQ=9)
                                             cTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActct                                               >  2:2262536/1‑151 (MQ=9)
                                              tACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCtt                                              <  2:6840292/151‑1 (MQ=9)
                                              tACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCtt                                              <  1:8269732/151‑1 (MQ=9)
                                              tACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCtt                                              <  1:11247739/151‑1 (MQ=9)
                                               aCAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtc                                             <  2:4473033/151‑1 (MQ=9)
                                                cAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcg                                            <  2:7605023/151‑1 (MQ=9)
                                                  gCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtc                                          >  2:6316744/1‑151 (MQ=9)
                                                   cGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcg                                         >  2:3543097/1‑151 (MQ=9)
                                                   cGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcg                                         >  1:1091771/1‑151 (MQ=9)
                                                    gTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                        <  1:10412150/151‑1 (MQ=9)
                                                     ttCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgtc                                       <  1:8665327/151‑1 (MQ=9)
                                                      tCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgtcg                                      >  1:7201719/1‑151 (MQ=9)
                                                        aaaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                    >  2:5332436/1‑151 (MQ=9)
                                                        aaaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                    >  2:6045659/1‑151 (MQ=9)
                                                        aaaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                    >  2:9861911/1‑151 (MQ=9)
                                                        aaaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                    <  1:6584668/151‑1 (MQ=9)
                                                         aaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                   <  2:9631457/151‑1 (MQ=9)
                                                         aaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                   >  1:2650743/1‑151 (MQ=9)
                                                         aaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                   <  1:458005/151‑1 (MQ=9)
                                                         aaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                   >  2:1549355/1‑151 (MQ=9)
                                                           ttCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCggg                                 >  1:2081047/1‑151 (MQ=9)
                                                              gCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgt                              >  2:5736950/1‑151 (MQ=9)
                                                                    aaGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTgcg                        >  1:787831/1‑151 (MQ=9)
                                                                     aGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTgcgc                       <  2:2261956/151‑1 (MQ=9)
                                                                      gCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCa                      >  1:8346223/1‑151 (MQ=9)
                                                                         ccTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATta                   <  2:7513814/151‑1 (MQ=9)
                                                                         ccTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATta                   <  1:4054086/151‑1 (MQ=9)
                                                                           tCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCCTTTACTCTTCGTCGTCGCCGGGAGTGGTGAGCATgata                 >  2:11105956/1‑151 (MQ=11)
                                                                           tCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATtata                 <  1:9913345/151‑1 (MQ=9)
                                                                             ggtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAgg               <  1:4351486/149‑1 (MQ=9)
                                                                              gtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAggg              <  1:9677537/150‑1 (MQ=14)
                                                                                 aaGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgaga           >  2:5866996/1‑151 (MQ=14)
                                                                                  atcatCGTTTAGAAACCCCATCGGGACCGAAACTCTGCGGCGAACAACCTGCCGGCGCGCTTCACCACTCCAGCGAGAGAATCCTTAAATCATTGCTTTTCAAGCATTTACTCTTCGTCGTCGCCGGGATTGGTGCGAATTATAGGGAGAt          <  2:2044161/147‑1 (MQ=11)
                                                                                    cTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAg        <  2:8933666/151‑1 (MQ=14)
                                                                                    cTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATGATAGGGAGATAg        >  2:10138968/1‑151 (MQ=9)
                                                                                     ttCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGa       <  1:5657358/151‑1 (MQ=14)
                                                                                        gTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACACCCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAAc    <  2:4326994/151‑1 (MQ=9)
                                                                                         tttCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACt   <  2:2003214/151‑1 (MQ=14)
                                                                                         tttCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACt   <  2:2839518/151‑1 (MQ=14)
                                                                                               ||                                                                                                                                                
TTGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGATGAAGCTTCGTTTCGAGCCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACTG  >  NC_008463/4957454‑4957694

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: