Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,566,788 A→G pseudogene (1139/1275 nt) PA14_RS16215 → fatty acid desaturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,566,7880AG99.8% 1881.5 / NA 474pseudogene (1139/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (223/250);  minor base T (0/1);  total (223/251)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 475 total aligned reads displayed.

CCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGA  >  NC_008463/3566639‑3566934
                                                                                                                                                     |                                                                                                                                                   
  tGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacca                                                                                                                                                 >  1:4672013/1‑151 (MQ=255)
   gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccac                                                                                                                                                >  2:10378137/1‑151 (MQ=255)
    cccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccacc                                                                                                                                               <  1:10057377/151‑1 (MQ=255)
    cccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccacc                                                                                                                                               >  2:2365385/1‑151 (MQ=255)
      cAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCCTGTGGGAGACCTTCCTGATCCTGCGCGACCACctc                                                                                                                                             >  2:6133943/1‑151 (MQ=255)
       aGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACctct                                                                                                                                            >  1:5837363/1‑151 (MQ=255)
         ccACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGc                                                                                                                                          <  2:11185760/151‑1 (MQ=255)
         ccACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGc                                                                                                                                          >  1:7340979/1‑151 (MQ=255)
          cACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGcc                                                                                                                                         <  1:1018234/151‑1 (MQ=255)
          cACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGcc                                                                                                                                         <  2:7308706/151‑1 (MQ=255)
          cACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGcc                                                                                                                                         <  2:10629115/151‑1 (MQ=255)
             gACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcgcg                                                                                                                                      >  2:5371315/1‑151 (MQ=255)
              aCTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCgg                                                                                                                                     <  2:9870647/151‑1 (MQ=255)
                ttCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGcc                                                                                                                                   <  1:6158042/151‑1 (MQ=255)
                 tCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa                                                                                                                                  >  1:110179/1‑151 (MQ=255)
                 tCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa                                                                                                                                  <  1:6725088/151‑1 (MQ=255)
                  ccaccaCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAg                                                                                                                                 <  2:1547467/151‑1 (MQ=255)
                   caccacCGCTAGACGGGGTTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCTAAGTTCGGCCCGATGACCGTGACCTGGGGAATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCAGCGGCCAGt                                                                                                                                <  2:4550531/151‑1 (MQ=255)
                      ccaccGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCg                                                                                                                              <  2:7125140/148‑1 (MQ=255)
                      ccaccGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCg                                                                                                                              <  1:5228317/148‑1 (MQ=255)
                          cGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGAcc                                                                                                                          <  1:6383487/151‑1 (MQ=255)
                           gCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGAccc                                                                                                                         <  1:539329/151‑1 (MQ=255)
                            cAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCg                                                                                                                        <  2:11038123/151‑1 (MQ=255)
                             aaGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGt                                                                                                                       <  1:1878996/151‑1 (MQ=255)
                              aGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGtt                                                                                                                      <  1:8410726/151‑1 (MQ=255)
                               gCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTc                                                                                                                     <  1:835530/151‑1 (MQ=255)
                                ccGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCg                                                                                                                    <  1:8173977/151‑1 (MQ=255)
                                 cGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCgg                                                                                                                   <  1:5052430/151‑1 (MQ=255)
                                  gggggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGa                                                                                                                 <  2:6714584/151‑1 (MQ=255)
                                     ggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaatc                                                                                                              <  1:2758134/151‑4 (MQ=255)
                                     ggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaatc                                                                                                              >  2:388469/1‑148 (MQ=255)
                                       tGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaatcgt                                                                                                            <  1:113239/151‑6 (MQ=255)
                                        gAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAGCGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtc                                                                                                            >  1:10292444/1‑151 (MQ=255)
                                        gAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtc                                                                                                            <  1:2396227/151‑1 (MQ=255)
                                        gAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtc                                                                                                            >  2:8989434/1‑151 (MQ=255)
                                          aCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGa                                                                                                          >  1:4267636/1‑151 (MQ=255)
                                            ttCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACt                                                                                                        >  1:10683560/1‑151 (MQ=255)
                                             tCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTg                                                                                                       >  2:7487098/1‑151 (MQ=255)
                                                  cgcATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTAcc                                                                                                  <  2:3744185/151‑1 (MQ=255)
                                                       cgcgCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCgag                                                                                             >  1:349850/1‑151 (MQ=255)
                                                          gcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                                                          >  2:7160334/1‑151 (MQ=255)
                                                          gcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                                                          <  2:2099219/151‑1 (MQ=255)
                                                             ggAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCAt                                                                                       <  1:9492738/151‑1 (MQ=255)
                                                              gAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATg                                                                                      >  2:5477968/1‑151 (MQ=255)
                                                               aGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGc                                                                                     <  2:3424335/151‑1 (MQ=255)
                                                                  gcgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCg                                                                                  <  2:2616619/151‑1 (MQ=255)
                                                                   cgcgGCCAACGAGGCGCGCCCCTCGACGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCgg                                                                                 <  2:4430472/151‑1 (MQ=255)
                                                                   cgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCgg                                                                                 <  1:3727543/151‑1 (MQ=255)
                                                                    gcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGATCTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgc                                                                                <  2:3657189/151‑1 (MQ=255)
                                                                      ggCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctc                                                                              >  1:328722/1‑151 (MQ=255)
                                                                       gCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcg                                                                             <  1:9126464/151‑1 (MQ=255)
                                                                         caacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgct                                                                           <  2:7987635/151‑1 (MQ=255)
                                                                         caacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgct                                                                           <  2:3199485/151‑1 (MQ=255)
                                                                          aacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgctc                                                                          <  2:9699686/151‑1 (MQ=255)
                                                                           acgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTcc                                                                         <  1:3608530/148‑1 (MQ=255)
                                                                            cgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGGGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGGCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCg                                                                        >  1:7669642/3‑151 (MQ=255)
                                                                              aGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGga                                                                      >  2:491788/1‑151 (MQ=255)
                                                                               ggCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGgag                                                                     <  2:1626377/151‑1 (MQ=255)
                                                                               ggCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGgag                                                                     >  2:10785057/1‑151 (MQ=255)
                                                                                      cccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATg                                                              <  2:1319377/151‑1 (MQ=255)
                                                                                        cTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgcc                                                            <  2:5219892/151‑1 (MQ=255)
                                                                                         tcgaaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccg                                                           <  2:1782167/148‑1 (MQ=255)
                                                                                           gaaGTTCGGCCCGATGACCGAGAACTGGGGCATGTGGGATAACTACCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGATCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccgcc                                                         >  2:890492/2‑151 (MQ=255)
                                                                                               ttCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAcc                                                     <  2:8675754/151‑1 (MQ=255)
                                                                                               ttCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAAcc                                                     <  1:1949089/151‑1 (MQ=255)
                                                                                                tCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCccacca                                                    >  1:9969569/1‑151 (MQ=255)
                                                                                                acGGCCCGATGACCGAGACATGGGGCATGAGGGAGAACTACCTGATCCTGCGCGACCACCTCAGCCGCGGCCACACGGACCCGTACGGAAACGTCCACTGGTACCGAGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCa                                                    >  2:8325393/2‑151 (MQ=255)
                                                                                                 cGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTTATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGACCGGCTCGCTCCGGAGACCGATGCCGCCAACCAg                                                   <  2:5793896/151‑1 (MQ=255)
                                                                                                 cGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTTGTACCGCGACAACCAAGCCCGGCTCGCTCCGGAGGCGTAAGCCGCCAACCac                                                   >  2:4509678/1‑150 (MQ=255)
                                                                                                 cGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAg                                                   >  1:1717814/1‑151 (MQ=255)
                                                                                                 cGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAg                                                   >  2:10513049/1‑151 (MQ=255)
                                                                                                  ggCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGc                                                  >  1:2987795/1‑151 (MQ=255)
                                                                                                     ccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGcccc                                               >  1:6591688/1‑151 (MQ=255)
                                                                                                     ccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGcccc                                               >  2:5084946/1‑151 (MQ=255)
                                                                                                       gATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCcac                                             >  1:4906419/1‑151 (MQ=255)
                                                                                                       gATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCcac                                             <  1:4190126/151‑1 (MQ=255)
                                                                                                              gagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaacc                                      <  2:6943151/148‑1 (MQ=255)
                                                                                                              gagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaacc                                      >  1:8034731/4‑151 (MQ=255)
                                                                                                              gagaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaacc                                      <  2:9195535/148‑1 (MQ=255)
                                                                                                                gaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGga                                                                      <  1:7822098/116‑1 (MQ=255)
                                                                                                                gaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGga                                                                      >  2:7822098/2‑117 (MQ=255)
                                                                                                                gaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACAccaaccaa                                    >  1:6998395/2‑151 (MQ=255)
                                                                                                                gaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACAccaaccaa                                    <  1:5699735/150‑1 (MQ=255)
                                                                                                                 aCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACcaaccaac                                   >  1:3395668/1‑151 (MQ=255)
                                                                                                                 aCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAACATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACcaaccaac                                   <  2:5140845/151‑1 (MQ=255)
                                                                                                                  ccTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCaaccaacc                                  >  1:9690711/1‑151 (MQ=255)
                                                                                                                       ggCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCAt                             <  2:8012823/151‑1 (MQ=255)
                                                                                                                         cATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATcc                           <  1:1477055/151‑1 (MQ=255)
                                                                                                                         cATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCAGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCTCTCCGGAGACCGATGCCGCCAACCATCCCCACACCAACCAACCCGCATcc                           >  2:8728355/1‑151 (MQ=255)
                                                                                                                             tgGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAgg                       >  2:6226826/1‑151 (MQ=255)
                                                                                                                               ggAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCt                     <  1:2197836/151‑1 (MQ=255)
                                                                                                                                   aaCTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCg                 <  2:8161153/151‑1 (MQ=255)
                                                                                                                                   aaCTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCg                 >  2:4682395/1‑151 (MQ=255)
                                                                                                                                           gATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATca         >  1:3451142/1‑151 (MQ=255)
                                                                                                                                            aTCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcac        <  1:8963035/151‑1 (MQ=255)
                                                                                                                                             tCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCCACCAACCCGCATCCCAGGCTACCGCCTGATcaca       <  2:984074/151‑1 (MQ=255)
                                                                                                                                             tCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcaca       >  2:4269566/1‑151 (MQ=255)
                                                                                                                                             tCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcaca       >  2:6277089/1‑151 (MQ=255)
                                                                                                                                                tGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAg    <  1:8772789/151‑1 (MQ=255)
                                                                                                                                                 gcgcGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAgg   <  1:4536521/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGA  >  NC_008463/3566639‑3566934

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: