Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,106 A→G N139N (AAT→AAC)  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,1060AG95.7% 205.8 / ‑0.9 70N139N (AAT→AAC) ettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (27/40);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

AGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCT‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGT  >  NC_000913/4629862‑4630350
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
aggaaCGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTGGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGa                                                                                                                                                                                                                                                 <  2:482962/248‑1 (MQ=255)
aGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGa                                                                                                                                                                                                                                                 <  1:448331/251‑1 (MQ=255)
     cGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGCGAGGTTAGCGCTTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCCGGTTCTGACcgtc                                                                                                                                                                                                                                            >  2:29368/1‑251 (MQ=255)
      ggTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGCAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCCTCCGCCGCACGCTCCAGCTGTACGTTCAGGTAAGGACcgtcg                                                                                                                                                                                                                                           >  1:6931/1‑251 (MQ=255)
      ggTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACcgtcg                                                                                                                                                                                                                                           >  1:44867/1‑251 (MQ=255)
          gTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGCTTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCGGGCGGCACAAC‑GCTACGCGCCGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGACCGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTg                                                                                                                                                                                                                                         >  1:146152/1‑249 (MQ=255)
                 aTCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAAt                                                                                                                                                                                                                                <  1:352222/251‑1 (MQ=255)
                          ccGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCTGCCAGGTGGTTGGTCGGTTCGCCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTc                                                                                                                                                                                                                       <  2:273595/250‑1 (MQ=255)
                            ggTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCa                                                                                                                                                                                                                     <  1:184762/251‑1 (MQ=255)
                                     gAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCtt                                                                                                                                                                                                            >  1:356292/1‑251 (MQ=255)
                                       aGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCCTGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGt                                                                                                                                                                                                          >  1:453283/1‑251 (MQ=255)
                                          cGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTGCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTc                                                                                                                                                                                                        <  2:250095/250‑1 (MQ=255)
                                           gTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCCCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGCATGATCTCTTCCAGCCGGCCTTTTTccg                                                                                                                                                                                                      >  2:235584/1‑249 (MQ=255)
                                              cAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCTAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGCCGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCgg                                                                                                                                                                                                   >  1:206908/1‑251 (MQ=255)
                                               aGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGc                                                                                                                                                                                                  >  1:39596/1‑251 (MQ=255)
                                                  aaGCGTTCCAGCCAGGCCACGGTTTCGGCAGCCAGGTGGTTGTGCGGTTCGGCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAg                                                                                                                                                                                               <  2:195615/251‑1 (MQ=255)
                                                      gTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGCGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTTATCAGCGGccgacttt                                                                                                                                                                                           >  2:386209/1‑245 (MQ=255)
                                                                  ccACGGGTTCGGCAGCCAGGTGGGTGGTCGGTTCGTCGAGCAGAAGCATGTCTGGTTTTTCCAGAAGCAGGCGGCCAAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCa                                                                                                                                                                               <  2:453283/251‑1 (MQ=255)
                                                                   cACGGTTTCGGCATCCAGGGGTTTGGTCGTTGCGTCGAGCAGAAGCATTCCTGGTTTTTCCAGCAGCAGGCGGCCCAAC‑GCTTCGCGACGACGCCCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCGTGTTCAGCGGCCAGCTTGTCAAAATCGGCAt                                                                                                                                                                              <  2:294896/251‑1 (MQ=255)
                                                                    ccggtgtCGGCATCCAGGTGGTTGTGCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTTCGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATc                                                                                                                                                                             <  2:391252/245‑1 (MQ=255)
                                                                      ggATTCGGCACCCAGGGGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGGCGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGAGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCg                                                                                                                                                                           <  2:191428/251‑1 (MQ=255)
                                                                       gATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCgg                                                                                                                                                                          >  1:264525/1‑251 (MQ=255)
                                                                           cGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGCTTTTCGCGTCCCAGTCCGGCAGAAGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGCGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCg                                                                                                                                                                      >  1:290784/1‑251 (MQ=255)
                                                                           cGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCg                                                                                                                                                                      >  1:149701/1‑251 (MQ=255)
                                                                                 ccatgtggtTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTa                                                                                                                                                                 <  1:40601/246‑1 (MQ=255)
                                                                                  cAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCCAAATCGGCATCCGGATCGGCGTAc                                                                                                                                                                >  1:364819/1‑250 (MQ=255)
                                                                                    ggtggtTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGTTCACCACCG‑AGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGAC‑‑AGCGCATCCGCCGCACGCTCCCGCTCTACGATCAGATTATGACCGTCGTGTGCCTGCATGCTCTCTTCCAGCCGCCCTTGTTCAGCGGCCAGCTTGTCACAATCGTCATCCGGATCGGCGttcctcccc                                                                                                                                                          >  2:114207/1‑243 (MQ=255)
                                                                                     gtggtTGGTCGGTTCGTAGAGCAGCCGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGGCGTTCACCACCGGAGAGGTTAGCGCTTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCCGACTGATCTCTTCCAGACGGCCTTGTTCAGAGGCCAGCTTGTCAACATAGGCATCCGGATCGGCGtcctgcg                                                                                                                                                            >  2:38275/1‑245 (MQ=255)
                                                                                       ggttggtCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAGTGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCTTACAgcgc                                                                                                                                                           >  1:67634/1‑250 (MQ=255)
                                                                                          tggtCGGTTCGTCGGGCAGCAGCATGCCTGGTTGTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGCTAGCGATTTTCGCGTCCCAGTCCGGCAGACGGAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATa                                                                                                                                                        <  2:54431/250‑1 (MQ=255)
                                                                                          tggtCGGTTCGTCGAGCAGCAGACGGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGAGACGACGCTCACCACCGGAGAGGTTAGCTATTTTCGCGTCCCAGTCCGGCAGACGGAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTGCAGCGGCCAGCTTGGCAAAATCGGCATCCGGATCGGCGTACAGCGCATa                                                                                                                                                        <  2:398856/250‑1 (MQ=255)
                                                                                                 ttCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCCTCCGGATCGGCGTACCGCGCATACACTTCAt                                                                                                                                                >  2:260506/1‑251 (MQ=255)
                                                                                                  tcgcggaggagcggcaTGGCTGGTTTTCCCAGCAGGCGGCGGCACAAC‑GCTTCGCGACGACGCTCCCCACCGGAGAGGTTAGCGAGTTGCGCGTCCCAGCCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCGTGAATGATCTCTTCCAGACGGCCGTGTTCAGCGGCCAGCTTGTCAAAAGCGGCAGCCGGATCGGCGTACAGCGCATactcgtcat                                                                                                                                                <  2:233067/237‑8 (MQ=255)
                                                                                                    gtcgAGCAGCAGCATGTCTGGGTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCGGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGTATCGGCGTACAGCGCATACACTTCATCCa                                                                                                                                             <  1:386209/251‑1 (MQ=255)
                                                                                                                   gTCTGTTTTTTCCAGGAGCGGGCGG‑ACAACCTCTACGCGACGCCGCCCACCCCCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTAAGACCGGCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCTGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCg                                                                                                                               <  2:41996/250‑1 (MQ=255)
                                                                                                                   gTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGt                                                                                                                              >  1:409157/1‑251 (MQ=255)
                                                                                                                       ggTTTTTCCTGCAGCGGGCGGCACAAC‑GCGACGCGCCGACTCTCCCCACCGGAGAGTTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGAAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTaa                                                                                                                           <  2:296444/250‑1 (MQ=255)
                                                                                                                       ggTTTTCCCAGCAGCAGGCGGCACAAC‑GCTACTCGACGACGCTCACAACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCTCATCCGCCGCACTCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAAc                                                                                                                          <  1:381376/251‑1 (MQ=255)
                                                                                                                              cccgcagatggcggcACAAC‑GCTACGCGTCGACGCTCACCACCGGAGAGGTTAGCGATTGTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTca                                                                                                                   <  2:478320/242‑2 (MQ=255)
                                                                                                                               cagcagcagGCGGCACAAC‑TCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTcag                                                                                                                  <  1:416181/251‑3 (MQ=255)
                                                                                                                                    gcagGCGGCACAAC‑GCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCT‑GAAAcc                                                                                                            <  1:112306/251‑1 (MQ=255)
                                                                                                                                    gcagGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCCGGCGTTTCAGGGAGTTAACCACTTCA‑GAAAcc                                                                                                            >  1:138832/1‑251 (MQ=255)
                                                                                                                                     cagGCGGCACAAC‑GCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGACCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTAAGGGCGTTACCCACTTcagcaac                                                                                                              >  2:21304/1‑243 (MQ=255)
                                                                                                                                          ggcACAAC‑GCTACGCGACGACGCTCACCACCGGCGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGcttct                                                                                                      >  1:194252/1‑251 (MQ=255)
                                                                                                                                            cacaAC‑GCTACGCGACCACGCTCCCCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCCCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCGTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCttcttc                                                                                                    <  2:223627/251‑1 (MQ=255)
                                                                                                                                                             gacgcccaccacCGGAGAGTTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTGCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTcacgcacg                                                                                    <  2:119039/245‑1 (MQ=255)
                                                                                                                                                                gcccaccacCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGgtg                                                                                 <  2:413141/248‑1 (MQ=255)
                                                                                                                                                                   caccacCGGGGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGCCGTGAGCCGGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGgtgtgt                                                                              <  2:501003/251‑1 (MQ=255)
                                                                                                                                                                     ccaccGGAGGGGTTTGCGTTTTTCTCGACCCAGTCCTGCATACGTAGCGAATCCGCCGCACGCTCCAGCTGGACGGTCAGGGTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGTTCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGgtgtgtgc                                                                            <  2:302951/251‑3 (MQ=255)
                                                                                                                                                                        ccGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGAGTGTTCCgg                                                                         >  1:5010/1‑251 (MQ=255)
                                                                                                                                                                            agagGTGAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTAAGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGCCGGACTTGTTCAGCGGCCAGCTTGTCAACATCGGCATCCGGCGCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGAGTGTTC‑GGGtttc                                                                    >  2:460676/1‑249 (MQ=255)
                                                                                                                                                                             gagGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCGCTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCa                                                                    <  1:194058/251‑1 (MQ=255)
                                                                                                                                                                                           gcgtgtccgataccggaagacgcgtcgcAGCCGCCGCACGCTCCAGCTTTACGGTAAGGTTTTGACCGTCGTGAGCCTGAATGATCTCGTCCAGACGGCATTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCTTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACT‑CTAGAAACCGCGTCTTCCATGGACTCACGCACGGTGTGGTCCGGGTTCAGCTGCGGTTCCTGc                                                      <  2:206876/226‑1 (MQ=255)
                                                                                                                                                                                             gcgTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCgg                                                    <  1:131761/251‑1 (MQ=255)
                                                                                                                                                                                               gcccccgTCCGGCAGACGTGCCGCATCCGCCGCCCGCCCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCCCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACGCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCa                                                  <  2:146152/245‑1 (MQ=255)
                                                                                                                                                                                                   cAGTCCGTCAGACGTAGCGCATCCGCCGCCCGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTTAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAggta                                              <  1:165094/251‑4 (MQ=255)
                                                                                                                                                                                                   cAGTCCGGCAGACGTAGCGCCTCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAggta                                              <  1:334214/251‑4 (MQ=255)
                                                                                                                                                                                                          gCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCacac                                       >  2:82740/1‑248 (MQ=255)
                                                                                                                                                                                                              aCGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGa                                   <  1:445951/251‑1 (MQ=255)
                                                                                                                                                                                                               cGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTCCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGGGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAgggcacccactttg                                    >  2:162622/1‑236 (MQ=255)
                                                                                                                                                                                                                  aGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCGGAATGATCTCTGCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGt                                <  1:44838/250‑1 (MQ=255)
                                                                                                                                                                                                                    cgcATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTg                             <  2:39596/251‑1 (MQ=255)
                                                                                                                                                                                                                      cATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCCGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTGACCAAACTTGATGtttgtc                           >  2:22853/1‑246 (MQ=255)
                                                                                                                                                                                                                       aTCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCt                          >  1:408631/1‑251 (MQ=255)
                                                                                                                                                                                                                           gatcgacgttCCAGCTGTCCGTTCAGGTTATGCCCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGGCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCgg                      <  2:416717/242‑1 (MQ=255)
                                                                                                                                                                                                                                aCGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGCAACCACTCTTTATGTCTGGCTGCGGaggtc                 >  1:416717/1‑247 (MQ=255)
                                                                                                                                                                                                                                         cTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAAt                                        <  1:391254/219‑1 (MQ=255)
                                                                                                                                                                                                                                             ccGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGAt    <  1:127803/250‑1 (MQ=255)
                                                                                                                                                                                                                                               gTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAGCCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGt  <  1:424949/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCA‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGt  <  1:235584/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
AGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAAC‑GCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCT‑GAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGT  >  NC_000913/4629862‑4630350

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