Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,653,350 C→A 6.8% intergenic (‑185/+28) PP_4959 ← / ← fba diguanylate cyclase/fructose‑bisphosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,653,3500CA 6.8% 644.5 / 21.8 196intergenic (‑185/+28)PP_4959/fbadiguanylate cyclase/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (103/79);  new base A (8/6);  total (111/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

GATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCA  >  NC_002947/5653206‑5653497
                                                                                                                                                |                                                                                                                                                   
gATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCgg                                                                                                                                                <  2:533931/150‑1 (MQ=255)
gATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCgg                                                                                                                                                <  2:1387627/150‑1 (MQ=255)
    cATTAGTACTCTAGACCTGAGGCTGCAACACGCGATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGttt                                                                                                                                            <  2:5724/150‑1 (MQ=255)
    cATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGttt                                                                                                                                            <  1:1159990/150‑1 (MQ=255)
      ttAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGttttt                                                                                                                                          >  1:4544957/1‑150 (MQ=255)
      ttAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGttttt                                                                                                                                          <  2:3137293/150‑1 (MQ=255)
        aGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCa                                                                                                                                        <  2:1968681/150‑1 (MQ=255)
        aGTACTCGAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCa                                                                                                                                        <  1:2513569/150‑1 (MQ=255)
         gTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAg                                                                                                                                       >  1:5358236/1‑150 (MQ=255)
          tACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAgg                                                                                                                                      <  2:890576/150‑1 (MQ=255)
          tACTCTAGACCAGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAgg                                                                                                                                      <  2:2050511/150‑1 (MQ=255)
             tctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGtt                                                                                                                                   >  1:4918700/1‑150 (MQ=255)
              ctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTa                                                                                                                                  >  1:3431789/1‑150 (MQ=255)
               tAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTAc                                                                                                                                 <  1:983400/150‑1 (MQ=255)
               tAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTAc                                                                                                                                 <  2:125698/150‑1 (MQ=255)
                aGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTAcc                                                                                                                                >  1:1830210/1‑150 (MQ=255)
                   ccagAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGt                                                                                                                             <  1:3731265/147‑1 (MQ=255)
                    cTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTc                                                                                                                            >  1:5270163/1‑150 (MQ=255)
                    cTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTc                                                                                                                            >  2:1635861/1‑150 (MQ=255)
                     tGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCg                                                                                                                           >  2:1822901/1‑150 (MQ=255)
                           tGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGtttttt                                                                                                                                         >  1:3917308/1‑129 (MQ=255)
                           tGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGtttttt                                                                                                                                         <  2:3917308/130‑2 (MQ=255)
                                cacGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGAcc                                                                                                                <  2:4865920/150‑1 (MQ=255)
                                cacGCAATGCAATCAGAACACTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGAcc                                                                                                                >  1:3175658/1‑150 (MQ=255)
                                   gcaatgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTg                                                                                                             >  2:4261174/1‑150 (MQ=255)
                                   gcaatgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTg                                                                                                             >  2:4290916/1‑150 (MQ=255)
                                    caatgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTgg                                                                                                            <  1:628503/150‑1 (MQ=255)
                                     aatgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTggc                                                                                                           <  1:3024018/150‑1 (MQ=255)
                                     aatgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTggc                                                                                                           <  2:2340534/150‑1 (MQ=255)
                                           atCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAg                                                                                                     >  1:3256965/1‑150 (MQ=255)
                                            tCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAg                                                                                                     >  2:4436252/1‑149 (MQ=255)
                                            tCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGt                                                                                                    >  2:3452914/1‑150 (MQ=255)
                                              aGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTATCAGGGTTACCTGTCGATCAGTTGACCTAGGCGGCCAGTTc                                                                                                  >  2:3004627/1‑150 (MQ=255)
                                               gAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCa                                                                                                 >  1:2956644/1‑150 (MQ=255)
                                                      aTCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCa                                                                                          <  1:2689265/150‑1 (MQ=255)
                                                          aTGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTa                                                                                      >  1:2634044/1‑150 (MQ=255)
                                                          aTGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTa                                                                                      >  1:5359578/1‑150 (MQ=255)
                                                            gATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAgc                                                                                    >  2:843948/1‑150 (MQ=255)
                                                            gATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAgc                                                                                    >  2:843998/1‑150 (MQ=255)
                                                              tCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAgcgc                                                                                  >  2:2507005/1‑150 (MQ=255)
                                                              tCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAgcgc                                                                                  <  1:354275/150‑1 (MQ=255)
                                                                 aaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTgg                                                                               >  1:5419264/1‑150 (MQ=255)
                                                                    tGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTAc                                                                            >  2:1554435/1‑150 (MQ=255)
                                                                     gTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTg                                                                                                             >  1:1127227/1‑116 (MQ=255)
                                                                     gTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTg                                                                                                             <  2:1127227/116‑1 (MQ=255)
                                                                     gTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCAGCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGCTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACa                                                                           <  1:2231285/150‑1 (MQ=255)
                                                                          gTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCt                                                                      >  2:3791809/1‑150 (MQ=255)
                                                                            aGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTc                                                                    <  1:216174/150‑1 (MQ=255)
                                                                            aGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTc                                                                    <  1:1697112/150‑1 (MQ=255)
                                                                              ccGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAATCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACGTTTCAGGTAGCGCTGGTACATGCCTTCCa                                                                  >  2:4747525/1‑150 (MQ=255)
                                                                              ccGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCa                                                                  >  2:284172/1‑150 (MQ=255)
                                                                              ccGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCa                                                                  <  2:1218734/150‑1 (MQ=255)
                                                                                gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                                >  2:2243871/1‑150 (MQ=255)
                                                                                gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                                >  1:3699838/1‑150 (MQ=255)
                                                                                gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                                >  2:1793625/1‑150 (MQ=255)
                                                                                gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                                <  2:1568147/150‑1 (MQ=255)
                                                                                gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                                <  2:4828019/150‑1 (MQ=255)
                                                                                 aCTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCagag                                                               >  2:4407993/1‑150 (MQ=255)
                                                                                        gtcggtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCTGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCgg                                                        >  2:2501902/1‑150 (MQ=255)
                                                                                        gtcggtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCgg                                                        >  1:879537/1‑150 (MQ=255)
                                                                                         tcggtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGc                                                       <  1:3144758/150‑1 (MQ=255)
                                                                                            gtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTg                                                    >  1:4932663/1‑150 (MQ=255)
                                                                                               tGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATc                                                 >  1:1575038/1‑150 (MQ=255)
                                                                                                gCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCt                                                >  1:702680/1‑150 (MQ=255)
                                                                                                 ccGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCtt                                               >  2:4113272/1‑150 (MQ=255)
                                                                                                 ccGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCtt                                               >  2:4112513/1‑150 (MQ=255)
                                                                                                  cGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTc                                              <  2:2566429/150‑1 (MQ=255)
                                                                                                    gAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGa                                            <  1:582876/150‑1 (MQ=255)
                                                                                                    gAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGa                                            <  1:5425550/150‑1 (MQ=255)
                                                                                                    gAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGa                                            <  2:2533221/150‑1 (MQ=255)
                                                                                                    gAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGa                                            <  2:3309361/150‑1 (MQ=255)
                                                                                                      aTTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAgg                                          <  1:2338231/150‑1 (MQ=255)
                                                                                                        tATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTCCCAGAGAGATCGGCTTGATCTTCGAGGCa                                        <  2:4936224/150‑1 (MQ=255)
                                                                                                         aTGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCAt                                       <  2:3583268/150‑1 (MQ=255)
                                                                                                           gTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTg                                     >  1:469672/1‑150 (MQ=255)
                                                                                                            tACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGc                                    >  2:5102860/1‑150 (MQ=255)
                                                                                                            tACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGc                                    >  1:1939471/1‑150 (MQ=255)
                                                                                                             aCAGGTTTCCATCTCCCATAAAAAAACACGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCAGGTACATGCCTTCAAGAGAGATAGGCTTGATCTTCTAGGCATTGcc                                   >  1:4262383/1‑150 (MQ=255)
                                                                                                              cAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGccc                                  >  1:4792635/1‑150 (MQ=255)
                                                                                                              cAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGccc                                  >  1:1461450/1‑150 (MQ=255)
                                                                                                                 gTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCAtt                                      >  2:1469255/1‑143 (MQ=255)
                                                                                                                  tttCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCgg                              <  2:514104/150‑1 (MQ=255)
                                                                                                                   ttCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGt                             >  2:4732053/1‑150 (MQ=255)
                                                                                                                         ctcCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGaa                       <  1:3502444/150‑1 (MQ=255)
                                                                                                                           cccATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAgg                     <  2:4544957/150‑1 (MQ=255)
                                                                                                                            ccATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGc                    >  2:1923328/1‑150 (MQ=255)
                                                                                                                              aTAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCtt                  >  1:164625/1‑150 (MQ=255)
                                                                                                                              aTAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCtt                  >  2:2716822/1‑150 (MQ=255)
                                                                                                                              aTAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTAGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCtt                  >  1:3466922/1‑150 (MQ=255)
                                                                                                                               tAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTc                 <  2:127045/150‑1 (MQ=255)
                                                                                                                               ctgaaaaaCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTc                 >  1:1529054/4‑150 (MQ=255)
                                                                                                                                aaaaaaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCg                <  1:212083/150‑1 (MQ=255)
                                                                                                                                 gaaaaaCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCa                                                                                          >  1:442676/2‑75 (MQ=255)
                                                                                                                                   aaaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCt                                                                                              >  2:281982/1‑69 (MQ=255)
                                                                                                                                   aaaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCt                                                                                              <  1:281982/69‑1 (MQ=255)
                                                                                                                                    aaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAc            <  1:1722189/150‑1 (MQ=255)
                                                                                                                                        ccGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAcgcgc        >  1:3989724/1‑150 (MQ=255)
                                                                                                                                          gCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGa      >  2:2092544/1‑150 (MQ=255)
                                                                                                                                             tCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGc   >  1:4946014/1‑150 (MQ=255)
                                                                                                                                                cTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAcgcgc        >  2:3969356/1‑142 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCA  >  NC_002947/5653206‑5653497

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: ‑