Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,341,704 Δ1 bp pseudogene (664/957 nt) PA14_RS10890 → DUF58 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,341,7020C.99.1% 2280.6 / 8.4 470pseudogene (662/957 nt)PA14_RS10890DUF58 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base . (216/247);  total (220/250)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 487 total aligned reads displayed.

GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGG  >  NC_008463/2341550‑2341848
                                                                                                                                                        |                                                                                                                                                  
gAGCTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAc                                                                                                                                                    <  1:1044158/151‑1 (MQ=255)
 aGCTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTccc                                                                                                                                                   >  2:2938985/1‑151 (MQ=255)
   cTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGTCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGCCTACCAGGGGCTGCGCGCCTAccag                                                                                                                                                 >  1:10134677/1‑149 (MQ=255)
   cTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccag                                                                                                                                                 <  1:11121928/151‑3 (MQ=255)
    tGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccagc                                                                                                                                                >  1:9081939/1‑148 (MQ=255)
       gTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCgg                                                                                                                                            <  1:9820303/151‑1 (MQ=255)
        tCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCggg                                                                                                                                           <  1:2047733/151‑1 (MQ=255)
          gACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCCGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTCCCAGCCGGGCg                                                                                                                                         >  1:5829106/1‑151 (MQ=255)
          gACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCg                                                                                                                                         >  1:2763850/1‑151 (MQ=255)
            ccTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGAc                                                                                                                                       >  1:6431876/1‑151 (MQ=255)
               tgcgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTc                                                                                                                                     >  2:6517280/5‑151 (MQ=255)
               tgcgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTc                                                                                                                                     <  2:7874268/147‑1 (MQ=255)
                gcgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGCCTGCGCGCCTACCAGCCGGGCGACTcg                                                                                                                                    >  2:1016485/4‑151 (MQ=255)
                 cgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTcgc                                                                                                                                   <  1:8356583/149‑1 (MQ=255)
                  gccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTcgcg                                                                                                                                  >  2:4708549/2‑151 (MQ=255)
                  gccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTcgcg                                                                                                                                  <  1:8595223/150‑1 (MQ=255)
                   ccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcgg                                                                                                                                 >  2:1876638/1‑151 (MQ=255)
                    cagcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggc                                                                                                                                <  2:6016234/151‑1 (MQ=255)
                     agcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcg                                                                                                                               <  2:2505536/151‑1 (MQ=255)
                       cagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcggc                                                                                                                             >  2:11048956/1‑151 (MQ=255)
                        agGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcggcg                                                                                                                            >  2:7643090/1‑151 (MQ=255)
                              tGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCa                                                                                                                      >  2:3156190/1‑151 (MQ=255)
                               ggTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGTGGAGGGGCGGCGACGGCTACCAGGGACTCCGCGCCTCGCGGACGGGCGGCTCGAGGAGGCGCCTGGAc                                                                                                                     >  2:5629564/1‑151 (MQ=255)
                               ggTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCAc                                                                                                                     >  2:1380431/1‑151 (MQ=255)
                               ggTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCAc                                                                                                                     <  2:11346674/151‑1 (MQ=255)
                                   tACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGa                                                                                                                 >  1:5508031/1‑151 (MQ=255)
                                   tACCCTCGGCCGCTGGCCGGCGAGCCGGTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCCGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACACGCGGCGGCGCCTGCACTGGa                                                                                                                 <  2:6661604/151‑1 (MQ=255)
                                    aCCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACAGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGaa                                                                                                                >  1:1279461/1‑151 (MQ=255)
                                     cccTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAg                                                                                                               >  1:1576989/1‑151 (MQ=255)
                                     cccTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGGCGACGACTACCAGGGACTGCGCGCCTACCAGCCGTGCGACTCGCGGCGGCGCCTGCACTGGAAg                                                                                                               >  2:379463/1‑151 (MQ=255)
                                       cTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGc                                                                                                             >  2:7343212/1‑151 (MQ=255)
                                        tCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCg                                                                                                            >  1:5948855/1‑151 (MQ=255)
                                          ggCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGtt                                                                                                          <  1:635279/151‑1 (MQ=255)
                                          ggCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGtt                                                                                                          >  2:3482247/1‑151 (MQ=255)
                                          ggCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGtt                                                                                                          >  1:1854403/1‑151 (MQ=255)
                                             cGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTc                                                                                                       >  1:9918721/1‑151 (MQ=255)
                                              gCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGAGCGCGCTGGCGAGCACGTGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTTGAAGGCGTTTTcc                                                                                                      <  2:9227346/151‑1 (MQ=255)
                                               cTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTccc                                                                                                     <  2:10878332/151‑1 (MQ=255)
                                                tGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCg                                                                                                    <  1:964951/151‑1 (MQ=255)
                                                 ggCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGt                                                                                                   >  1:5118561/1‑151 (MQ=255)
                                                 ggCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGt                                                                                                   >  1:7104268/1‑151 (MQ=255)
                                                    cGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGc                                                                                                >  1:10756722/1‑151 (MQ=255)
                                                    cGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCGGGCGATCAAGGGGTTCGTAGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACCCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGc                                                                                                <  1:2435861/151‑1 (MQ=255)
                                                        cgagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAg                                                                                             <  2:2775814/147‑1 (MQ=255)
                                                        cgagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAg                                                                                             <  1:9486414/147‑1 (MQ=255)
                                                          agccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAggg                                                                                           >  2:10672073/3‑151 (MQ=255)
                                                              ggtCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgctg                                                                                       >  1:309705/3‑151 (MQ=255)
                                                              gTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgctg                                                                                       <  1:8102182/151‑1 (MQ=255)
                                                              gTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTGTTCCCGTGGCCAGGGgctg                                                                                       >  1:6641081/1‑151 (MQ=255)
                                                                 ccTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGctgctg                                                                                    >  1:9283163/1‑151 (MQ=255)
                                                                 ccTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGctgctg                                                                                    <  1:531692/151‑1 (MQ=255)
                                                                   tCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGGCTACCAGCCGGGCGACTCGCGGCGGCGCGTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGt                                                                                  <  2:1651020/151‑1 (MQ=255)
                                                                     gcgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGa                                                                                <  1:4476776/151‑1 (MQ=255)
                                                                          gggCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGag                                                                           >  1:7829889/1‑150 (MQ=255)
                                                                             cgcgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCCGGGGATGCTGGTGAAGGACTTc                                                                        >  2:8160953/1‑151 (MQ=255)
                                                                             cgcgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTc                                                                        <  1:11450123/151‑1 (MQ=255)
                                                                              gcgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcg                                                                       >  1:10556919/1‑151 (MQ=255)
                                                                                 gaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgccg                                                                    <  2:8667823/151‑1 (MQ=255)
                                                                                     aTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCt                                                                <  2:2589888/151‑1 (MQ=255)
                                                                                      tGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTc                                                               <  2:8618433/151‑1 (MQ=255)
                                                                                      tGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTc                                                               <  1:4659643/151‑1 (MQ=255)
                                                                                         cGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggc                                                            >  2:9580349/1‑151 (MQ=255)
                                                                                           aGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggcgg                                                          <  1:3575253/151‑1 (MQ=255)
                                                                                            ggACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCggg                                                         <  2:6941738/151‑1 (MQ=255)
                                                                                                   gggCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGAcc                                                  >  1:252244/1‑151 (MQ=255)
                                                                                                     gcgcgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCtg                                                <  2:3671125/151‑1 (MQ=255)
                                                                                                      cgcgcgcTGGCGAGCACGGGCGTCGTCGACTACCAGTGAATGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTCTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCtgt                                               <  1:2344755/151‑1 (MQ=255)
                                                                                                        cgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTgg                                             <  1:2897794/151‑1 (MQ=255)
                                                                                                         gcgcTGGCGAGCACGGGCGTCGCCGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGGGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGc                                            <  1:6782033/151‑1 (MQ=255)
                                                                                                          cgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCt                                           <  1:3850804/151‑1 (MQ=255)
                                                                                                             tGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGa                                        <  1:4365176/151‑1 (MQ=255)
                                                                                                             tGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCAGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGa                                        <  2:123218/151‑1 (MQ=255)
                                                                                                             gggCGAGCACGGGCGTCGACGACGACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGa                                        <  2:156073/150‑1 (MQ=255)
                                                                                                                 gAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTc                                    <  2:5166852/151‑1 (MQ=255)
                                                                                                                 gAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTc                                    >  2:7180683/1‑151 (MQ=255)
                                                                                                                   gCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGa                                  >  2:3371264/1‑151 (MQ=255)
                                                                                                                   gCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGa                                  <  1:7343212/151‑1 (MQ=255)
                                                                                                                    cACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAg                                 >  1:10272018/1‑151 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGACGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgc                                >  1:3330971/1‑151 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgc                                >  2:978844/1‑151 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgc                                >  1:10941249/1‑151 (MQ=255)
                                                                                                                      cGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgcg                               <  1:4266340/151‑1 (MQ=255)
                                                                                                                       gggCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgcgc                              <  1:606041/151‑1 (MQ=255)
                                                                                                                           gTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcg                          >  1:621443/1‑151 (MQ=255)
                                                                                                                           gTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcg                          >  2:4980434/1‑151 (MQ=255)
                                                                                                                           gTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcg                          <  1:7541215/151‑1 (MQ=255)
                                                                                                                           gTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcg                          >  2:10437162/1‑151 (MQ=255)
                                                                                                                           gTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcg                          >  2:4258278/1‑151 (MQ=255)
                                                                                                                               acgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcgg                      >  1:8997584/1‑151 (MQ=255)
                                                                                                                                cgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCggg                     <  2:8496313/151‑1 (MQ=255)
                                                                                                                                cgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCggg                     >  2:880811/1‑151 (MQ=255)
                                                                                                                                  acTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGa                   <  1:3908270/151‑1 (MQ=255)
                                                                                                                                      ccAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAtctc               <  2:4415034/151‑1 (MQ=255)
                                                                                                                                        aGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGa             <  2:2366030/151‑1 (MQ=255)
                                                                                                                                        aGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGa             >  2:8996361/1‑151 (MQ=255)
                                                                                                                                           gACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAg          <  2:2191859/151‑1 (MQ=255)
                                                                                                                                               gcgcgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGACAGTCGc      >  2:9360122/1‑151 (MQ=255)
                                                                                                                                                  cgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTg   <  2:6255594/151‑1 (MQ=255)
                                                                                                                                                  cgcCTACCAGCCGGGCGACTCGCGGCGGCGCCAGCACTGGCAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCACGAAAGTCGCCTg   <  2:5780286/151‑1 (MQ=255)
                                                                                                                                                   gcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTgg  <  1:8773855/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGG  >  NC_008463/2341550‑2341848

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: