Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,407,236 A→G intergenic (+107/‑205) PA14_RS20115 → / → toxA ATP‑grasp domain‑containing protein/exotoxin A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,407,2360AG99.8% 1568.8 / NA 395intergenic (+107/‑205)PA14_RS20115/toxAATP‑grasp domain‑containing protein/exotoxin A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (172/222);  minor base T (1/0);  total (173/222)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 398 total aligned reads displayed.

TGGCCCGCCGCACGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGAAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCTCGCCCGAC  >  NC_008463/4407088‑4407383
                                                                                                                                                    |                                                                                                                                                   
tGGCCCGCCGCACGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGgaa                                                                                                                                                   <  1:4500042/151‑1 (MQ=255)
tGGCCCGCCGCACGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGgaa                                                                                                                                                   <  1:6298224/151‑1 (MQ=255)
 ggCCCGCCGCACGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGACACTCGCCCAGAGCCTCGACCACTTACTCTAATCCGAAGCCCAGgaag                                                                                                                                                  >  2:7216178/1‑150 (MQ=255)
    ccgccgCACGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGt                                                                                                                                               >  2:10263351/1‑151 (MQ=255)
     cgccgcACGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGtt                                                                                                                                              <  1:615235/151‑1 (MQ=255)
        cgcACGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGACCGGTTgtc                                                                                                                                           >  1:11470915/1‑151 (MQ=255)
         gcACGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTgtcc                                                                                                                                          >  1:7887455/1‑151 (MQ=255)
          cACGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTgtccg                                                                                                                                         <  2:11532718/151‑1 (MQ=255)
           aCGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTgtccgt                                                                                                                                        <  1:7556278/151‑1 (MQ=255)
            cGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTgtccgtc                                                                                                                                       >  1:4104034/1‑151 (MQ=255)
             gCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTgtccgtcc                                                                                                                                      >  1:5753855/1‑151 (MQ=255)
              cTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTccc                                                                                                                                     >  2:10679505/1‑151 (MQ=255)
               ttCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTcccc                                                                                                                                    <  2:6305133/151‑1 (MQ=255)
                tCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCt                                                                                                                                   <  2:9677344/151‑1 (MQ=255)
                tCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGAACAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCGCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCAGTCCCCt                                                                                                                                   >  2:9435514/1‑151 (MQ=255)
                 cgcgGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTc                                                                                                                                  <  2:3037921/151‑1 (MQ=255)
                 cgcgGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTc                                                                                                                                  <  1:1853220/151‑1 (MQ=255)
                 cgcgGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTc                                                                                                                                  <  2:10339496/151‑1 (MQ=255)
                        aTGCTGGCACAGGCATAGTCGGACCCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATg                                                                                                                           >  1:3119214/1‑151 (MQ=255)
                         tGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATgg                                                                                                                          >  1:5013829/1‑151 (MQ=255)
                         tGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATgg                                                                                                                          <  1:5328242/151‑1 (MQ=255)
                           cTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGat                                                                                                                        <  1:2568079/151‑1 (MQ=255)
                           cTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGat                                                                                                                        <  2:10838594/151‑1 (MQ=255)
                             ggCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGatat                                                                                                                      <  2:6213294/151‑1 (MQ=255)
                                acaGGCATAGTCGGAACCGCTGCACTCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATttc                                                                                                                   <  1:1655972/151‑1 (MQ=255)
                                 caGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATttct                                                                                                                  <  1:8766057/151‑1 (MQ=255)
                                 caGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATttct                                                                                                                  <  1:5947262/151‑1 (MQ=255)
                                   ggCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATttcttc                                                                                                                >  1:8886319/1‑151 (MQ=255)
                                   ggCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATttcttc                                                                                                                >  1:5184357/1‑151 (MQ=255)
                                    gCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTcc                                                                                                               <  1:10981093/151‑1 (MQ=255)
                                       tAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCAct                                                                                                            <  1:2313016/151‑1 (MQ=255)
                                       tAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCAct                                                                                                            <  2:5924505/151‑1 (MQ=255)
                                        aGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCActc                                                                                                           >  1:8274011/1‑151 (MQ=255)
                                         gTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCActct                                                                                                          <  2:558067/151‑1 (MQ=255)
                                            ggAACCGCTGCAATCCCAACCCGACCCAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACGTACTCTCATCCGAAGCCCAGGAACGGGTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCa                                                                                                       >  2:2315107/1‑151 (MQ=255)
                                            ggAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCa                                                                                                       <  1:4226543/151‑1 (MQ=255)
                                            ggAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACACTGCa                                                                                                       >  1:2242648/1‑151 (MQ=255)
                                              aaCCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCaag                                                                                                     >  2:2567977/1‑150 (MQ=255)
                                               aCCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATc                                                                                                    <  1:8949614/151‑1 (MQ=255)
                                                  gCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCATTCCAg                                                                                                 >  1:3806945/1‑151 (MQ=255)
                                                  gCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAg                                                                                                 <  1:10677623/151‑1 (MQ=255)
                                                  gCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAg                                                                                                 <  2:4600495/151‑1 (MQ=255)
                                                  gCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAg                                                                                                 <  1:3436232/151‑1 (MQ=255)
                                                   cTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGt                                                                                                <  2:11575803/151‑1 (MQ=255)
                                                    tGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGtt                                                                                               >  2:2046937/1‑151 (MQ=255)
                                                    tGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGtt                                                                                               <  2:516451/151‑1 (MQ=255)
                                                    tGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGtt                                                                                               >  2:8000239/1‑151 (MQ=255)
                                                          cccAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATaaa                                                                                         <  2:2723924/151‑1 (MQ=255)
                                                           ccAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGGGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCagatat                                                                                        >  2:8742441/1‑146 (MQ=255)
                                                              aCCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATccc                                                                                     >  1:6963514/1‑151 (MQ=255)
                                                              aCCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACGTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCCTAAATccc                                                                                     <  2:613092/151‑1 (MQ=255)
                                                                ccGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCCGTTCATAAATCCCag                                                                                   >  2:7494895/1‑150 (MQ=255)
                                                                ccGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCAt                                                                                   >  1:3916582/1‑151 (MQ=255)
                                                                ccGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCAt                                                                                   >  1:6687030/1‑151 (MQ=255)
                                                                 cGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATa                                                                                  >  1:10133042/1‑151 (MQ=255)
                                                                 cGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATa                                                                                  <  1:4664857/151‑1 (MQ=255)
                                                                  gACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATaa                                                                                 <  2:6885425/151‑1 (MQ=255)
                                                                   aCCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATaaa                                                                                <  1:3513613/151‑1 (MQ=255)
                                                                       aGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGcc                                                                            >  2:7706572/1‑151 (MQ=255)
                                                                        gCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGccc                                                                           <  1:7431190/151‑1 (MQ=255)
                                                                           ccccTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCctct                                                                        >  1:3268849/1‑151 (MQ=255)
                                                                             ccTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTc                                                                      >  2:764481/1‑151 (MQ=255)
                                                                                tGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGc                                                                   >  1:5907503/1‑151 (MQ=255)
                                                                                 gatcgatcCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCt                                                                  <  2:4788136/151‑1 (MQ=255)
                                                                                 gatcgatcCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCt                                                                  <  1:926013/151‑1 (MQ=255)
                                                                                   tcgatcCCGATGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAATCCCAGGAACGGATGTCCGTCCCCTCGCCCATTGATATTTATTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTcc                                                                >  2:3357251/1‑151 (MQ=255)
                                                                                   tcgatcCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTcc                                                                >  1:7767670/1‑151 (MQ=255)
                                                                                     gatcCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTcccc                                                              <  2:3893161/151‑1 (MQ=255)
                                                                                         ccGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCa                                                          <  1:3687099/151‑1 (MQ=255)
                                                                                         ccGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCa                                                          <  1:11202720/151‑1 (MQ=255)
                                                                                         ccGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCa                                                          >  2:5789096/1‑151 (MQ=255)
                                                                                         ccGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCa                                                          <  1:844753/151‑1 (MQ=255)
                                                                                           gAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGc                                                        >  1:11116297/1‑151 (MQ=255)
                                                                                            cggggggCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGcc                                                       <  2:9282108/150‑1 (MQ=255)
                                                                                                      gTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATc                                             >  1:1990963/1‑151 (MQ=255)
                                                                                                         aCTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCg                                          >  2:1515199/1‑151 (MQ=255)
                                                                                                             gCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCAcc                                      <  2:10113848/151‑1 (MQ=255)
                                                                                                              cccAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCAccc                                     <  2:648864/151‑1 (MQ=255)
                                                                                                                    gCCTCGACCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACg                               <  1:2494365/151‑1 (MQ=255)
                                                                                                                          aCCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCccgc                         >  1:10252218/1‑151 (MQ=255)
                                                                                                                          aCCACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCccgc                         >  2:1341709/1‑151 (MQ=255)
                                                                                                                            cACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCccgccg                       >  2:3694091/1‑151 (MQ=255)
                                                                                                                            cACTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCccgccg                       >  1:980332/1‑151 (MQ=255)
                                                                                                                             aCTTACTCTCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCcgccgc                      >  2:1110557/1‑151 (MQ=255)
                                                                                                                                  ctctCATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTcc                 >  2:8438798/1‑151 (MQ=255)
                                                                                                                                   tctcATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTccg                <  1:9722147/151‑1 (MQ=255)
                                                                                                                                     taaTCCGAAGCCGAGGAACGGTTGTCCGTACCCAAGCCCATCGATATTTCTTCCACTCTGCAATCCAGGTCATAAAACCAATAAAATCCCTCTACCGCTCCCCGCCATCCTCCCCGCATCCCGCACCCTAGACGCACCGCAGCTCAccgcc              >  2:897583/3‑151 (MQ=255)
                                                                                                                                      cATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTccgccg             <  2:4293092/151‑1 (MQ=255)
                                                                                                                                      cATCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTccgccg             <  1:9940744/151‑1 (MQ=255)
                                                                                                                                       aTCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCgg            >  2:5741712/1‑151 (MQ=255)
                                                                                                                                       aTCCGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCgg            >  1:8582805/1‑151 (MQ=255)
                                                                                                                                         ccGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCt          >  1:4750287/1‑151 (MQ=255)
                                                                                                                                          cGAAGCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCTc         <  2:179193/151‑1 (MQ=255)
                                                                                                                                          cGAAGCACAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTACTTCCACTCTGCAATCCAGTTCATAAAGCCCATAATAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCTc         <  2:7955538/151‑1 (MQ=255)
                                                                                                                                              gCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCCTAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCTCGccc     <  1:9165780/151‑1 (MQ=255)
                                                                                                                                              gCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCTCGccc     <  1:5690505/151‑1 (MQ=255)
                                                                                                                                              gCCCAGGAACGGTTGTCCGTCCCCTCGCCCATGGAGATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCTCGccc     <  2:2129215/151‑1 (MQ=255)
                                                                                                                                               cccAGGAACGGTAGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCAAAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCCCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCTCGCCCg    <  1:2840487/151‑1 (MQ=255)
                                                                                                                                                 cAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCTCGCCCGAc  <  1:7306582/151‑1 (MQ=255)
                                                                                                                                                 cAGGAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCTCGCCCGAc  <  1:1479899/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGCCCGCCGCACGCTTCGCGGCGATGCTGGCACAGGCATAGTCGGAACCGCTGCAATCCCAACCCGACCAAGCTCCCCTTGATCGATCCCGAGGGGGGCTCGTCACTCGCCCAGTGCCTCGACCACTTACTCTCATCCGAAGCCCAGAAACGGTTGTCCGTCCCCTCGCCCATGGATATTTCTTCCACTCTGCAATCCAGTTCATAAATCCCATAAAAGCCCTCTTCCGCTCCCCGCCAGCCTCCCCGCATCCCGCACCCTAGACGCCCCGCCGCTCTCCGCCGGCTCGCCCGAC  >  NC_008463/4407088‑4407383

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: