Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,485,090 T→G 7.7% intergenic (+33/+161) moaE → / ← rhlB molybdopterin synthase large subunit/ATP‑dependent RNA helicase RhlB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,485,0900TG 7.7% 812.1 / 7.3 297intergenic (+33/+161)moaE/rhlBmolybdopterin synthase large subunit/ATP‑dependent RNA helicase RhlB
Reads supporting (aligned to +/- strand):  ref base T (114/151);  new base G (30/2);  total (144/153)
Fisher's exact test for biased strand distribution p-value = 1.52e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

TTCATCATGGACTACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCCACAGGGCGCGTGTAG  >  NC_002947/1484941‑1485239
                                                                                                                                                     |                                                                                                                                                     
ttCATCATGGACTACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGt                                                                                                                                                       >  2:3727648/1‑150 (MQ=255)
ttCATCATGGACTACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGt                                                                                                                                                       >  2:1363294/1‑150 (MQ=255)
ttCATCATGGACTACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGt                                                                                                                                                       >  2:7215023/1‑150 (MQ=255)
  catcatGGACTACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGttt                                                                                                                                                     >  2:2274875/1‑150 (MQ=255)
    tcatGGACTACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTAc                                                                                                                                                   <  1:5718457/150‑1 (MQ=255)
    tcatGGACTACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAAGACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTAc                                                                                                                                                   <  2:4777939/150‑1 (MQ=255)
      atGGACTACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTAccc                                                                                                                                                 >  1:100358/1‑150 (MQ=255)
        ggACTACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCcgc                                                                                                                                               <  2:919357/150‑1 (MQ=255)
            tACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTGTACCCGCGaaa                                                                                                                                           >  2:4748214/1‑150 (MQ=255)
              ccTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAgg                                                                                                                                         <  2:5842662/150‑1 (MQ=255)
               cTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTTGGAGCGGGTTTACCCGCGCAAggg                                                                                                                                        >  1:1273709/1‑150 (MQ=255)
               cTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGGTTTCCCGCGAAAggg                                                                                                                                        >  1:6884417/1‑150 (MQ=255)
                tGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGggc                                                                                                                                       <  2:6879325/150‑1 (MQ=255)
                 gAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGggcc                                                                                                                                      >  2:4185869/1‑150 (MQ=255)
                 gAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGggcc                                                                                                                                      <  1:4641083/150‑1 (MQ=255)
                  aaGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGCAAGggccg                                                                                                                                     >  2:7111311/1‑150 (MQ=255)
                     aCCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCACGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGggccggcc                                                                                                                                  <  1:7412976/150‑1 (MQ=255)
                     aCCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGggccggcc                                                                                                                                  >  2:2197966/1‑150 (MQ=255)
                      cccGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGccc                                                                                                                                 >  2:7747356/1‑150 (MQ=255)
                       ccGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCct                                                                                                                                >  2:3678420/1‑150 (MQ=255)
                       ccGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGGTTACCCGCGAAAGGGCCGGCCct                                                                                                                                >  2:73824/1‑150 (MQ=255)
                           gcgcCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCctgctg                                                                                                                            <  1:8690987/150‑1 (MQ=255)
                            cgcCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTgg                                                                                                                           <  2:5107299/150‑1 (MQ=255)
                             gcCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCGGCTGGc                                                                                                                          <  1:2786755/150‑1 (MQ=255)
                                 ttCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGCAAGGGCCGGCCCGGCTGGCCaaa                                                                                                                      >  1:8591773/1‑149 (MQ=255)
                                    tGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTACTCAAGGCGCCACCGGTCGCCGGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCGGGCCAAactc                                                                                                                   >  2:5569425/1‑150 (MQ=255)
                                     ggAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTcc                                                                                                                  >  2:2820886/1‑150 (MQ=255)
                                      gaagaaGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCa                                                                                                                 >  2:660354/1‑150 (MQ=255)
                                      gaagaaGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGACGGTTTAACCGCGCAAAGGCCCGCCCGGCAGGCCaaacaaaa                                                                                                                 >  1:5224794/1‑144 (MQ=255)
                                       aagaagGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCAc                                                                                                                <  2:8632161/150‑1 (MQ=255)
                                       aagaagGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCAc                                                                                                                <  2:1094732/150‑1 (MQ=255)
                                         gaagGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGGTTGCCCGCGAAAGGGCTGGCCCCGCGGGCCAAGCAcccccc                                                                                                              >  1:8493714/1‑150 (MQ=255)
                                            ggAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACccgcc                                                                                                           <  1:2258410/150‑1 (MQ=255)
                                             gagaATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACccgccg                                                                                                          >  2:7053698/1‑150 (MQ=255)
                                             gagaATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGGTTACCCGCGAAAGGGCCGGCCCCGCAGGCCCAACTCCAccccccc                                                                                                          >  1:7654450/1‑149 (MQ=39)
                                               gaATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAATGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCtt                                                                                                        >  1:4228067/1‑150 (MQ=255)
                                                  tACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCAc                                                                                                     <  2:1578953/150‑1 (MQ=255)
                                                   aCTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCAcc                                                                                                    <  1:5889861/150‑1 (MQ=255)
                                                    cTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCg                                                                                                   >  2:2447925/1‑150 (MQ=255)
                                                     tCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGCAAGGGCCGGCCCGGCTGGCCAATCTCCCCCCGCCGTTCACCGt                                                                                                  >  1:4978460/1‑150 (MQ=255)
                                                     tCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGt                                                                                                  >  2:8036608/1‑150 (MQ=255)
                                                     tCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGt                                                                                                  >  2:5683932/1‑150 (MQ=255)
                                                      cAGGACGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGtt                                                                                                 <  1:3349156/150‑1 (MQ=255)
                                                      cAGGAAGGGCCAAGGTGGTTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGGTTAACCGCGAAAGGGCCCGCCCTGCTGGCCAATCTCCACCCGCCGGTCACCgat                                                                                                 >  2:1851742/1‑148 (MQ=38)
                                                         gAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGTGCCGGGTTAACCGCGCAAGGGCCGGCCCAGCTGGCCAAACTCCACCCCCCCTTCAcccgggcc                                                                                              >  1:5958334/1‑144 (MQ=255)
                                                         gAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTgcc                                                                                              >  2:3353857/1‑150 (MQ=255)
                                                         gAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTgcc                                                                                              <  2:5348239/150‑1 (MQ=255)
                                                         gAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGGTTTCCCGCGGAAGGGCCCGCCCGGCTGGCCCAAAACCACCCGCCCGTCACCGATgcc                                                                                              >  1:3040982/1‑150 (MQ=255)
                                                             ggCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCa                                                                                          <  2:2612827/150‑1 (MQ=255)
                                                               ccAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGCAAGGGCCGGCCCTGCTGGCCAATCGCCACCCGCCGTTACCCGTTgccgccccg                                                                                        >  1:4321365/1‑147 (MQ=37)
                                                                cAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCAcaa                                                                                       >  1:3577787/1‑148 (MQ=255)
                                                                  aGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGGTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAAACTCCACCCCCCCGTCACCGTTGCCCCCACAAcc                                                                                     >  1:638906/1‑150 (MQ=37)
                                                                           gAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTcgcg                                                                            <  2:8762009/150‑1 (MQ=255)
                                                                            aGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCgg                                                                           <  1:7747356/150‑1 (MQ=255)
                                                                                 cggcAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGCAAGGGCCGGCCCGGCTGGCCACTCTCCACCCGCCGTTGACCGGTGCCGCCCCAACCCAATTCGCGGGTaaa                                                                      >  1:5339486/1‑150 (MQ=25)
                                                                                 cggcAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTaaa                                                                      >  2:1014645/1‑150 (MQ=255)
                                                                                  ggcAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAAc                                                                     <  2:5904406/150‑1 (MQ=255)
                                                                                       cgTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCt                                                                <  1:16232/149‑1 (MQ=255)
                                                                                           aCCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTTGGGGCGGGGTTTGCCGCCCAACGGCCGGCCCGGCTGGCCAAGCACCACCCCCCCGTCACCCGTGGCGCCCCAACCCGATTCGCGGGGAAAACCGCGccga                                                            >  1:5170098/1‑148 (MQ=255)
                                                                                           aCCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCAGGCCAAACTCCACCCGCCGTTCACCGTTGCCGCCACAACCCTATTCGCGGGTAAACCCGCTccaa                                                            >  2:3885699/1‑148 (MQ=255)
                                                                                           aCCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGGTTAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGGTCACCGTTGCCGCCACAACCGTATTCGCGGGTAAACCCGCAccga                                                            >  1:1448609/1‑148 (MQ=21)
                                                                                             cAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTAACCGTTGCCGCCACTACCGTATTCGCGGGTAAAACCGCTCCTaaa                                                          >  2:6273504/1‑148 (MQ=255)
                                                                                                gATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGGTCACCGTTGCCGCCACTACCGTATTCGCGGGGAAACCCGCTCCTACAggg                                                       >  2:4716818/1‑150 (MQ=255)
                                                                                                 ttGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTTTTCGCGGGTAAACCCGCTCCTACAgggg                                                      <  1:2626730/149‑1 (MQ=255)
                                                                                                  tGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAggggg                                                     <  1:5661673/150‑1 (MQ=255)
                                                                                                  tGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAggggg                                                     <  2:7286704/150‑1 (MQ=255)
                                                                                                   gcggcgGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAgggggg                                                    <  2:304235/150‑1 (MQ=255)
                                                                                                     ggcggGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTTGGAGCGGGTTTACCCCCGAAAGGGCCGGCCCAGCTGGCCAATATACAACCCCCGTTAACCCTTTCCCCCACAACCGGATTAGCGGGTAAACCCCCACCAACAAGGGGGCt                                                  >  1:5081505/1‑150 (MQ=255)
                                                                                                     ggcggGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGGTGCCGCCACTACCGGATTCGCGGGTAAACCCGCTCCTCCAGGGGGGCt                                                  >  2:6038436/1‑150 (MQ=255)
                                                                                                          ggCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTtgtgt                                             <  1:1335167/150‑1 (MQ=255)
                                                                                                          ggCGCTGGTAATCAAGGCGCCACCGGTCGCCGGTAGGAGCGGGTTTCCCCGCGACCGGGCCGGCCCTGCTGGCCAAGCCCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTtgtgt                                             <  1:1131727/150‑1 (MQ=18)
                                                                                                           gcgcTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTtgtgtg                                            <  1:4825538/150‑1 (MQ=255)
                                                                                                           gcgcTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGGTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATAACCACCCGCCGGTCACCCTTGCCCCCACTACCCTATTCGCGGGTAAACCCGCGCCAACAGGGGGGCTtgtgtg                                            >  1:2121551/1‑150 (MQ=21)
                                                                                                            cgcTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTgg                                           <  1:2510841/150‑1 (MQ=255)
                                                                                                             ggtGGTAATCAAGGCGCCACCGGTCGCGTGTTGGAGGGGGTTTCCCCGCGACAGGGCCGGCCCTGCTGGCCAATCTCCCCCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAACCCCGCTCCTACAGGGGGGCTTGTGTGGc                                          <  1:4748214/148‑1 (MQ=18)
                                                                                                             gcTGGTAATCACGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGc                                          <  2:1278315/150‑1 (MQ=255)
                                                                                                               tGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGGTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGa                                        >  2:5004655/1‑150 (MQ=255)
                                                                                                                 gTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGAtt                                      <  2:4265194/150‑1 (MQ=255)
                                                                                                                 gTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGCCCTGCTGGCCAATATCCACCCGCCGTTCACCGGTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCGTGTGTGGCGAtt                                      >  2:7588540/1‑150 (MQ=255)
                                                                                                                  tAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGTTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTg                                     <  1:8385676/150‑1 (MQ=255)
                                                                                                                       aaGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCg                                <  2:1131781/150‑1 (MQ=255)
                                                                                                                        aGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGc                               <  2:400221/150‑1 (MQ=255)
                                                                                                                         ggCGCCCCCGGTCGCCTGTAGGAGCGGGTTTCCCCGCGACAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGca                              <  1:8269238/150‑1 (MQ=255)
                                                                                                                         ggCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAACGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGca                              <  1:2949293/150‑1 (MQ=255)
                                                                                                                         ggCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGca                              <  1:668271/150‑1 (MQ=255)
                                                                                                                            gcCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACaa                           <  1:3859751/150‑1 (MQ=255)
                                                                                                                            accACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCCATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTGGTGTGGCGATTGGGCCGCACaa                           >  2:4352264/2‑150 (MQ=255)
                                                                                                                               aCCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGc                        <  2:3277613/150‑1 (MQ=255)
                                                                                                                               aCCGGTCGCCGGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGc                        <  1:3214885/150‑1 (MQ=255)
                                                                                                                                        cTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCcac               >  2:117275/1‑150 (MQ=255)
                                                                                                                                         tGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCcaca              >  1:5403462/1‑150 (MQ=255)
                                                                                                                                         tGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCcaca              >  2:5117639/1‑150 (MQ=255)
                                                                                                                                         tGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCcaca              <  2:8507430/150‑1 (MQ=255)
                                                                                                                                          gTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCCACAg             >  2:5959583/1‑150 (MQ=255)
                                                                                                                                           tcggAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCCACAgg            <  1:2274875/148‑1 (MQ=255)
                                                                                                                                           tAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCCACAgg            <  2:6803256/150‑1 (MQ=255)
                                                                                                                                                gCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCCACAGGGcgcg       >  2:8319147/1‑150 (MQ=255)
                                                                                                                                                 cGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCCCCCGCCGTTCCCCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCCACAGGGCGCgt      <  2:4324295/150‑1 (MQ=255)
                                                                                                                                                 cGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCCACAGGGCGCgt      <  1:7197739/150‑1 (MQ=255)
                                                                                                                                                    gTTTACCCGCGACCGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAACCCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCCACAGGGCGCGTGTa   <  1:73824/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCATCATGGACTACCTGAAGACCCGGGCGCCGTTCTGGAAGAAGGAGAATACTCAGGAAGGGCCAAGGTGGGTGGAGGGGCGGCAGAGTGACCAGGATGCGGCGGGGCGCTGGTAATCAAGGCGCCACCGGTCGCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACTACCGTATTCGCGGGTAAACCCGCTCCTACAGGGGGGCTTGTGTGGCGATTGGGCCGCACAAGGCCGCTCCCACAGGGCGCGTGTAG  >  NC_002947/1484941‑1485239

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: