Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,357,202 C→T intergenic (‑92/+221) sapA ← / ← ymjA putative periplasmic binding protein SapA/DUF2543 domain‑containing protein YmjA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,357,2020CT97.4% 121.2 / NA 39intergenic (‑92/+221)sapA/ymjAputative periplasmic binding protein SapA/DUF2543 domain‑containing protein YmjA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/21);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGCAATTAAATCATTACTCATTACCCCATTGATTCAA  >  NC_000913/1356969‑1357446
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
tAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTg                                                                                                                                                                                                                                      >  1:53590/1‑250 (MQ=255)
 aaaGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCCGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGa                                                                                                                                                                                                                                     >  2:249775/1‑250 (MQ=255)
         tGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGa                                                                                                                                                                                                                             <  1:249775/250‑1 (MQ=255)
         tGACTTGCCAGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGTATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGa                                                                                                                                                                                                                             <  2:113508/250‑1 (MQ=255)
                          cAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAACGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAAc                                                                                                                                                                                                           >  1:343680/1‑251 (MQ=255)
                          cAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAAc                                                                                                                                                                                                           >  1:407521/1‑251 (MQ=255)
                                aCAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGt                                                                                                                                                                                                     <  2:100766/251‑1 (MQ=255)
                                 aaaaaCCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTAGGTAAAGAGGCTATCTTACTAATATTTAATGGCAGTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTGGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTa                                                                                                                                                                                                    <  1:406159/250‑1 (MQ=255)
                                            ccgctgatttcggCATGCGGGGGATAGTCCGTCGCGGCGATTTCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATTCCTTGCGCATAATCGTAAGGGCTTTTATGGAAATAGTCTATCTTACTAATATTTAATGCCTTTTGCCATTACCTTTTGTGTTCAGGGGTAGTGATGAGGTCTGTGTGGTTGGTCTTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTTCCCGTGATTtc                                                                                                                                                                                         <  2:89379/239‑1 (MQ=255)
                                                    atCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                                                                                                                 >  1:327328/1‑251 (MQ=255)
                                                      cAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctgtaacta                                                                                                                                                                               <  1:124189/251‑1 (MQ=255)
                                                          aTGCGGGGGAGATTCGGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCtataactataac                                                                                                                                                                           >  1:38728/1‑251 (MQ=255)
                                                                                 cgATTACCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGCTAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTACTGACATTTGCAATTACCGTTTGTGTTCAGGGGTCTTAATGAGGGCTGTGTGGTTGGTCGTTGAAAAGGAGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgcgg                                                                                                                                                      >  2:362222/1‑246 (MQ=255)
                                                                                       ccggACCACGCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGTCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGGTCCCGTGATTTCTCGATAACTATAACGCACAGAACAACTTAGCGAGGTGCAGAGc                                                                                                                                              <  2:45529/248‑1 (MQ=255)
                                                                                          gACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGt                                                                                                                                           <  2:52597/251‑1 (MQ=255)
                                                                                            tcccTCACACGTCCCGCAATCACCAACATAGCCGTTAATACTTGGCGCCTAATATTAAGGGATTTTATGTAAAAAGGCTATCTTCCTAATATTTAATGACATTTGCCATTACCGCTTGTGTTCAGGGGTCGTAGTGAGGGCTGTGAGGTTCGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTTCCCGTGATTTCTCTTTAACTATAACTCAACGAAAAACTTAGCGAGGAGCAGAGCgtta                                                                                                                                          >  2:290063/4‑247 (MQ=255)
                                                                                            ccACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGACGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTAc                                                                                                                                         >  2:290096/1‑251 (MQ=255)
                                                                                             cACTCACCAGACCAGAAATCACAAAAAGAGACGATAATAACTGGCGCATAATATTAAGGGATTTAATGTAAAGAGGCAATGTTACTAATATTTAATGACAGTTGACATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGACGAGAGGTTGGTCGTTGAACAGGAGACGTCGCCATCTGTTCAGCATACGTGCCCGCGATTTCTCTATAAATCTAACTCACAGTACACCTTAGCGAGACGCAGAGCGGAACg                                                                                                                                        >  2:157238/1‑251 (MQ=255)
                                                                                                 cacaAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATc                                                                                                                                     <  1:414875/250‑1 (MQ=255)
                                                                                                            cAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTgaga                                                                                                                               >  2:303952/1‑245 (MQ=255)
                                                                                                            cAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTgaga                                                                                                                               <  1:303952/245‑1 (MQ=255)
                                                                                                                cACCAAAAGAGACGATAATTCCTGGCGCATAATATTAAGTGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGGCATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGTTGTGTGGTTGGTCGTTGAACAGGTGACGGCGCCATCTGTGCATCATTCGTACCCGTGATTTCTCTATAACTATAACGCACAGAACAACTTAGCGAGGATCAGATCGGTGCGATCGGGAGGAACCCATTa                                                                                                                      <  2:358879/250‑1 (MQ=255)
                                                                                                                cACCAAAAGAAACGTTAATACCTGGCGCCTAATATTCAGGGATTTTATGTAAAGCGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTa                                                                                                                      >  2:139326/1‑250 (MQ=255)
                                                                                                                              ataataCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTg                                                                                                       >  2:312147/1‑251 (MQ=255)
                                                                                                                                       ggCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCAtt                                                                                              <  2:86473/251‑1 (MQ=255)
                                                                                                                                                   tAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATg                                                                                  >  1:302154/1‑251 (MQ=255)
                                                                                                                                                           tttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTg                                                                           <  1:286529/250‑1 (MQ=255)
                                                                                                                                                                   aaGAGGCTATCTTACTAATATTTAATGACATTTGCCCGGACCGTTTGAGTTCAGAGGTCGTAATGAGGGATGTGTGGTTGGCCGTTGAACAGGTGACGTCGCCATCTGTTCAACAGTCGTACCAGTGATTTCTATCTAACTATAAATCAAAGCACAAATTAGAGAGGAGCAGAGCGGTACGATCGTGCGAAACCCACTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATTCGGCGTGAAAGCttt                                                                   >  2:66394/1‑250 (MQ=255)
                                                                                                                                                                       ggCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCg                                                              <  1:181069/251‑1 (MQ=255)
                                                                                                                                                                               tACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACa                                                          <  2:222443/247‑1 (MQ=255)
                                                                                                                                                                                    atatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTg                                                 <  1:276378/251‑1 (MQ=255)
                                                                                                                                                                                        ttAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGAcc                                              <  2:118226/250‑1 (MQ=255)
                                                                                                                                                                                             gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGt                                        <  1:87611/251‑1 (MQ=255)
                                                                                                                                                                                                  ttGCCATTACCGTTTGGGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGc                                    <  2:232536/250‑1 (MQ=255)
                                                                                                                                                                                                    gCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGCAAt                                 <  2:38728/251‑1 (MQ=255)
                                                                                                                                                                                                      cATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGAGTCTGACCGATCGTCTGCAAtt                                >  2:84372/1‑250 (MQ=255)
                                                                                                                                                                                                                       aGGGGACGCAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTAAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGCCAACGATACTGACCGATCGTCTGCAATTAAATCATTACTCATTAcc              >  2:269664/1‑251 (MQ=255)
                                                                                                                                                                                                                           gTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGGGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTAGCCGGACAACGATACTGGCCGATCGTCTGCAATTAAATCAGTACTCATTACCCCAt          <  2:162748/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGCACTTAAATCATTACTCATTACCCCATTGATTCaa  >  1:124019/1‑250 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
TAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGCAATTAAATCATTACTCATTACCCCATTGATTCAA  >  NC_000913/1356969‑1357446

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