Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,238 T→A S95S (TCA→TCT 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,2380TA97.4% 233.2 / ‑2.0 76S95S (TCA→TCTettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (31/43);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGAC  >  NC_000913/4629991‑4630467
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      
gcagcggGCGGTACAACGCTACGCGACGACGCTCACCCCCGGGACGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGGTTCAGGGCGTTAACCACTTCAGaa                                                                                                                                                                                                                                     <  2:497194/251‑1 (MQ=255)
     agGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGc                                                                                                                                                                                                                                >  1:605671/1‑251 (MQ=255)
      ggcggcACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGct                                                                                                                                                                                                                               <  1:80757/251‑1 (MQ=255)
       gcggcACAACGCTACGCGACGACGCTCCCCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACTTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGctt                                                                                                                                                                                                                              <  1:494548/251‑1 (MQ=255)
        cggcACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCCGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGAACGTCGTGAGCCTGCATGCTCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGCACAGCGCATACACTTCATCCAGGCGTTTCTGGGCGTTAACCACTTCAGAAACCGcttc                                                                                                                                                                                                                             >  2:618337/1‑251 (MQ=255)
         ggcACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGcttct                                                                                                                                                                                                                            <  1:324613/251‑1 (MQ=255)
         ggcACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGcttct                                                                                                                                                                                                                            >  1:675243/1‑251 (MQ=255)
            acaACGCTACGCGACGACGCGCACCCCCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCa                                                                                                                                                                                                                         <  1:70242/251‑1 (MQ=255)
            aaaaagCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCa                                                                                                                                                                                                                         <  1:654461/246‑1 (MQ=255)
             caACGCTACGCGACGACGCTCACCACCGGAGAGGGTAGCGCTTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGATATGACCGTCGTGCGCCTGAATGATCTCTTCCAAGCTGCCTTGTTCAACGGCCAGCTTGTCAAAATTGGCATCCGGATCGGCGGACAGCGCATCCACTGCTTCCAGGCGTTTCAGGGCTTTAACCCCTTCAGAACCCTCTTCTTCac                                                                                                                                                                                                                        >  2:328420/1‑250 (MQ=255)
                 gCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATgg                                                                                                                                                                                                                     >  1:74999/1‑250 (MQ=255)
                   caCGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCTTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACt                                                                                                                                                                                                                  <  2:202063/250‑1 (MQ=255)
                    ccgcgGCGACGCTCACCACCGGGTAGCTTAGCGGTTTGCGCTGCCCAGTCCGGCAGACGGTGCGAATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCGTAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTACCCACTTCAGAAACCGCTTCTTCAATGGAc                                                                                                                                                                                                                   <  2:416434/248‑1 (MQ=255)
                     cgcgACGACGCGCACCACCGGAGAGGATAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCCGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTACAGACGGCCTTGTTCAGCGGCCAGCTTGTCAACATCGGCATCCAGATCGGCGAACAGCGCATACACTTCCTCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTc                                                                                                                                                                                                                 >  2:60363/1‑250 (MQ=255)
                            acgctCACCACCGGAGAGGTCAGCGATTCTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGAACgg                                                                                                                                                                                                         <  1:248475/247‑1 (MQ=255)
                                ccaccacCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTGCAGAAACCGCGTCTGCAATGGACTCACGCACGgtgtg                                                                                                                                                                                                     <  1:529548/250‑1 (MQ=255)
                                 gaccccaggAGGGGTTAGCGTGTTTCGCGTCCCAGTCCGGCAGACGTGGCGATCCCGCCGCCCACTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCGCTGCCAGACGGCCGTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTGCAATGGACTCACGCACGgtgt                                                                                                                                                                                                      <  2:204135/242‑1 (MQ=255)
                                    caccGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTcc                                                                                                                                                                                                 <  1:190723/251‑1 (MQ=255)
                                     accGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAACTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCCTCCAGGCGTTTCAGGGCGGTAACCACTTCAGAAACCGCTTCTTCAATGGACTAAAGCACGGTGTGTTc                                                                                                                                                                                                  >  2:656151/1‑249 (MQ=255)
                                      ccGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGAACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACCCTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTACgg                                                                                                                                                                                               >  1:68022/1‑251 (MQ=255)
                                      ccGGAGAGGGTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCgg                                                                                                                                                                                               >  1:144671/1‑251 (MQ=255)
                                       cGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCggg                                                                                                                                                                                              <  1:28041/251‑1 (MQ=255)
                                        ggAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGt                                                                                                                                                                                             >  1:360108/1‑251 (MQ=255)
                                          agagGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGAATACCCTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTACCGCACggggacgtccggtgtc                                                                                                                                                                                           >  2:257040/1‑237 (MQ=255)
                                            gatgttatcggttttCGCGTCCCCGTCCGGCAGACGTGGCGCAGCCGCCGCACGCGCCATCTGTGCGTTCTGGTTAGGCCCGTCGTGAGCCTGAAGTAGCTCTTCCATACGGCCGTGTTCAGTGTCGAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATGCACTTCATCCAGGCGTTTCAGGGCGTGATCCTCTTCAGAAACCGCTTCTTCAATGGACTTACGCACGGTGTGTTCCTGGTTCa                                                                                                                                                                                          <  2:606293/240‑1 (MQ=255)
                                               gTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCt                                                                                                                                                                                       >  1:584174/1‑251 (MQ=255)
                                                ttAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGATATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTCACCACTTCAGCAACCGCTTCTTCAATGGACTCTCACACGGTGTGTTCCGGGTTCAACTg                                                                                                                                                                                      >  1:391097/1‑251 (MQ=255)
                                                  aGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCg                                                                                                                                                                                    <  1:440871/251‑1 (MQ=255)
                                                      atttttccgTCCCCGTCCGGCAGACTTAGCGCATCCGCCGCACGCTCCAGCTGTACGGTTAGTTTATGACCGTAGTGAGCCTGAAGTATCGCTTCCATACGGCTTTGTTCAGCGGCCAGCTTGTCAAAATCGGCAGCCGGATCGGCGTACAGTGCATTCACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTCCTTCAATGGACTCACGCTCGGTGTGTTCCGGTTTCAGCTGGGGTTc                                                                                                                                                                                <  1:517269/244‑1 (MQ=255)
                                                          tCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTg                                                                                                                                                                             <  1:640754/250‑1 (MQ=255)
                                                            gcgTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCgg                                                                                                                                                                          <  1:68692/251‑1 (MQ=255)
                                                                cccAGGCCTGCAGACGTAGCGCATCCGCCGCCCGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCGCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGCTCGGCGTACAGCGCATACACTTCATCCAGGCTTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAgg                                                                                                                                                                     <  2:311188/251‑2 (MQ=255)
                                                                 ccAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAggg                                                                                                                                                                    >  1:266821/1‑249 (MQ=255)
                                                                 ccAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACAGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAACCCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAggg                                                                                                                                                                    >  1:401763/1‑249 (MQ=255)
                                                                    gTCCGGCAGACGTGGCGGATCCGCCGCTCGCCCCAGCTGTACGTTCAGGTTATGACCGTCGTGGGCCTGAATTGTCCTTTCCAGACGGCCTTGTTCAGCTGCCAGCTTGGCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCTGGTTCAGCTGCGGTTCCTGCTGCAggttac                                                                                                                                                                 <  2:459420/251‑6 (MQ=255)
                                                                    gTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAAc                                                                                                                                                                 <  1:173131/251‑1 (MQ=255)
                                                                       cgtcAGACGTTGCGCTTCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGAACTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTGCAGCTGCGGTTCCTGCGGCAGGTAACCaa                                                                                                                                                              <  2:675243/248‑1 (MQ=255)
                                                                         gCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACAACTTCAGAAACCGCTTCTTCACTGGACTAACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCACGAAACCAATc                                                                                                                                                            >  2:551516/1‑251 (MQ=255)
                                                                          cAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCt                                                                                                                                                           <  1:12626/251‑1 (MQ=255)
                                                                              cGTAGCGCATCCGCCGCACGCTCCAGCTGTCCGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACAGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAgggaaccacttctttt                                                                                                                                                       >  2:542660/1‑236 (MQ=255)
                                                                                 aGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTc                                                                                                                                                    <  1:504556/251‑1 (MQ=255)
                                                                                   cgcATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTg                                                                                                                                                  <  1:58905/251‑1 (MQ=255)
                                                                                     cATCCGCCGCACGCTCCAGCTGTACGTTCATGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCTCTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTGCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGc                                                                                                                                                <  2:563724/251‑1 (MQ=255)
                                                                                             cgcACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGAc                                                                                                                                         >  2:113674/1‑251 (MQ=255)
                                                                                                  catcctgCTGTACCTTCCGGTTATGACCGTCGTGCGCCTGAATGATCTCTTCCAGACGGGCTTGTTCAGCGGACATATTTTCAAAAACGTCTACCGGATCGGCGTACAGCGCATACTCTTCATCCAGGCGTTTCAGGGCGTTAACCACCTCAGAAACCGCTTCTTCAATGGACGCACGCTCGGTGTGTTCCGGGTTCAGCTGCGGTCCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTTCGGTCGCGCt                                                                                                                                    <  2:245645/245‑1 (MQ=255)
                                                                                                        gCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCt                                                                                                                              <  1:51222/251‑1 (MQ=255)
                                                                                                           gTACGTTCAGGTTATGACCGTCGTGAGCCTGAATTATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGTCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCtt                                                                                                                             <  1:430502/249‑1 (MQ=255)
                                                                                                                ttCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTc                                                                                                                      >  1:605704/1‑251 (MQ=255)
                                                                                                                  caggttATGACCGTCGTGATCTTTAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACCCGTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCCCGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCtt                                                                                                                    <  2:124439/247‑1 (MQ=255)
                                                                                                                  cagggtttgacctcccggAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCTGCCATATTGTCAAAATCGGCATCCGGCTCGGCGTACAGCGCATTCACTTAATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGGTTCTTCAATGGACTTACGCTCGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCtt                                                                                                                    <  2:74999/234‑1 (MQ=255)
                                                                                                                        aTGATCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCa                                                                                                               >  1:36638/1‑250 (MQ=255)
                                                                                                                        aTGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCa                                                                                                               <  1:380473/250‑1 (MQ=255)
                                                                                                                         tGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGATCCCATCTTGATGTCTTGCTGCGAGCGCGCtgtccctttggattctttcacca                                                                                                              >  2:626187/1‑228 (MQ=255)
                                                                                                                            ccgccgtGTGCCTTAATGATCTCTCCCAGACGGCTTGGTTCAGCGGCCAGCTTGTCAAAATCTGCATCCGGATCGGCGTACAGCGCATACCCTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACTCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGcc                                                                                                          <  2:331740/247‑1 (MQ=255)
                                                                                                                             cggcctgatcccgAATGATCGCTTCCAGACGGCCTTGTTCAGCGGCCATCTTGTCAAAATCTGCATCCGTATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACGTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCTGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGccc                                                                                                         <  2:244712/239‑1 (MQ=255)
                                                                                                                                 gtGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATgg                                                                                                                                                                                                                     >  1:627347/1‑140 (MQ=255)
                                                                                                                                 gtGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATgg                                                                                                                                                                                                                     <  2:627347/140‑1 (MQ=255)
                                                                                                                                            tGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCgcagc                                                                                          >  1:280636/1‑251 (MQ=255)
                                                                                                                                                       aGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCt                                                                                                                                                   >  1:556361/1‑183 (MQ=255)
                                                                                                                                                       aGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAACCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCt                                                                                                                                                   <  2:556361/183‑1 (MQ=255)
                                                                                                                                                                          ccaccTTGTCAAACCCGGCACCCGGACCGGCGTACGGCGCATCCCCTCCACCCGGGCGTTCCAGGGGTTTACCCACTCCAGAACCCGCTTCTTCAAGGGCCTCACGCAGGGTGTGTTCCGGGTTCAGCTGCGTTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGCCGCGCTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGAccc                                                            <  1:531106/247‑1 (MQ=255)
                                                                                                                                                                                 gTCAAAATCGGCCTCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCATTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAAACAATCTTGATGTCTGGCTGCGGCCGCGCCTCACCTTCGATGTCTTTCTCAATCCCCCCCATAATGCGACGCAGGGTAGAATTGCCCGCGCCGTTCAGAcaccgggccc                                                     >  2:202344/1‑242 (MQ=255)
                                                                                                                                                                                                 ggtgcGGCGTACAGCGCATACCCTTCACCCAGGCTTTTCAGGGCTTTAACCACTTCAAAACCCGTTTCTTCAATGGACTCCCGCACGGTGTGTTCCGGGTTCAGCGGCGTTTCCTGCGGCAGTTAACCAATCTTTATGTCTGGCGGCGGCCCCGTTTCCCCTTCGATTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAgg                                     <  1:202344/247‑1 (MQ=255)
                                                                                                                                                                                                    ccGGCGTACAGCGCAACCCCTTCATCCAGGCGTTCCAGGGCGTTACCCCCTCCAAAACCCTCTTCTTCAATGGACCCCCGCAGGGTGTGTTCCGGGTTCAGCTGCGTTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAAGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGgaa                                  <  1:618337/250‑1 (MQ=255)
                                                                                                                                                                                                      ggCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGt                                                                                                                                                                                             >  2:431174/1‑95 (MQ=255)
                                                                                                                                                                                                      ggCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGt                                                                                                                                                                                             <  1:431174/95‑1 (MQ=255)
                                                                                                                                                                                                         gTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGAAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAAc                             >  1:55247/1‑251 (MQ=255)
                                                                                                                                                                                                          tACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCATTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGCAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGACGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTCGACATTCCCCCACCCATCAAACCCAGGCCCCCAAATTTTACCCCAggggaaaaacc                            >  2:363447/1‑243 (MQ=255)
                                                                                                                                                                                                            cAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAACCAATCTTGATGTCTGGCTGCGGAAGCGCTTCACCTTCGATGTCTTTATTAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCCCCGTTCAGACCCAGGACACCCATTTTTTCCCCCGGGGAGAAACTcc                          >  2:329570/1‑250 (MQ=255)
                                                                                                                                                                                                                  aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGGCTTTTTCAATGCCCGCCATAATGCGCACCAGGGTAGGCTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCAAGGGAAGAAACTCcgag                       >  1:275853/1‑248 (MQ=255)
                                                                                                                                                                                                                  aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAAACAATCTTGATGTCTGGCTGCGGACGCGCCTCACCTTCGGTGGCTTTATCAATGCCCGCCATAATGCGCCGCAGGGTAGACTTACCCCCGCCGTTCAGGCCCAGGACACCAATTTTTGTCCCAGGGAAGaaaatccgagggg                    >  2:342870/1‑241 (MQ=255)
                                                                                                                                                                                                                       cTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt               >  1:141965/1‑251 (MQ=255)
                                                                                                                                                                                                                         ccATCCAGGCGTTTAAGGGCGTTAACCCCTACAGAACCCGCTTCTTCATGGGCCTCACGCGCGGTGTGTTCGGGTTTCGGCTGCGGTTCCTGCGGCAGATAACCATTCTTGAGGTCTGGCTGCGGACGCGCTTCACCTTGGATGTTTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGtttt              <  1:328420/249‑1 (MQ=255)
                                                                                                                                                                                                                         ccATCCAGGCGTTACAGGGCGTAAACCGCTTAAGAAACCTCTTCTTCAGGGGCCTCACGCACGGTGTGTCGCGGGTTCAGCTGGGGTTCCTGGGGCAGGTAACCAATCTGTATGTCGGGCTGCGGCCGCGCTTCACCTTCGTTGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt               <  1:367624/248‑1 (MQ=255)
                                                                                                                                                                                                                                    tttCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCCATGCCCGCCATAATGCgcagca                                                                                         >  1:517519/1‑164 (MQ=255)
                                                                                                                                                                                                                                     ttcggggCGTAAACCACTTAAGAACCCGTTTCTTCAATGGACTCCCGCCGGGTTTGTTCCGGTTTCAGCTGCGTTCCCGGCGGCAGGTAACAAATTTTGGTGTCGGGCTGCGGACGCGCTCCCCCTCCGATGTCTTTATCAAGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGAc  <  1:631452/246‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      
GCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGAC  >  NC_000913/4629991‑4630467

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