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,5490GA100.0% 413.8 / NA 256intergenic (‑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);  new base A (123/133);  total (123/133)
*NC_0084634,957,5500CT99.6% 419.6 / NA 256intergenic (‑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 (122/133);  minor base A (1/0);  total (123/133)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

TTGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGATGAAGCTTCGTTTCGAGCCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACTG  >  NC_008463/4957454‑4957694
                                                                                               ||                                                                                                                                                
   cgcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCt                                                                                         >  2:1885344/1‑151 (MQ=9)
    gcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTg                                                                                        <  2:8734083/151‑1 (MQ=9)
     ccGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTCCCGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTgg                                                                                       <  2:6081012/151‑1 (MQ=255)
     ccGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTgg                                                                                       <  1:4960374/151‑1 (MQ=9)
     ccGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTgg                                                                                       >  1:1122551/1‑151 (MQ=9)
      cGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTggg                                                                                      <  2:5971643/151‑1 (MQ=9)
       gACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGga                                                                                     >  1:1886154/1‑151 (MQ=9)
       gACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACAGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGga                                                                                     >  1:124101/1‑151 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    <  1:10398923/151‑1 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    <  1:5775788/151‑1 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    <  1:5221879/151‑1 (MQ=9)
         ccTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgaga                                                                                   >  1:5090881/1‑151 (MQ=9)
          cTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagag                                                                                  >  1:1512621/1‑151 (MQ=9)
            tCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGaa                                                                                >  2:6971907/1‑151 (MQ=9)
              caaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAACTTCGCGGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATc                                                                              <  2:3166265/150‑1 (MQ=11)
              aaaCCATCTTCCGATCAGCGGGATGCGCATTACACATCGGTCAAATTCGGAGTCAAGCACCTCGGTGCAGCTTCGTTTCGTATACCGTTCGGGTCCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATc                                                                              <  2:5362490/151‑1 (MQ=11)
              aaaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATc                                                                              <  1:9804069/151‑1 (MQ=9)
               aaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATcc                                                                             <  1:5439915/151‑1 (MQ=9)
               aaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATcc                                                                             <  1:5640323/151‑1 (MQ=9)
                aCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCg                                                                            >  1:4359340/1‑151 (MQ=9)
                aCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCg                                                                            <  1:5377584/151‑1 (MQ=9)
                 ccATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGt                                                                           <  1:7392221/151‑1 (MQ=9)
                  cATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTa                                                                          >  1:404730/1‑151 (MQ=9)
                  cATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTa                                                                          >  2:4649803/1‑151 (MQ=9)
                    tCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTaaa                                                                        <  2:9465152/151‑1 (MQ=9)
                    tCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTaaa                                                                        <  2:3255734/151‑1 (MQ=9)
                     cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAt                                                                       <  1:2224882/151‑1 (MQ=9)
                      ttCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATc                                                                      >  1:6038726/1‑151 (MQ=9)
                        ccGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCt                                                                    >  1:31101/1‑151 (MQ=9)
                         cGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCtt                                                                   <  2:5073767/151‑1 (MQ=9)
                          gATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTg                                                                  >  1:1780546/1‑151 (MQ=9)
                          gATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTg                                                                  >  1:826593/1‑151 (MQ=9)
                          gATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTg                                                                  <  1:408641/151‑1 (MQ=9)
                           aTCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGa                                                                 >  1:6879234/1‑151 (MQ=9)
                           aTCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGa                                                                 >  1:5514924/1‑151 (MQ=9)
                            tCAGCGGGAGGCGCATTCTACAGCGTTCCAATTCGCTGGCAAGCCCCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGACCAACCTGCCGGAGCGCTTCTCGACTCCTGGGAGAGTATCCGTTAATCCTTGAt                                                                >  2:1416555/1‑151 (MQ=11)
                            tCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAt                                                                <  1:775976/151‑1 (MQ=9)
                            tCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAt                                                                <  2:355556/151‑1 (MQ=9)
                             cAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtt                                                               >  1:10620370/1‑151 (MQ=9)
                              aGCGGGAGGCGCATTCTACAGCGTTCAAAGTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTGCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAttt                                                              <  1:2580969/151‑1 (MQ=255)
                                 gggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCa                                                           >  2:2769020/1‑151 (MQ=9)
                                 gggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCa                                                           >  1:2354346/1‑151 (MQ=9)
                                  ggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCaa                                                          >  1:5728247/1‑151 (MQ=9)
                                  ggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCaa                                                          >  2:10489008/1‑151 (MQ=9)
                                    aGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGc                                                        <  1:10706286/151‑1 (MQ=9)
                                       cgcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAtt                                                     <  2:8690307/151‑1 (MQ=9)
                                        gcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGCGAGAATCCGTAAATCCTTGATTTTCAAGCAttt                                                    >  2:6415641/1‑151 (MQ=9)
                                        gcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAttt                                                    >  2:338557/1‑151 (MQ=9)
                                         cATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTa                                                   <  1:8546012/151‑1 (MQ=9)
                                         cATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTCAAGCATTTAc                                                  <  1:4508143/151‑1 (MQ=255)
                                          aTTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTAc                                                  >  2:7395847/1‑151 (MQ=9)
                                           ttCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTAct                                                 >  2:2410338/1‑151 (MQ=9)
                                           ttCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTAct                                                 <  1:6941835/151‑1 (MQ=9)
                                           ttCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTAct                                                 >  2:6758917/1‑151 (MQ=9)
                                              tACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCtt                                              <  2:4065760/151‑1 (MQ=9)
                                               aCAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtc                                             <  1:7256909/151‑1 (MQ=9)
                                                cAGCGTTCACATTCGCTGTCAAGCACCTCGGTGAAGCTTCGATTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGACCAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGAAAATCCTTGATTTTCAAGCATTTACTCTtcg                                            <  1:6302067/151‑1 (MQ=11)
                                                cAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcg                                            <  1:6476569/151‑1 (MQ=9)
                                                cAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcg                                            >  1:6698810/1‑151 (MQ=9)
                                                 aGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTgcgt                                           >  2:1140623/1‑151 (MQ=9)
                                                  gCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtc                                          >  2:7027030/1‑151 (MQ=9)
                                                   cGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcg                                         <  2:1685681/151‑1 (MQ=9)
                                                   cGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcg                                         >  2:778038/1‑151 (MQ=9)
                                                   cGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcg                                         >  1:8612999/1‑151 (MQ=9)
                                                    gTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                        >  1:8884828/1‑151 (MQ=9)
                                                    gTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                        >  2:3133749/1‑151 (MQ=9)
                                                     ttCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgtc                                       <  1:8993692/151‑1 (MQ=9)
                                                      tCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgtcg                                      >  2:67939/1‑151 (MQ=9)
                                                        aaaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                    >  2:10159438/1‑151 (MQ=9)
                                                          aTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCgg                                  >  2:4512843/1‑151 (MQ=9)
                                                            tCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGa                                <  1:147818/151‑1 (MQ=9)
                                                            tCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGa                                >  2:2097058/1‑151 (MQ=9)
                                                            tCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTGTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGa                                <  2:4821610/151‑1 (MQ=9)
                                                              gCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgt                              >  1:4892493/1‑151 (MQ=9)
                                                              gCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgt                              <  2:9295618/151‑1 (MQ=9)
                                                               cTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtg                             >  1:9344464/1‑151 (MQ=9)
                                                                 gTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtggt                           >  2:9352428/1‑151 (MQ=9)
                                                                 gTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTGCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAACTCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtggt                           <  2:10116044/151‑1 (MQ=255)
                                                                  tCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtggtg                          <  2:4968465/151‑1 (MQ=9)
                                                                  tCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtggtg                          <  1:8286980/151‑1 (MQ=9)
                                                                   cAAGCACCTCGGTGATGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCAATCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTgc                         <  2:850105/151‑1 (MQ=255)
                                                                   cAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTgc                         <  1:3493314/151‑1 (MQ=9)
                                                                       cACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCAt                     <  1:1187174/151‑1 (MQ=9)
                                                                        aCCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCAtt                    >  2:4195133/1‑151 (MQ=9)
                                                                        aCCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCAtt                    >  2:1000496/1‑151 (MQ=9)
                                                                        aCCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCAtt                    <  2:8182223/151‑1 (MQ=9)
                                                                          cTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATtat                  <  1:9140566/151‑1 (MQ=9)
                                                                               tGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAGTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGga             >  2:10741649/1‑151 (MQ=9)
                                                                               tGAAGCTTCGTTTCGAATCACTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGTGCGATTCACCTCTCCTGGGAGAGAATCAGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGATTGGTGCGCATTATATGga             <  2:7418405/151‑1 (MQ=11)
                                                                                gAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgag            >  2:4715785/1‑151 (MQ=14)
                                                                                gAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgag            <  1:10451433/151‑1 (MQ=14)
                                                                                gAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgag            <  1:10069167/151‑1 (MQ=14)
                                                                                 aaGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgaga           <  1:10237624/151‑1 (MQ=14)
                                                                                  aGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGAt          >  1:4556161/1‑151 (MQ=14)
                                                                                   gCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATa         <  2:2630547/151‑1 (MQ=14)
                                                                                     ttCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGa       >  2:1482435/1‑151 (MQ=14)
                                                                                     ttCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATGATAGGGAGATAGa       >  2:5695248/1‑151 (MQ=9)
                                                                                      tCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGaa      <  1:7090652/151‑1 (MQ=14)
                                                                                          ttCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACTg  <  1:9440384/151‑1 (MQ=14)
                                                                                          ttCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGGCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACTg  <  2:4033733/151‑1 (MQ=9)
                                                                                               ||                                                                                                                                                
TTGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGATGAAGCTTCGTTTCGAGCCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACTG  >  NC_008463/4957454‑4957694

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