Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,535,073 A→G intergenic (‑1020/‑1524) topAI ← / → ythB toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,535,0730AG94.9% 109.7 / ‑6.6 39intergenic (‑1020/‑1524)topAI/ythBtoxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (20/17);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TCTGGCAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGT  >  NC_000913/4534827‑4535319
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
tCTGACAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGATCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAATGAGGGATGGAAAAAGGACGATTTTCCCAGCTATGGGATATTTCCATATTGATGGTAGAGGCATGAAGCCACATGGAGGGAGCAGTCATAGCCct                                                                                                                                                                                                                                                     >  2:252865/1‑251 (MQ=255)
ccTGGCAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCct                                                                                                                                                                                                                                                     <  2:46097/250‑1 (MQ=255)
    gCAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACAAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGCCTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCATTGGGATATTTACATATTGATGGAAGAGGCATGAAGCCAAATGGAGGGAACCGTCATCACCCTCttt                                                                                                                                                                                                                                                 >  2:278542/1‑250 (MQ=255)
         gTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCt                                                                                                                                                                                                                                            <  2:142538/251‑1 (MQ=255)
                aGGTAGCATAGTTCAACTCTCCCTACCGGTATGATCCCAAGATCGAGCGCGCAATGTATCTGAATGCGCATCGTCGCTGTTTCAATAACATCATCGCCTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGCTGGACAAAGGACGATTTTACCCGCAATGGGATATTTCCATCTTGCTGATAGAGGCATGAAGCCAAATGGGGGGAGCAGTCATAACCCTCTTGAAGCTGAAACa                                                                                                                                                                                                                                      >  2:437525/1‑250 (MQ=255)
                     gCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATagc                                                                                                                                                                                                                                 >  1:477201/1‑250 (MQ=255)
                      cccAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATagc                                                                                                                                                                                                                                 >  2:488313/1‑249 (MQ=255)
                       ccAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATagcagc                                                                                                                                                                                                                              >  1:298680/1‑251 (MQ=255)
                                                           cGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCGCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCa                                                                                                                                                                                         <  2:24408/251‑1 (MQ=255)
                                                               cgcgGAATGTATCTGTATGAGCATCGTCGCTGATTCAATAACATCATCGCCTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAATGAGGGATGGAAACAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGAAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGGGGCAATTAAGgatgattatgagctcccttac                                                                                                                                                                                     >  2:376521/1‑236 (MQ=255)
                                                                  ggAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTAACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCAATATGGTGGGAGCCGTCGTAGCCATCTTGACGCTGTAACAATAGCAGCATGGATGGAGGACACTAAGGATGATTTAGATCGACAtttaggg                                                                                                                                                                                  >  2:221085/1‑246 (MQ=255)
                                                                               aTGAGCATCGGCGCTGGTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGTGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGCTGGTAGAGGCATGAAGCCAAATGGAGGGAGACGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATAAAGAGCGCCATTACGGGGAACAACTGTAc                                                                                                                                                                      >  2:324410/1‑250 (MQ=255)
                                                                                     aTCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAgttg                                                                                                                                                               >  2:27953/1‑251 (MQ=255)
                                                                                        gtcgCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGttgttg                                                                                                                                                            <  1:206439/251‑1 (MQ=255)
                                                                                        gtcgCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGttgttg                                                                                                                                                            <  1:474066/251‑1 (MQ=255)
                                                                                              gTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGtgtt                                                                                                                                                       <  1:488313/250‑1 (MQ=255)
                                                                                                         atcatcGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACgg                                                                                                                                                                                   >  2:65319/1‑211 (MQ=255)
                                                                                                         atcatcGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACgg                                                                                                                                                                                   <  1:65319/211‑1 (MQ=255)
                                                                                                                  tactGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGTTGTTGTGACGCCATTTTCGgcgt                                                                                                                                    >  1:111589/1‑248 (MQ=255)
                                                                                                                    ctGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGt                                                                                                                                >  1:443769/1‑251 (MQ=255)
                                                                                                                             acgttGCTATCACGGTAAGTTGCAGCCTAAAAGATGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTAt                                                                                                                       <  1:372792/248‑1 (MQ=255)
                                                                                                                                  gCTATTACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCaaa                                                                                                                   <  2:298680/250‑1 (MQ=255)
                                                                                                                                     aTCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGt                                                                                                               >  1:389022/1‑251 (MQ=255)
                                                                                                                                          ggTAAGTTGCAGCCTAAAAGTGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGCGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAACAGCATGGATGGTGGCCCATACGGATGATATAGCGCGCCATTACGGGGAACCACCGTACCAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTCTCAAAATGTCATTg                                                                                                          >  2:152962/1‑251 (MQ=255)
                                                                                                                                           gTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGCTTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGaa                                                                                                                                                                                                                                         >  1:160420/1‑124 (MQ=255)
                                                                                                                                                     gCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGa                                                                                               <  1:358400/251‑1 (MQ=255)
                                                                                                                                                                ggATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGGTGATATAGAGCGCCATTACGGGGAACCACTGTCCAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGc                                                                                     >  2:131661/1‑250 (MQ=255)
                                                                                                                                                                      aaaaaGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACg                                                                              <  1:4804/251‑1 (MQ=255)
                                                                                                                                                                               gATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGCTGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGtt                                                                     >  1:308085/1‑251 (MQ=255)
                                                                                                                                                                                     ccccGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGa                                                               <  1:252865/251‑1 (MQ=255)
                                                                                                                                                                                           aaTGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCAt                                                                                                            >  2:92200/1‑200 (MQ=255)
                                                                                                                                                                                           aaTGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCAt                                                                                                            <  1:92200/200‑1 (MQ=255)
                                                                                                                                                                                                   atTTACAAATTGAGGGTAGATGCATGAAGCCAAGTGGAGGGAGCTGTCATAGCCCTCTTGAAGCGGAAACAATAGCTGCATTGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAgatgcttcact                                                  <  2:4504/250‑11 (MQ=255)
                                                                                                                                                                                                           atTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGGTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCtt                                         <  1:118460/251‑1 (MQ=255)
                                                                                                                                                                                                              gATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGc                                       >  1:413754/1‑250 (MQ=255)
                                                                                                                                                                                                              gATGGTAGAGGCAGGAAGCCAAATGGAGGGGGCCGTCATAGCCCGCTTGAAGCTGAAACAATAGCAGCATTGCGGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGACATTTGCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTTAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTAATGACAGTGTGGACTGTCCGCTCCtcgc                                       <  2:19447/250‑4 (MQ=255)
                                                                                                                                                                                                                        gCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGa                             >  1:294159/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCt    <  1:131661/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  cATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGt  >  1:183396/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
TCTGGCAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGT  >  NC_000913/4534827‑4535319

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