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

TCTGGCAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGT  >  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)
                                                           cGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCGCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCa                                                                                                                                                                                         <  2:24408/251‑1 (MQ=255)
                                                               cgcgGAATGTATCTGTATGAGCATCGTCGCTGATTCAATAACATCATCGCCTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAATGAGGGATGGAAACAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGAAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGGGGCAAT‑TAAGgatgattatgagctcccttac                                                                                                                                                                                     >  2:376521/1‑236 (MQ=255)
                                                                  ggAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTAACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCAATATGGTGGGAGCCGTCGTAGCCATCTTGACGCTGTAACAATAGCAGCATGGATGGAGG‑ACACTAAGGATGATTTAGATCGACAtttaggg                                                                                                                                                                                  >  2:221085/1‑246 (MQ=255)
                                                                               aTGAGCATCGGCGCTGGTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGTGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGCTGGTAGAGGCATGAAGCCAAATGGAGGGAGACGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATAAAGAGCGCCATTACGGGGAACAACTGTAc                                                                                                                                                                      >  2:324410/1‑250 (MQ=255)
                                                                                     aTCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAgttg                                                                                                                                                               >  2:27953/1‑251 (MQ=255)
                                                                                        gtcgCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGttgttg                                                                                                                                                            <  1:206439/251‑1 (MQ=255)
                                                                                        gtcgCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGttgttg                                                                                                                                                            <  1:474066/251‑1 (MQ=255)
                                                                                              gTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGtgtt                                                                                                                                                       <  1:488313/250‑1 (MQ=255)
                                                                                                         atcatcGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACgg                                                                                                                                                                                   >  2:65319/1‑211 (MQ=255)
                                                                                                         atcatcGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACgg                                                                                                                                                                                   <  1:65319/211‑1 (MQ=255)
                                                                                                                  tactGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGTTGTTGTGACGCCATTTTCGgcgt                                                                                                                                    >  1:111589/1‑248 (MQ=255)
                                                                                                                    ctGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGt                                                                                                                                >  1:443769/1‑251 (MQ=255)
                                                                                                                             acgttGCTATCACGGTAAGTTGCAGCCTAAAAGATGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTAt                                                                                                                       <  1:372792/248‑1 (MQ=255)
                                                                                                                                  gCTATTACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCaaa                                                                                                                   <  2:298680/250‑1 (MQ=255)
                                                                                                                                     aTCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGt                                                                                                               >  1:389022/1‑251 (MQ=255)
                                                                                                                                          ggTAAGTTGCAGCCTAAAAGTGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGCGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAACAGCATGGATGGTGGCCCA‑TACGGATGATATAGCGCGCCATTACGGGGAACCACCGTACCAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTCTCAAAATGTCATTg                                                                                                          >  2:152962/1‑251 (MQ=255)
                                                                                                                                           gTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGCTTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGaa                                                                                                                                                                                                                                         >  1:160420/1‑124 (MQ=255)
                                                                                                                                                     gCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGa                                                                                               <  1:358400/251‑1 (MQ=255)
                                                                                                                                                                ggATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGGTGATATAGAGCGCCATTACGGGGAACCACTGTCCAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGc                                                                                     >  2:131661/1‑250 (MQ=255)
                                                                                                                                                                      aaaaaGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACg                                                                              <  1:4804/251‑1 (MQ=255)
                                                                                                                                                                               gATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGCTGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGtt                                                                     >  1:308085/1‑251 (MQ=255)
                                                                                                                                                                                     ccccGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGa                                                               <  1:252865/251‑1 (MQ=255)
                                                                                                                                                                                           aaTGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCAt                                                                                                            >  2:92200/1‑200 (MQ=255)
                                                                                                                                                                                           aaTGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCAt                                                                                                            <  1:92200/200‑1 (MQ=255)
                                                                                                                                                                                                   atTTACAAATTGAGGGTAGATGCATGAAGCCAAGTGGAGGGAGCTGTCATAGCCCTCTTGAAGCGGAAACAATAGCTGCATTGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAgatgcttcact                                                  <  2:4504/250‑11 (MQ=255)
                                                                                                                                                                                                           atTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGGTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCtt                                         <  1:118460/251‑1 (MQ=255)
                                                                                                                                                                                                              gATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGc                                       >  1:413754/1‑250 (MQ=255)
                                                                                                                                                                                                              gATGGTAGAGGCAGGAAGCCAAATGGAGGGGGCCGTCATAGCCCGCTTGAAGCTGAAACAATAGCAGCATTGCGGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGACATTTGCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTTAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTAATGACAGTGTGGACTGTCCGCTCCtcgc                                       <  2:19447/250‑4 (MQ=255)
                                                                                                                                                                                                                        gCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGa                             >  1:294159/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCt    <  1:131661/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  cATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGt  >  1:183396/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
TCTGGCAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACA‑TAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGT  >  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: ‑