Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,022 T→C E167E (GAA→GAG ettA ← energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,0220TC92.3% 158.5 / 3.3 65E167E (GAA→GAGettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (31/29);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

CAGGAGGAGTAGTTACCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTC  >  NC_000913/4629783‑4630252
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      
cagGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCgg                                                                                                                                                                                                                              >  1:176659/1‑250 (MQ=255)
   gaggagTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGCGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGCCGACGCTCACCACCGgag                                                                                                                                                                                                                            >  1:484363/1‑249 (MQ=255)
     ggagTAGTTGCCTTCCCTCGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTACCGGGCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGTAAGCGTTCCAGCCAGGCCGCGGATTCGGCAGCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGGTTTTCCAGCAGCAGGCGGCACAACTCTACGCGACGACGCTCACCACCGGAGAGGt                                                                                                                                                                                                                        <  1:223913/251‑1 (MQ=255)
       agtagtTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTGCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTa                                                                                                                                                                                                                      <  1:373898/251‑1 (MQ=255)
        gtagtTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAg                                                                                                                                                                                                                     <  2:126820/251‑1 (MQ=255)
                  tCCCACGGAATACCTTCCCCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTcg                                                                                                                                                                                                            <  2:271274/250‑1 (MQ=255)
                  tCCCACGGAATACATTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTcg                                                                                                                                                                                                            <  1:5988/250‑1 (MQ=255)
                    ccACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGt                                                                                                                                                                                                         >  1:391530/1‑251 (MQ=255)
                                        ggTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAgc                                                                                                                                                                                     >  1:342631/1‑251 (MQ=255)
                                          tCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAgcgc                                                                                                                                                                                   >  1:208378/1‑251 (MQ=255)
                                              gTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCCCCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGGCCGGCAGACGTAGCGCATcc                                                                                                                                                                               >  2:184762/1‑251 (MQ=255)
                                                          aGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGcacc                                                                                                                                                                   >  1:247040/1‑248 (MQ=255)
                                                                       aTCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTATACGTTc                                                                                                                                                       <  1:443941/250‑1 (MQ=255)
                                                                             gAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTc                                                                                                                                                                                                        <  1:222667/195‑1 (MQ=255)
                                                                             gAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTc                                                                                                                                                                                                        >  2:222667/1‑195 (MQ=255)
                                                                             gAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGACAGGCAGCCGTAGCGCATCCGCCGCACGCTCCAGCTGTAAGTTCAggttct                                                                                                                                                >  2:466462/1‑246 (MQ=255)
                                                                               aggaaCGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTGGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGa                                                                                                                                              <  2:482962/248‑1 (MQ=255)
                                                                               aGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGa                                                                                                                                              <  1:448331/251‑1 (MQ=255)
                                                                                    cGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGCGAGGTTAGCGCTTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCCGGTTCTGACcgtc                                                                                                                                         >  2:29368/1‑251 (MQ=255)
                                                                                     ggTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGCAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCCTCCGCCGCACGCTCCAGCTGTACGTTCAGGTAAGGACcgtcg                                                                                                                                        >  1:6931/1‑251 (MQ=255)
                                                                                     ggTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACcgtcg                                                                                                                                        >  1:44867/1‑251 (MQ=255)
                                                                                         gTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGCTTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCGGGCGGCACAACGCTACGCGCCGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGACCGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTg                                                                                                                                      >  1:146152/1‑249 (MQ=255)
                                                                                                aTCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAAt                                                                                                                             <  1:352222/251‑1 (MQ=255)
                                                                                                   gCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCt                                                                                                                                                               >  1:85050/1‑214 (MQ=255)
                                                                                                         ccGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCTGCCAGGTGGTTGGTCGGTTCGCCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTc                                                                                                                    <  2:273595/250‑1 (MQ=255)
                                                                                                           ggTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCa                                                                                                                  <  1:184762/251‑1 (MQ=255)
                                                                                                                    gAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCtt                                                                                                         >  1:356292/1‑251 (MQ=255)
                                                                                                                      aGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCCTGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGt                                                                                                       >  1:453283/1‑251 (MQ=255)
                                                                                                                         cGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTGCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTc                                                                                                     <  2:250095/250‑1 (MQ=255)
                                                                                                                          gTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCCCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGCATGATCTCTTCCAGCCGGCCTTTTTccg                                                                                                   >  2:235584/1‑249 (MQ=255)
                                                                                                                             cAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCTAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGCCGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCgg                                                                                                >  1:206908/1‑251 (MQ=255)
                                                                                                                              aGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGc                                                                                               >  1:39596/1‑251 (MQ=255)
                                                                                                                              aGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCGGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATc                                                                                                                                                                                <  1:177837/170‑1 (MQ=255)
                                                                                                                              aGGAAGCGTTCCAGCCAGGCCACGAATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACACCGCTACGCGACGACGCTCACCACCGGAGAGGGTAGCGATTTTCGCGTCCCAGTACGGCAGACGTAGCGCATc                                                                                                                                                                                >  2:177837/1‑170 (MQ=255)
                                                                                                                               ggAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATc                                                                                                                                                                                <  1:170161/169‑1 (MQ=255)
                                                                                                                                 aaGCGTTCCAGCCAGGCCACGGTTTCGGCAGCCAGGTGGTTGTGCGGTTCGGCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAg                                                                                            <  2:195615/251‑1 (MQ=255)
                                                                                                                                     gTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGCGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTTATCAGCGGccgacttt                                                                                        >  2:386209/1‑245 (MQ=255)
                                                                                                                                                 ccACGGGTTCGGCAGCCAGGTGGGTGGTCGGTTCGTCGAGCAGAAGCATGTCTGGTTTTTCCAGAAGCAGGCGGCCAAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCa                                                                            <  2:453283/251‑1 (MQ=255)
                                                                                                                                                  cACGGTTTCGGCATCCAGGGGTTTGGTCGTTGCGTCGAGCAGAAGCATTCCTGGTTTTTCCAGCAGCAGGCGGCCCAACGCTTCGCGACGACGCCCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCGTGTTCAGCGGCCAGCTTGTCAAAATCGGCAt                                                                           <  2:294896/251‑1 (MQ=255)
                                                                                                                                                   ccggtgtCGGCATCCAGGTGGTTGTGCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTTCGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATc                                                                          <  2:391252/245‑1 (MQ=255)
                                                                                                                                                     ggATTCGGCACCCAGGGGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGGCGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGAGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCg                                                                        <  2:191428/251‑1 (MQ=255)
                                                                                                                                                      gATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCgg                                                                       >  1:264525/1‑251 (MQ=255)
                                                                                                                                                          cGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGCTTTTCGCGTCCCAGTCCGGCAGAAGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGCGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCg                                                                   >  1:290784/1‑251 (MQ=255)
                                                                                                                                                          cGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCg                                                                   >  1:149701/1‑251 (MQ=255)
                                                                                                                                                                ccatgtggtTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTa                                                              <  1:40601/246‑1 (MQ=255)
                                                                                                                                                                 cAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCCAAATCGGCATCCGGATCGGCGTAc                                                             >  1:364819/1‑250 (MQ=255)
                                                                                                                                                                   ggtggtTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGAC‑‑AGCGCATCCGCCGCACGCTCCCGCTCTACGATCAGATTATGACCGTCGTGTGCCTGCATGCTCTCTTCCAGCCGCCCTTGTTCAGCGGCCAGCTTGTCACAATCGTCATCCGGATCGGCGttcctcccc                                                       >  2:114207/1‑243 (MQ=255)
                                                                                                                                                                    gtggtTGGTCGGTTCGTAGAGCAGCCGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGGCGTTCACCACCGGAGAGGTTAGCGCTTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCCGACTGATCTCTTCCAGACGGCCTTGTTCAGAGGCCAGCTTGTCAACATAGGCATCCGGATCGGCGtcctgcg                                                         >  2:38275/1‑245 (MQ=255)
                                                                                                                                                                      ggttggtCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAGTGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCTTACAgcgc                                                        >  1:67634/1‑250 (MQ=255)
                                                                                                                                                                         tggtCGGTTCGTCGGGCAGCAGCATGCCTGGTTGTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGCTAGCGATTTTCGCGTCCCAGTCCGGCAGACGGAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATa                                                     <  2:54431/250‑1 (MQ=255)
                                                                                                                                                                         tggtCGGTTCGTCGAGCAGCAGACGGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGAGACGACGCTCACCACCGGAGAGGTTAGCTATTTTCGCGTCCCAGTCCGGCAGACGGAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTGCAGCGGCCAGCTTGGCAAAATCGGCATCCGGATCGGCGTACAGCGCATa                                                     <  2:398856/250‑1 (MQ=255)
                                                                                                                                                                                ttCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCCTCCGGATCGGCGTACCGCGCATACACTTCAt                                             >  2:260506/1‑251 (MQ=255)
                                                                                                                                                                                 tcgcggaggagcggcaTGGCTGGTTTTCCCAGCAGGCGGCGGCACAACGCTTCGCGACGACGCTCCCCACCGGAGAGGTTAGCGAGTTGCGCGTCCCAGCCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCGTGAATGATCTCTTCCAGACGGCCGTGTTCAGCGGCCAGCTTGTCAAAAGCGGCAGCCGGATCGGCGTACAGCGCATactcgtcat                                             <  2:233067/237‑8 (MQ=255)
                                                                                                                                                                                   gtcgAGCAGCAGCATGTCTGGGTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCGGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGTATCGGCGTACAGCGCATACACTTCATCCa                                          <  1:386209/251‑1 (MQ=255)
                                                                                                                                                                                                  gTCTGTTTTTTCCAGGAGCGGGCGGACAACCTCTACGCGACGCCGCCCACCCCCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTAAGACCGGCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCTGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCg                            <  2:41996/250‑1 (MQ=255)
                                                                                                                                                                                                  gTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGt                           >  1:409157/1‑251 (MQ=255)
                                                                                                                                                                                                      ggTTTTTCCTGCAGCGGGCGGCACAACGCGACGCGCCGACTCTCCCCACCGGAGAGTTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGAAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTaa                        <  2:296444/250‑1 (MQ=255)
                                                                                                                                                                                                      ggTTTTCCCAGCAGCAGGCGGCACAACGCTACTCGACGACGCTCACAACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCTCATCCGCCGCACTCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAAc                       <  1:381376/251‑1 (MQ=255)
                                                                                                                                                                                                             cccgcagatggcggcACAACGCTACGCGTCGACGCTCACCACCGGAGAGGTTAGCGATTGTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTca                <  2:478320/242‑2 (MQ=255)
                                                                                                                                                                                                              cagcagcagGCGGCACAACTCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTcag               <  1:416181/251‑3 (MQ=255)
                                                                                                                                                                                                                   gcagGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAAcc          <  1:112306/251‑1 (MQ=255)
                                                                                                                                                                                                                   gcagGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCCGGCGTTTCAGGGAGTTAACCACTTCAGAAAcc          >  1:138832/1‑251 (MQ=255)
                                                                                                                                                                                                                    cagGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGACCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTAAGGGCGTTACCCACTTcagcaac           >  2:21304/1‑243 (MQ=255)
                                                                                                                                                                                                                         ggcACAACGCTACGCGACGACGCTCACCACCGGCGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGcttct    >  1:194252/1‑251 (MQ=255)
                                                                                                                                                                                                                           cacaACGCTACGCGACCACGCTCCCCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCCCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCGTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCttcttc  <  2:223627/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      
CAGGAGGAGTAGTTACCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTC  >  NC_000913/4629783‑4630252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: