Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,947,641 2 bp→GC coding (455‑456/1416 nt) PA14_RS22670 → DNA polymerase Y family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,947,6410CG99.2% 1009.1 / ‑4.2 264T152R (ACA→AGA) PA14_RS22670DNA polymerase Y family protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base G (139/123);  total (141/123)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084634,947,6420AC99.2% 1024.1 / ‑1.4 264coding (456/1416 nt)PA14_RS22670DNA polymerase Y family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (139/123);  minor base . (2/0);  total (141/123)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

CCTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAAGGTCATGCGCACTACACAGGACAACATTTCCATCGAGGACGAGACAGCGCTGCAACACCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTGC  >  NC_008463/4947522‑4947760
                                                                                                                         ||                                                                                                                      
ccTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGAc                                                                                           <  1:7169172/150‑1 (MQ=255)
 cTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACa                                                                                          >  2:1700956/1‑150 (MQ=255)
  tGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAg                                                                                         >  1:10327928/1‑150 (MQ=255)
  tGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAg                                                                                         <  2:6672776/150‑1 (MQ=255)
  tGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAg                                                                                         <  1:8108616/150‑1 (MQ=255)
    ttCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAgcg                                                                                       >  2:2317671/1‑150 (MQ=255)
     tCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAgcgc                                                                                      <  1:5996284/150‑1 (MQ=255)
      cGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCt                                                                                     >  2:9845575/1‑150 (MQ=255)
       ggCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTg                                                                                    >  2:2600982/1‑150 (MQ=255)
       ggCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTg                                                                                    <  2:3204186/150‑1 (MQ=255)
       ggCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTg                                                                                    >  2:6883148/1‑150 (MQ=255)
        gCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGc                                                                                   <  1:3701094/150‑1 (MQ=255)
         ccccTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCa                                                                                  <  1:3440930/150‑1 (MQ=255)
           ccTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAac                                                                                >  1:2657942/1‑150 (MQ=255)
           ccTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAac                                                                                <  1:8155935/150‑1 (MQ=255)
             tGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAacgc                                                                              <  1:5035803/150‑3 (MQ=255)
              ggCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGcc                                                                             >  1:6919683/1‑150 (MQ=255)
              ggCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGcc                                                                             <  1:7938275/150‑1 (MQ=255)
              ggCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGcc                                                                             <  2:9295107/150‑1 (MQ=255)
              ggCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGcc                                                                             <  1:1583630/150‑1 (MQ=255)
               gCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCa                                                                            >  1:9317878/1‑150 (MQ=255)
                ccACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAg                                                                           >  2:10531220/1‑150 (MQ=255)
                ccACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAg                                                                           <  1:1019726/150‑1 (MQ=255)
                   cGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTg                                                                        <  2:3366327/150‑1 (MQ=255)
                      ttCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTcc                                                                     >  2:1472787/1‑150 (MQ=255)
                      ttCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTcc                                                                     <  2:5549949/150‑1 (MQ=255)
                          aGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCGTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGccc                                                                 >  2:5319674/1‑150 (MQ=255)
                           ggCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCt                                                                >  2:492420/1‑150 (MQ=255)
                           ggCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCt                                                                <  1:4211574/150‑1 (MQ=255)
                           ggCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCt                                                                >  2:9239715/1‑150 (MQ=255)
                           ggCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCt                                                                >  1:4887378/1‑150 (MQ=255)
                            gCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCtt                                                               >  2:3902489/1‑150 (MQ=255)
                            gCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCtt                                                               >  2:5506412/1‑150 (MQ=255)
                             cccGCTTGCGCAAGGAACTCACGGCGCTGGGCTTGCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTc                                                              >  1:9880927/1‑150 (MQ=255)
                             cccGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTc                                                              >  2:8201710/1‑150 (MQ=255)
                             cccGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTc                                                              >  2:2940096/1‑150 (MQ=255)
                             cccGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTc                                                              >  1:6555518/1‑150 (MQ=255)
                               cGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCg                                                            <  1:7154693/150‑1 (MQ=255)
                               cGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCg                                                            >  2:5985539/1‑150 (MQ=255)
                                    gcgcAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGa                                                       >  2:272761/1‑150 (MQ=255)
                                    gcgcAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGa                                                       >  2:2087646/1‑150 (MQ=255)
                                      gcAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGAcc                                                     >  1:8399663/1‑150 (MQ=255)
                                      gcAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCAGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGAcc                                                     <  2:9414701/150‑1 (MQ=255)
                                       cAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCg                                                    >  1:10465527/1‑150 (MQ=255)
                                        aaGGAACTCACGGCGCTGGGCTACCGCCCTCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGCGGGTCCGCGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGc                                                   <  2:6746742/150‑1 (MQ=255)
                                          ggAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCt                                                 >  1:10396653/1‑150 (MQ=255)
                                           gAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCACCGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTg                                                >  1:5602381/1‑150 (MQ=255)
                                             aCTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGACCAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCACCGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTggg                                              >  2:6186605/1‑150 (MQ=255)
                                              cTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGa                                             <  1:996878/150‑1 (MQ=255)
                                               tCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGAc                                            <  2:4253761/150‑1 (MQ=255)
                                                  cGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGActct                                         <  1:3083897/150‑1 (MQ=255)
                                                  cGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGActct                                         <  1:10631461/150‑1 (MQ=255)
                                                  cGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGActct                                         <  2:2116488/150‑1 (MQ=255)
                                                   ggCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACtctc                                        >  1:725427/1‑150 (MQ=255)
                                                   ggCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACtctc                                        <  2:858332/150‑1 (MQ=255)
                                                   ggCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACtctc                                        <  2:1438895/150‑1 (MQ=255)
                                                    gcgcTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTcg                                       >  2:813926/1‑150 (MQ=255)
                                                    gcgcTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTcg                                       <  2:3548148/150‑1 (MQ=255)
                                                      gcTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTcgcg                                     <  2:9635398/150‑1 (MQ=255)
                                                       cTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCgg                                    <  1:3209387/150‑1 (MQ=255)
                                                         gggCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGa                                  <  2:4823706/150‑1 (MQ=255)
                                                              tCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTgg                             >  2:4411513/1‑150 (MQ=255)
                                                              tCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTgg                             >  1:8998305/1‑150 (MQ=255)
                                                               ccGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGc                            <  1:2333382/150‑1 (MQ=255)
                                                               ccGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGc                            <  2:1306929/150‑1 (MQ=255)
                                                               ccGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGc                            <  1:7716352/150‑1 (MQ=255)
                                                                cGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGcc                           >  2:279662/1‑150 (MQ=255)
                                                                 gCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCg                          <  1:4247094/150‑1 (MQ=255)
                                                                  ccATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGt                         <  1:443434/150‑1 (MQ=255)
                                                                  ccATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGt                         >  2:250495/1‑150 (MQ=255)
                                                                  ccATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGt                         <  1:6416618/150‑1 (MQ=255)
                                                                   cATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTc                        >  2:1988390/1‑150 (MQ=255)
                                                                    aTCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCg                       <  1:5265626/150‑1 (MQ=255)
                                                                    aTCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCg                       >  1:5781852/1‑150 (MQ=255)
                                                                      cgcctgcgcctgGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGcc                     <  2:2249568/150‑1 (MQ=255)
                                                                      cGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGcc                     <  2:2252460/150‑1 (MQ=255)
                                                                       gCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGccc                    <  2:5717566/150‑1 (MQ=255)
                                                                        ccTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCt                   <  2:2208938/150‑1 (MQ=255)
                                                                        ccTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCt                   <  2:9091680/150‑1 (MQ=255)
                                                                         cTGTGCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTc                  <  2:1175389/150‑1 (MQ=255)
                                                                          tgtgCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCg                 >  2:7452380/1‑150 (MQ=255)
                                                                          tgtgCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGAAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCg                 >  2:6547075/1‑150 (MQ=255)
                                                                           gtgCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCgg                >  2:3781104/1‑150 (MQ=255)
                                                                            tgCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGTCCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGa               >  2:2270963/1‑150 (MQ=255)
                                                                             gCCTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGAc              >  2:3806475/1‑150 (MQ=255)
                                                                               cTGGCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACTGCCACGACGGCGGGGCCATCGGGGACGAGACAGCGCTGCAACGCCGGTGGTCCGCCCTGCGGGTGGACCGGCTGGGGGGCTCGCGGGTGGGGGCCGGCGGCCGCggcggc            >  2:7098564/1‑146 (MQ=255)
                                                                                 ggCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCACCGCCAGTTGTCCGCCCTTCCGCTGGACCGCCGGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCAt          >  2:888467/1‑150 (MQ=255)
                                                                                 ggCGCCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCAt          <  2:307205/150‑1 (MQ=255)
                                                                                  gcgcCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATg         <  1:1849742/150‑1 (MQ=255)
                                                                                    gcCGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATggg       >  2:6264994/1‑150 (MQ=255)
                                                                                      cGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGAc     >  2:7585234/1‑150 (MQ=255)
                                                                                      cGAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGAc     >  2:4640919/1‑150 (MQ=255)
                                                                                       gAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACt    <  1:9108624/150‑1 (MQ=255)
                                                                                       gAACCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACt    <  2:5097871/150‑1 (MQ=255)
                                                                                        aaCCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTg   >  2:8395303/1‑150 (MQ=255)
                                                                                        aaCCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTg   >  2:10387096/1‑150 (MQ=255)
                                                                                         aCCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGTGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTgc  >  2:3510224/1‑150 (MQ=255)
                                                                                         aCCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTgc  >  2:1546412/1‑150 (MQ=255)
                                                                                         aCCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTgc  >  2:953644/1‑150 (MQ=255)
                                                                                         aCCCGGTGGCGGCGAGGGTCCTGGCGAACAGCCACGACGGCCTGGCCATCGAGGACGAGACAGCGCTGCAACGCCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTgc  >  2:10168321/1‑150 (MQ=255)
                                                                                                                         ||                                                                                                                      
CCTGTTCGGCCCCTGGCCACGCTTCGAGGCCCGCTTGCGCAAGGAACTCACGGCGCTGGGCTTCCGCCATCGCCTGTGCCTGGCGCCGAACCCGGTGGCGGCGAAGGTCATGCGCACTACACAGGACAACATTTCCATCGAGGACGAGACAGCGCTGCAACACCAGTTGTCCGCCCTTCCGCTGGACCGCCTGGGACTCTCGCGGGAGCTGGCCGTCGCCCTCGGACGCATGGGACTGC  >  NC_008463/4947522‑4947760

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: