Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,953,687 T→A 5.7% intergenic (‑63/+19) PP_mr62 ← / ← trxA Pseudomon‑Rho regulatory region/thioredoxin I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,953,6870TA 5.7% 940.9 / 31.0 278intergenic (‑63/+19)PP_mr62/trxAPseudomon‑Rho regulatory region/thioredoxin I
Reads supporting (aligned to +/- strand):  ref base T (134/128);  new base A (10/6);  total (144/134)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

Only 100 of 279 total aligned reads displayed.

CTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCGGGG  >  NC_002947/5953543‑5953833
                                                                                                                                                 |                                                                                                                                                  
 tGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCgg                                                                                                                                              <  1:2499130/150‑1 (MQ=255)
       gCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGGCCCCGCAAATGCGGGGctt                                                                                                                                         <  2:2451076/150‑3 (MQ=255)
        cGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCttt                                                                                                                                       <  1:109636/150‑1 (MQ=255)
        cGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCttt                                                                                                                                       <  2:2632720/150‑1 (MQ=255)
        cGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCttt                                                                                                                                       <  1:4986447/150‑1 (MQ=255)
         gACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCtttt                                                                                                                                      >  2:3003063/1‑150 (MQ=255)
          aCAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAAGGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTtc                                                                                                                                     <  2:3458312/150‑1 (MQ=255)
            aGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGGCCCCGCAAATGCGGGGCttttctt                                                                                                                                    <  2:2190959/150‑4 (MQ=255)
            aGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTtcac                                                                                                                                   >  2:3519325/1‑150 (MQ=255)
             ggATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTtcaca                                                                                                                                  <  1:355647/150‑1 (MQ=255)
              gATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTtcacat                                                                                                                                 >  1:79273/1‑150 (MQ=255)
                tCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTtcacatca                                                                                                                               >  2:5356791/1‑150 (MQ=255)
                  gcgTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTtcacatcaca                                                                                                                             <  1:712862/150‑1 (MQ=255)
                     tACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGt                                                                                                                          <  2:2788617/150‑1 (MQ=255)
                     tACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGt                                                                                                                          <  2:3783656/150‑1 (MQ=255)
                     tACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCATTTCACATCACAGGt                                                                                                                          <  2:1747571/150‑1 (MQ=255)
                      aCTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTg                                                                                                                         <  1:2034814/150‑1 (MQ=255)
                      aCTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTg                                                                                                                         <  1:4267564/150‑1 (MQ=255)
                       cTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCg                                                                                                                                               <  1:4438690/127‑1 (MQ=255)
                       cTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCg                                                                                                                                               >  2:4438690/1‑127 (MQ=255)
                          aTTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCg                                                                                                                     >  1:4626682/1‑150 (MQ=255)
                            tCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTc                                                                                                                   <  1:2075228/150‑1 (MQ=255)
                              gAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGa                                                                                                                 <  2:5488428/150‑1 (MQ=255)
                              gAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCCTTTCACATCACAGGTGGGCGTCGa                                                                                                                 >  1:3807684/1‑150 (MQ=255)
                               aaTTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAg                                                                                                                <  2:2025148/150‑1 (MQ=255)
                               aaTTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAg                                                                                                                >  1:2977077/1‑150 (MQ=255)
                                     ggATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGc                                                                                                          >  2:3676674/1‑150 (MQ=255)
                                       aTGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCgg                                                                                                        >  1:3658784/1‑150 (MQ=255)
                                       aTGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCgg                                                                                                        <  1:2724680/150‑1 (MQ=255)
                                          cTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCa                                                                                                     >  2:1681576/1‑150 (MQ=255)
                                           ttCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAg                                                                                                    >  2:804290/1‑150 (MQ=255)
                                           ttCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAg                                                                                                    >  1:1129680/1‑150 (MQ=255)
                                           ttCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCATTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAg                                                                                                    <  1:4687753/150‑1 (MQ=255)
                                            tCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGc                                                                                                   >  2:3995469/1‑150 (MQ=255)
                                             cGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCt                                                                                                  <  1:1756642/150‑1 (MQ=255)
                                              ggCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTg                                                                                                 >  1:2941191/1‑150 (MQ=255)
                                                    ggCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGAttt                                                                                           <  1:162872/150‑1 (MQ=255)
                                                    ggCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGAttt                                                                                           <  1:25816/150‑1 (MQ=255)
                                                     gCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTg                                                                                          <  2:4536319/150‑1 (MQ=255)
                                                     gCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTg                                                                                          >  2:4901277/1‑150 (MQ=255)
                                                     gCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTg                                                                                          <  2:4382672/150‑1 (MQ=255)
                                                       acgTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGa                                                                                        <  2:3072873/148‑1 (MQ=255)
                                                         gTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACa                                                                                      >  2:3608557/1‑150 (MQ=255)
                                                          tAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAg                                                                                     >  2:2665970/1‑150 (MQ=255)
                                                            gagGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTg                                                                                   <  1:4028594/150‑1 (MQ=255)
                                                               ggCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTgcgc                                                                                >  2:39509/1‑150 (MQ=255)
                                                                  aCTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTgg                                                                                                                        <  2:5625720/107‑1 (MQ=255)
                                                                  aCTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTgg                                                                                                                        >  1:5625720/1‑107 (MQ=255)
                                                                       aaaaGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg                                                                        >  1:1965357/1‑150 (MQ=255)
                                                                         aaGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt                                                                      <  2:1654959/150‑1 (MQ=255)
                                                                         aaGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt                                                                      <  2:2489397/150‑1 (MQ=255)
                                                                         aaGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtggt                                                                      <  2:1141192/150‑1 (MQ=255)
                                                                                         aaTGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGGGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCg                                                      >  1:5528201/1‑150 (MQ=255)
                                                                                             tAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCt                                                  >  1:1298623/1‑150 (MQ=255)
                                                                                             tAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCATTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCt                                                  <  2:1602147/150‑1 (MQ=255)
                                                                                              aaCACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCtt                                                 <  1:2406417/150‑1 (MQ=255)
                                                                                              aaCACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCtt                                                 <  1:3467910/150‑1 (MQ=255)
                                                                                              aaCACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCtt                                                 <  1:3939046/150‑1 (MQ=255)
                                                                                               acacTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTtg                                                                        <  1:5083187/126‑1 (MQ=255)
                                                                                                 acTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGaa                                              <  2:2973817/150‑1 (MQ=255)
                                                                                                    tGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAg                                           >  2:357070/1‑150 (MQ=255)
                                                                                                      cTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCa                                         <  1:3047137/150‑1 (MQ=255)
                                                                                                       ttttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCCTTTGCGGGGCCTTTTCACCTCCCAGGTGGGCGTCGAGGTATGCGGCCAGCTGCGATTTGGCCGGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCAt                                        >  2:5422250/1‑150 (MQ=255)
                                                                                                       ttttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAAGCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCAt                                        >  1:1452468/1‑150 (MQ=255)
                                                                                                       ttttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCAt                                        >  2:390955/1‑150 (MQ=255)
                                                                                                       ttttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCAt                                        >  2:1741464/1‑150 (MQ=255)
                                                                                                       ttttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCAt                                        >  1:4193472/1‑150 (MQ=255)
                                                                                                        tttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATc                                       <  1:3386375/150‑1 (MQ=255)
                                                                                                         ttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCa                                      <  1:4935622/150‑1 (MQ=255)
                                                                                                         ttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCa                                      >  2:779712/1‑150 (MQ=255)
                                                                                                         ttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCa                                      <  2:1378039/150‑1 (MQ=255)
                                                                                                         ttttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCa                                      <  1:4612888/150‑1 (MQ=255)
                                                                                                           ttCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGc                                    <  2:5093889/150‑1 (MQ=255)
                                                                                                              tGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTGCTTGAACAGCATCAGCGTc                                 >  1:4148557/1‑150 (MQ=255)
                                                                                                               gCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCATTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCg                                <  1:2175919/150‑1 (MQ=255)
                                                                                                                       tGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCa                        >  2:3233193/1‑150 (MQ=255)
                                                                                                                          tctGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCCTCAGCGTCGGGATACCACGc                     >  1:4484789/1‑150 (MQ=255)
                                                                                                                          tctGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGc                     >  1:5225153/1‑150 (MQ=255)
                                                                                                                          tctGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGc                     >  1:3252244/1‑150 (MQ=255)
                                                                                                                            tGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGcac                   >  2:2734846/1‑150 (MQ=255)
                                                                                                                             gAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGcaca                  >  1:1789247/1‑150 (MQ=255)
                                                                                                                              aaaaaGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCacac                 >  2:5017992/1‑150 (MQ=255)
                                                                                                                               aaaagaaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGATCAGCATCAGCGTCGGGATACCACGCACAcc                <  1:3003063/150‑1 (MQ=255)
                                                                                                                               aaaagaaaagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACAcc                >  2:2265433/1‑150 (MQ=255)
                                                                                                                                gtgaaaaggCCCCGCAAATGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCg               >  1:5599887/5‑150 (MQ=255)
                                                                                                                                   aaaaggCCCCGCAAATGCGGGGCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCt           >  1:4329750/2‑150 (MQ=255)
                                                                                                                                    aaaggCCCCGCAAATGCGGGGCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCtt          >  2:106402/1‑150 (MQ=255)
                                                                                                                                      aagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTg         >  2:4333012/1‑150 (MQ=255)
                                                                                                                                      aagCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTg         >  1:4819179/1‑150 (MQ=255)
                                                                                                                                      aagCCCCGCATTTGCGGGGCCATTTCACATCACAGGAGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTAGCCGTGCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGAGCtgg         >  2:5238694/1‑148 (MQ=255)
                                                                                                                                       agCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTgg        <  1:3109232/150‑1 (MQ=255)
                                                                                                                                         gCCCCGCAAATGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGc       >  1:3333270/1‑150 (MQ=255)
                                                                                                                                          ccccGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCt                                                  >  1:3855449/1‑106 (MQ=255)
                                                                                                                                          ccccGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCt                                                  <  2:3855449/106‑1 (MQ=255)
                                                                                                                                          ccccGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGcc      >  2:4614633/1‑150 (MQ=255)
                                                                                                                                           cccGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCg     >  1:1022764/1‑150 (MQ=255)
                                                                                                                                            ccGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCt                                                  <  1:619625/104‑1 (MQ=255)
                                                                                                                                            ccGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCgg    >  1:2320824/1‑150 (MQ=255)
                                                                                                                                              gCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCgggg  <  2:4626744/150‑1 (MQ=255)
                                                                                                                                              gCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCgggg  <  2:2355909/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGATGCTTCGGCGACGGCATGTAGAGGGCACTGGAAAAGCAGTGCGAGGCCGAATGTAACACTTGCTTTTTTCTGCGTCTAGTGCTCTGAAAAAGAAAAGCCCCGCATTTGCGGGGCCTTTTCACATCACAGGTGGGCGTCGAGGAATGCGGCCAGCTGCGATTTGGACAGTGCGCCAACCTTGGTGGCTTCGACGTTGCCGTTCTTGAACAGCATCAGCGTCGGGATACCACGCACACCGTGCTTGGCCGGGG  >  NC_002947/5953543‑5953833

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: