Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,304 A→G Y73Y (TAT→TAC)  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,3040AG98.6% 185.4 / ‑6.0 69Y73Y (TAT→TAC) ettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (42/26);  minor base T (0/1);  total (42/27)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

TCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTA  >  NC_000913/4630059‑4630545
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
cccGGCAGACGCAGCGCAGGCGCCGCACGCCCCAGCTGTACTTATAGGTTAGGACCGTCGTGAGCCTGAAGGATCTCTTCCAGACTGCCTTGTGCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACGTCAGAATCCGCTTCTGCAATTGACTCACGCACTGTGTGCGCCGGGTTCAGCTGCGGTTCCTGATGCAGGTAAcc                                                                                                                                                                                                                                                 <  2:93264/250‑1 (MQ=255)
    gCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGCCCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAAcccatt                                                                                                                                                                                                                                             >  1:598401/1‑247 (MQ=255)
      aGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCtt                                                                                                                                                                                                                                           >  2:493475/1‑251 (MQ=255)
      aGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGAGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAACCCATCtt                                                                                                                                                                                                                                           >  2:178694/1‑251 (MQ=255)
        aCGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACCCTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGa                                                                                                                                                                                                                                         <  1:539317/251‑1 (MQ=255)
            aGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTc                                                                                                                                                                                                                                     >  1:363921/1‑251 (MQ=255)
               gcATCCGCCGCACGCGCCAGCTGTACGTTCAGGTTAGGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCGTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAATCCGCTTCTTCAATGGACTCACGCACGGGGTTTTGCGGGTTCAGCTGCGGTGCCTGCGGCAGGTAACCAATCTTGATGTCTgg                                                                                                                                                                                                                                  <  2:65667/251‑1 (MQ=255)
                 aTCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGTCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCACTGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGCTGTCTGG‑cc                                                                                                                                                                                                                               >  1:146700/1‑250 (MQ=255)
                         cACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGAcgc                                                                                                                                                                                                                       <  1:546776/251‑1 (MQ=255)
                          aCGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCACAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGAcgcg                                                                                                                                                                                                                      >  2:585617/1‑251 (MQ=255)
                          aCGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGCTCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGACTGG‑CTGCGGAcgcg                                                                                                                                                                                                                      >  1:623969/1‑251 (MQ=255)
                          aCGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCACTCTTGGAGTCTGG‑CTGCgggcggg                                                                                                                                                                                                                      >  2:131637/1‑246 (MQ=255)
                            gCTCCAGCTGTACTTTGAGGTTAGGACCGTCGTGAGCCTGAATGAGCTCTTCCAGACGGCCGTGGTCAGCGGCCAGCTTGGCAAAATCGGCATCCTGATCGGCGTACAGCGCATACACTTCCTCCGGGCTTTTCAGGGCGTTAACCTCTTCAGAAGCCGCTGCGTCAAGGGACTCACGCACGGTGTGTTCCGGGTTCAGCGGCGGCTCCGGCGTCAGGTAACCAATCTTGATGTCTTA‑CTGCGGTcgctct                                                                                                                                                                                                                    <  2:555692/251‑4 (MQ=255)
                              tCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCCGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGAGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATATTGATGTCTTG‑CTGCGGACGCGCTT‑c                                                                                                                                                                                                                 >  2:606053/1‑251 (MQ=255)
                                 aGCTGTACGTGCAGGTTATGACCGTCGTGAGCCGGAATGATCTCTTCCAGACGGCATTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGGACAGCGAATACACTTCATCCAGGCGTGTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCTGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACTCTCTT‑CAcc                                                                                                                                                                                                              <  2:356828/251‑1 (MQ=255)
                                   cTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCCATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCtt                                                                                                                                                                                                            >  2:124364/1‑251 (MQ=255)
                                     gTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCTCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCg                                                                                                                                                                                                          <  1:27854/251‑1 (MQ=255)
                                     gTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCg                                                                                                                                                                                                          >  1:232457/1‑251 (MQ=255)
                                             aggttATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑Cttt                                                                                                                                                                                                 <  1:286501/248‑1 (MQ=255)
                                                ttATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTTATGTCTGG‑CTGCGGTCGCTCTT‑CCCCTTCGATGT‑CTTTATc                                                                                                                                                                                              >  1:404461/1‑251 (MQ=255)
                                                   tGACCGTAGTGAGCCTGAATGATCTCTTCCAGCCGGACTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCCT‑CACCTTCGCTGT‑CTTTATccat                                                                                                                                                                                           >  2:336537/1‑248 (MQ=255)
                                                    gACCGTCGTGAGCCTGAATGAACTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCGT‑CACCTTCTAGGT‑CTTTATCAAt                                                                                                                                                                                           <  2:72931/250‑1 (MQ=255)
                                                           gtGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑ATTTTTCAATGCCCGcc                                                                                                                                                                                    >  2:188495/1‑250 (MQ=255)
                                                           gtGAACCTGAATGATCTCTTCCAGACGGCCTTGTTCACAGGCCAGCTTGTCACAATCGGCATCCGGATCGGCGGACAGCGCATAAACTTCATCCCGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCCGCTGCGGTTCCTGCGGCAGGGAACCCATCTTTATGCCTGG‑CTACGGAGGCGCTTTCAC‑TTCTATGG‑CTTTATCAAtcgccgccc                                                                                                                                                                                   >  2:425866/1‑243 (MQ=255)
                                                            tGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCAt                                                                                                                                                                                  >  2:473122/1‑251 (MQ=255)
                                                                   gAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATgcgc                                                                                                                                                                           >  1:188267/1‑251 (MQ=255)
                                                                   gAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATgcgc                                                                                                                                                                           >  1:524479/1‑251 (MQ=255)
                                                                    aaTGATCCCTCCCAGACGGCCTTGTTCAGCGGCCAGCTTGTAAAAATCGCCATCCGGATCGGCGTCCAGCGCATACCCTCCATCCAGGCGTTTCAGGGCGTTAACCCCTTCAAAAACCGCTTCTTCAAGGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CGGCGTACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCgca                                                                                                                                                                          <  1:606053/251‑1 (MQ=255)
                                                                        aTCTCTTCCAGACGGCCTTGTTCAGCGGCCCGCTTGACAAAATCGTCATCCGAATAGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCCGAAACCGCTTCTTCAATGGACTCACGCACGGAGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAAGCTTGATGTCTGG‑CTGCGGACGCCCTT‑CACCTTCGCTGT‑ATTTCTCAATTCCCGCCATAATgcacctc                                                                                                                                                                        >  2:55064/1‑244 (MQ=255)
                                                                           tctTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGt                                                                                                                                                                   >  2:178578/1‑251 (MQ=255)
                                                                                          tGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGcc                                                                                                                                                    >  1:321023/1‑251 (MQ=255)
                                                                                           gTTCGGCGGCCAGCTTGACAAAATCGGCCTCCGGAGCGCCTTACAGCCCATCCCCTTCACCCCGGCTTTTCAGGGCTTTAACCTCTCCAAAACCCGCTTCTTCAATGGACTCCCGCACGTTGTTTTCCGGGTTCATCTGCGGTTCCTGCGGCGGGTACCCATTCTTGATGTCTGG‑CTGCGGACGCGTTT‑CACCTCCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGcc                                                                                                                                                    <  1:24686/250‑1 (MQ=255)
                                                                                               cgcggcaagtttgaaaaaaTCGGCATCCGGATCGGCGTACAGCGCATCCACTTCCTCCAGGCGTTCCGGGGCGTTAACCACTTCAGAAACCTCTTCTTCAAGGGACTCACGCGCGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTGGATGTCGGG‑CTGCGGACGCGCTC‑CACCTTCGATGT‑CTTTAACAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCa                                                                                                                                               <  2:188267/236‑1 (MQ=255)
                                                                                                 cgtccAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCCGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCACTGCCCGCCATAATGCTCAGCAGGGTAGACTTACCCGCGCCGTTccg                                                                                                                                              >  2:161929/4‑248 (MQ=255)
                                                                                                          tgccaaaaccgccaccgggaccGGCGTACAGCGCATACACTCCACCCAGGCGTTTCAGGGCGTAAACCACTTCAGAAACCGCTTCTTCAATGGCCTCCCGCACGGTGTTTTCCGGGTTCAGCTGCGTTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGCCGCGCTT‑CACCTGCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCTCCGTTCAGACCCAGGa                                                                                                                                      <  1:505387/228‑1 (MQ=255)
                                                                                                          tGTCACAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGACCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCCGACCCAGGac                                                                                                                                     >  2:625231/1‑250 (MQ=255)
                                                                                                           gTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCa                                                                                                                                                                                             <  1:91476/193‑1 (MQ=255)
                                                                                                           gTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCa                                                                                                                                                                                             >  2:91476/1‑193 (MQ=255)
                                                                                                                   cGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAAttttt                                                                                                                           >  1:559731/1‑251 (MQ=255)
                                                                                                                        cccGGATCGGCGTCCAGCGCCTACCTTTCCCCCAGGCGTTTCAGGGCGTTAACCACTTCAGAACCCTCTTCTTCAATGGACCCCCCCCCGGTGTGTTCCGGTTCCAGCTGCGGTTCCTCCGCCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑TTTTATCACTGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCTCGCCGTTCAGACCCAGGACACCAATTTTTGcccc                                                                                                                      <  1:178694/250‑1 (MQ=255)
                                                                                                                             aTCGGCGTACAGCGCATACACTTCATCCATGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCATGCGGCAGGTAACCAATATTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCCATGACCGCCATAAAGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGAACCAGGGCACCAATTTTTgagatcggc                                                                                                                  >  2:102333/1‑242 (MQ=255)
                                                                                                                                  cGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGgaagaa                                                                                                             >  2:476388/1‑250 (MQ=255)
                                                                                                                                  cGTACAGCGCATACACTTCATACAGACGGTTCAGGGCGTTAACCAATTCAATAACCGCTTCTTCATTGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTTTGG‑ATGCGGACGCGCTT‑CACCTTCGCAGT‑CTTTATCAATGACCGCCATAATGCGCCGCAGGGGAGACTTACCCGCGCCGTTCATACACAGGACAACACTTTTTGTCCCAGGgcagaa                                                                                                             >  2:42303/1‑250 (MQ=255)
                                                                                                                                    tACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGCTGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGCCACCAATTTTTGCCCCAGGGAAGCAACt                                                                                                          >  2:164643/1‑251 (MQ=255)
                                                                                                                                            aTACACTTCATCCAGGCGTTTCAGGGCGTTCACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCagag                                                                                                     >  2:86699/1‑248 (MQ=255)
                                                                                                                                                    cATCCAGGCTTTTCAGGGCGTTACCCACTTCAAAAACCGCTCCTTAAATGGACTCCCGCACGGTGTGTTCCGGGTCCAGCGGCGGTCCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGCCGCGCTT‑CCCCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTc                                                                                          <  1:457412/251‑1 (MQ=255)
                                                                                                                                                        cAGGCGTATCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGAATTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCaaa                                                                                       >  2:572452/1‑250 (MQ=255)
                                                                                                                                                           gCGTTTCAGGGCGTTAACCCCTTCAGAAACCGCTTCTTAAATGGACTCACGCGCGGTGTGTGCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATAAAGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAatat                                                                                   <  1:473122/251‑1 (MQ=255)
                                                                                                                                                                  aGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGttt                                                                            >  2:113107/1‑251 (MQ=255)
                                                                                                                                                                   gggCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGtttt                                                                           >  2:213968/1‑250 (MQ=255)
                                                                                                                                                                     gCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCCATCTTGCTTTTTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATACCCCCCATAATGCGCAGCAGGGGAGACTTGCCCGCGCCGTTCAGACCCCATACACCAATTTTTGCCCCAGGGAAGAAAATCAGAGAGATGTTTTTCAAAATATGGAGTTTcgg                                                                         >  2:336752/1‑251 (MQ=255)
                                                                                                                                                                       gTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑ATGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCACGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTcggcg                                                                       >  2:96131/1‑251 (MQ=255)
                                                                                                                                                                            cccctgcagAAACCGTTTTTTCAATGTACTCCCGCACGGTGGGTTCCGGGTTCAGCGCCGGTCCCGGCGGCAGGTAACCAATTTTGATGTCTGG‑CGGCGGACGCGCTC‑CCCCTTCGATGT‑CTTTACCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGaaca                                                                  <  2:488999/243‑1 (MQ=255)
                                                                                                                                                                              aCTTCAGAAACCGTTTCTTCAAGGGCCCCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CCCCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGaacaac                                                                <  1:493475/251‑1 (MQ=255)
                                                                                                                                                                                 tcagAAACCGCTTCTTCAATGGACTCACGCACGGCGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTCACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACCCCAATTTTTGCCCCCGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACttt                                                             >  1:231638/4‑251 (MQ=255)
                                                                                                                                                                                   agAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAAAAACTTTGc                                                           >  1:399653/2‑251 (MQ=255)
                                                                                                                                                                                        ccgttttttcAATGGCCTCACGCTCGTTGTGTTCGGGGTTCGCCGGCGGTTCCGGCGGCAGTTAACCACT‑TGGGTTGCTGGGCGGCGGACGCGCTT‑CACCTGCGTTGT‑CTTTACAAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGTCACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATTACGTTTCGGCGGAACAACTTTGCCGaca                                                      <  2:429461/244‑1 (MQ=255)
                                                                                                                                                                                         cGCTTCTTATATGGACTCACTCTGGGTGTGTTCCGGGTTCAGCTGCGGTTCCCGCGGCAGGTAACCATGCTTTATGTCTGG‑CTGCGCACGCGCTC‑CCCCCTCG‑TGTTCTTTCTAAATGTCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGAGGTTTTGCAAAAGATGACGTTTCGGCGGAACAACTTTGCCGacac                                                     <  2:248148/251‑1 (MQ=255)
                                                                                                                                                                                             tcttcaggggccccccgcacgGTGTGTTCGGGTTTCAGCTGCGTGTCCTGCGGCAGGTAACCAATCTTGATGTCGGG‑CTGCGGACGCGCTC‑CACCTTCGATTT‑CTTTACCAGGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGAt                                                  <  1:86699/235‑1 (MQ=255)
                                                                                                                                                                                                tcAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCAt                                              >  1:23611/1‑251 (MQ=255)
                                                                                                                                                                                                                         cggggTTCAGGTGTGGTCCCTGCGGCAGGTAACAAATCTTGAGGTCGGG‑CGGGGGACGCGCTT‑CTCCTTCGATGT‑TTTTAAAAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTg                     <  1:178578/249‑1 (MQ=255)
                                                                                                                                                                                                                          cGGGTTCAGCTTCGGTTCCTGCGGCAGGTAACCAATCTTGGTGTCTGG‑CTGCGGACTCGCTT‑CCCCTTCGATGT‑CTTTCTCAATTCCCGCCATAATGCGCAGCAGGGTAGACTTACCCCCGCCGTTCAGACCCAGGCCACCAATTTTTGCCCCAGGGAAGAAACTCCGAGAGATGTTTTTCAAAATATACCGTTTCGGCCGAACACCTTTGCCGACACGATGCATGGTATAAACGAATTTAGCCACGtttg                    >  2:42/1‑249 (MQ=255)
                                                                                                                                                                                                                           gggTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACg                                    >  1:548369/1‑234 (MQ=255)
                                                                                                                                                                                                                              ttCAGTTGGGGTTCCTGCGGCGTGTAACCAATCTTTATGTCTGG‑CTGCGGACGCGCTT‑CACCTGCGGTGT‑CTTTCTCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACt                 <  1:249573/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gCTGCGGTTCCGGCGGAAGGTAACCAATCTGGATGTCTGG‑CGGCGGCCGCGCTC‑CCCCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGcc            <  1:353146/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gCGGTTCCGGCGGCAGGTACCCAATCTTGTTGTCGGG‑CGGCGGCCGCGCTC‑CCCCTGCGATGT‑TTTTACCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCct           <  1:585617/249‑1 (MQ=255)
                                                                                                                                                                                                                                      cGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CCCCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCCTGGTATAAACGTATTGAGCCACGTTGGACTTCGCctctt        >  1:822/1‑250 (MQ=255)
                                                                                                                                                                                                                                          tCCTGCGGCAGGTAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGtt    >  1:321708/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gCGGCAGGTACCCAATCTTGATGTCTGG‑CTGCGGCCGCGCTT‑CCCCTGCTGTGT‑CTTAATAAAGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTa  <  1:444012/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
TCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGG‑CTGCGGACGCGCTT‑CACCTTCGATGT‑CTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTA  >  NC_000913/4630059‑4630545

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