Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,281,477 G→C intergenic (‑25/+44) katE ← / ← PA14_RS14985 catalase HPII/general stress protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,281,4770GC99.5% 1417.5 / ‑2.0 368intergenic (‑25/+44)katE/PA14_RS14985catalase HPII/general stress protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base C (159/206);  total (162/206)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 373 total aligned reads displayed.

GGCGTCCTCGCGATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGGGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTGCCGCCGCTGGCCTGGCCGCCCTTCTTG  >  NC_008463/3281330‑3281623
                                                                                                                                                   |                                                                                                                                                  
ggCGTCCTCGCGATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATggggcgg                                                                                                                                                  <  1:9534079/150‑1 (MQ=255)
ggCGTCCTCGCGATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCAACGGTTTACGCTCCTCATCCTGGGCATGAAAATggggcgg                                                                                                                                                  >  2:7728369/1‑150 (MQ=255)
 gCGTCCTCGCGATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGTGGCGGa                                                                                                                                                 >  2:5163111/1‑150 (MQ=255)
 gCGTCCTCGCGATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGa                                                                                                                                                 <  1:8851074/150‑1 (MQ=255)
  cGTCCTCGCGATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGAt                                                                                                                                                <  2:1884883/150‑1 (MQ=255)
   gTCCTCGCGATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATg                                                                                                                                               <  1:4635742/150‑1 (MQ=255)
       tCGCGATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTGGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTc                                                                                                                                           >  2:1415657/1‑150 (MQ=255)
         gcgATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTccc                                                                                                                                         <  2:2750311/150‑1 (MQ=255)
         gcgATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTccc                                                                                                                                         <  2:807320/150‑1 (MQ=255)
          cgATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCt                                                                                                                                        <  2:9852289/150‑1 (MQ=255)
             tAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCt                                                                                                                                     <  1:4230366/150‑1 (MQ=255)
             tAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCt                                                                                                                                     >  2:5300097/1‑150 (MQ=255)
              aGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTc                                                                                                                                    <  1:9781927/150‑1 (MQ=255)
               ggCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTcc                                                                                                                                   <  2:1089733/150‑1 (MQ=255)
               ggCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTcc                                                                                                                                   <  1:10461680/150‑1 (MQ=255)
               ggCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTcc                                                                                                                                   >  1:7098804/1‑150 (MQ=255)
                  cTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGcc                                                                                                                                >  2:3348290/1‑150 (MQ=255)
                   tCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGccc                                                                                                                               >  1:588702/1‑150 (MQ=255)
                    ccAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGcccc                                                                                                                              >  1:3280270/1‑150 (MQ=255)
                     cAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCa                                                                                                                             >  2:10216004/1‑150 (MQ=255)
                          gTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCAttcac                                                                                                                        <  2:8277038/150‑1 (MQ=255)
                           ttCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCCTCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCAttcact                                                                                                                       <  2:7834386/150‑1 (MQ=255)
                           ttCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCAttcact                                                                                                                       >  1:1610094/1‑150 (MQ=255)
                            tCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCAttcactt                                                                                                                      <  2:4656511/150‑1 (MQ=255)
                               aGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCAttcacttcac                                                                                                                   <  1:1118298/150‑1 (MQ=255)
                               aGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCAttcacttcac                                                                                                                   <  2:2351177/150‑1 (MQ=255)
                               aGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCAttcacttcac                                                                                                                   >  1:409734/1‑150 (MQ=255)
                                  ttGGCTTTGCTGTTCCCACGGTCGACGTAGACGGTGACGGCAGCATGTCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGTTCATGTAAATGGGGCGTAAGCCTCCCTCCTCCGCCCCATACACTTCACTgc                                                                                                                >  2:6823694/1‑150 (MQ=255)
                                   tGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTgcg                                                                                                               <  2:1262117/150‑1 (MQ=255)
                                     gCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGcgcg                                                                                                             <  2:7256789/150‑1 (MQ=255)
                                      cGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCgg                                                                                                            <  1:5470190/150‑1 (MQ=255)
                                       ggtGCTGGTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGt                                                                                                           <  2:10172261/148‑1 (MQ=255)
                                       gTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGt                                                                                                           <  1:1704076/150‑1 (MQ=255)
                                          gCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCAGTCACTTCACTGCGCGGTTAt                                                                                                        >  2:3045802/1‑150 (MQ=255)
                                           cTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATc                                                                                                       <  2:2704790/150‑1 (MQ=255)
                                           cTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATc                                                                                                       <  1:6714312/150‑1 (MQ=255)
                                           cTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATc                                                                                                       <  1:7980797/150‑1 (MQ=255)
                                           cTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATc                                                                                                       >  1:390303/1‑150 (MQ=255)
                                           cTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATc                                                                                                       <  1:3036474/150‑1 (MQ=255)
                                            tGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCa                                                                                                      <  2:6323752/150‑1 (MQ=255)
                                            tGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCa                                                                                                      <  1:3428821/150‑1 (MQ=255)
                                             gTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTTCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCAAGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAg                                                                                                     >  2:9183839/1‑150 (MQ=255)
                                               tGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCCGtt                                                                                                   <  2:4377303/150‑1 (MQ=255)
                                               tGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGGTTCCTCTCCACCTCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGtt                                                                                                   <  2:7623643/150‑1 (MQ=255)
                                                 ccACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCCTCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTcc                                                                                                 <  2:8314830/150‑1 (MQ=255)
                                                  cACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTccc                                                                                                >  1:5084369/1‑150 (MQ=255)
                                                     ggTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCccgcc                                                                                             <  2:9102837/150‑1 (MQ=255)
                                                       tCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGt                                                                                           <  2:3765837/150‑1 (MQ=255)
                                                       tCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGt                                                                                           <  2:1325011/150‑1 (MQ=255)
                                                       tCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGt                                                                                           <  1:8279806/150‑1 (MQ=255)
                                                         gACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTgg                                                                                         >  1:9908026/1‑150 (MQ=255)
                                                          aCGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGc                                                                                        <  1:9156670/150‑1 (MQ=255)
                                                          aCGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGc                                                                                        <  2:4365229/150‑1 (MQ=255)
                                                            ggTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTg                                                                                      <  2:4152928/150‑1 (MQ=255)
                                                             gtgtCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTgc                                                                                     >  2:4063199/1‑150 (MQ=255)
                                                             gtgtCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTgc                                                                                     <  2:4895683/150‑1 (MQ=255)
                                                               gtCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTgcgc                                                                                   >  2:3135481/1‑150 (MQ=255)
                                                                tCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCt                                                                                  >  1:1423702/1‑150 (MQ=255)
                                                                 cGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCCTCGGTTTCCTCTCCTCAACCTGGGCATGAAAATGGGGCGGAGGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTGCCCGCCGTGGCTGCGCTg                                                                                 <  2:2294309/150‑1 (MQ=255)
                                                                  ggTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTgg                                                                                >  1:8384785/1‑150 (MQ=255)
                                                                   gTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGc                                                                               >  2:7159922/1‑150 (MQ=255)
                                                                   gTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGc                                                                               <  1:9024264/150‑1 (MQ=255)
                                                                      ccggccgCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCAc                                                                            <  1:8521255/150‑1 (MQ=255)
                                                                       cggccgCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCAcc                                                                           >  2:3854534/1‑150 (MQ=255)
                                                                               tGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGc                                                                   >  2:9727392/1‑150 (MQ=255)
                                                                               tGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGc                                                                   >  1:10206815/1‑150 (MQ=255)
                                                                               tGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGc                                                                   >  1:1982737/1‑150 (MQ=255)
                                                                                  ctgcgctgcTCGTTGTTCTGCTCGCTCATCGGTGTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAg                                                                <  2:8871293/150‑1 (MQ=255)
                                                                                    gcgctgcTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGcc                                                              <  2:450691/150‑1 (MQ=255)
                                                                                     cgctgcTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCt                                                             >  1:10650868/1‑150 (MQ=255)
                                                                                     cgctgcTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCt                                                             >  1:7845541/1‑150 (MQ=255)
                                                                                      gctgctCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTc                                                            <  2:8450111/150‑1 (MQ=255)
                                                                                           tCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAgg                                                       >  1:7381050/1‑150 (MQ=255)
                                                                                             gttgttCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGcc                                                     >  1:6137012/1‑150 (MQ=255)
                                                                                              ttgttCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGccc                                                    <  1:4591319/150‑1 (MQ=255)
                                                                                               tgttCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCg                                                   <  1:6502379/150‑1 (MQ=255)
                                                                                                gttCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGc                                                  <  2:1762521/150‑1 (MQ=255)
                                                                                                gttCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGc                                                  <  2:5724483/150‑1 (MQ=255)
                                                                                                 ttCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCt                                                 >  1:2927590/1‑150 (MQ=255)
                                                                                                     gctcgctcATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCgg                                             >  1:1560108/1‑150 (MQ=255)
                                                                                                       tcgctcATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGAt                                           <  1:2242260/150‑1 (MQ=255)
                                                                                                          ctcATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGt                                        >  2:6193456/1‑150 (MQ=255)
                                                                                                            cATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTc                                      <  2:1154660/150‑1 (MQ=255)
                                                                                                            cATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTc                                      >  1:1205958/1‑150 (MQ=255)
                                                                                                            cATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTc                                      <  1:2773842/150‑1 (MQ=255)
                                                                                                            cATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTc                                      >  1:559985/1‑150 (MQ=255)
                                                                                                                ggTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGa                                  >  1:10077510/1‑150 (MQ=255)
                                                                                                                 gTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGaa                                 <  1:9964397/150‑1 (MQ=255)
                                                                                                                    tcctctcctgATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGtt                              <  2:3126213/150‑1 (MQ=255)
                                                                                                                      ctctcctcATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTgc                            <  1:3568627/150‑1 (MQ=255)
                                                                                                                         tcctcATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTgccgc                         <  1:9383645/150‑1 (MQ=255)
                                                                                                                          cctcATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTgccgcc                        >  2:5589999/1‑150 (MQ=255)
                                                                                                                            tcATCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTGCcgccgc                      >  1:2607123/1‑150 (MQ=255)
                                                                                                                              aTCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTGCCGCCGctg                    <  2:5835555/150‑1 (MQ=255)
                                                                                                                              aTCCTGGGCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTGCCGCCGctg                    >  2:3529034/1‑150 (MQ=255)
                                                                                                                                ccTGGGCATGAAAATGGGGGGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTGCCGCCGctggc                  >  2:2071803/1‑150 (MQ=255)
                                                                                                                                    ggCATGTAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCTCTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTGCAGCCGctggcctgg              <  2:8201384/150‑1 (MQ=255)
                                                                                                                                     gCATGAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTGCCGCCGctggcctggc             >  2:10610467/1‑150 (MQ=255)
                                                                                                                                         gAAAATGGGGCGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTGCCGCCGCTGGCCTGgccgcc         <  2:5173934/150‑1 (MQ=255)
                                                                                                                                                gggcggATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTGCCGCCGCTGGCCTGGCCGCCCTTCTTg  >  2:5952091/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCGTCCTCGCGATAGGCCTCCAGCTGTTCCAGCTTGGCGTTGCTGTTGCCACGGTCGACGGTGTCGGTGCCGGCCGCCTGGCTGCGCTGCTCGTTGTTCTGCTCGCTCATCGGTTTCCTCTCCTCATCCTGGGCATGAAAATGGGGGGGATGCCTCCCTGCTCCGCCCCATTCACTTCACTGCGCGGTTATCAGTTCCCGCCGTGGCTGCGCTGGCCACCTTTCTTGCCAGCCTCGGAGGCCCGCTGCGGATCGTTCTTGAAGTTGCCGCCGCTGGCCTGGCCGCCCTTCTTG  >  NC_008463/3281330‑3281623

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: