Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,631,245 C→T S171S (AGC→AGT)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,2450CT97.1% 105.0 / NA 36S171S (AGC→AGT) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/20);  minor base G (1/0);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTT  >  NC_000913/4631001‑4631478
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         
gcgcTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCTTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAGGCGCAATTTAATTACTACTAGGCAAAATGGAACACC‑GGGAAGGGGGAAGAAGCCTGGCAAGGGGCCAAAGAGCTGTGGGTAACCGGCAAGAGCCAGCCTAAAGGCTGTGACAAGGTATTTAGGGTCTgg                                                                                                                                                                                                                                     >  2:200962/1‑251 (MQ=255)
             gCTTCCTCCCGCTCGAACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGCAGCCTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGCACTACCGAAGCGCAATGTTATTACTACTATGCAAAATGGA‑CACCCGGGCAGCGTTGAGAAGCCTGGCCAGGGGCGAAAGAGCTGTGGCTAACCGGCACGAGCCAGCCTAACGCCTGTGACCAGTTATTTAGTGGCTGGCGTGCGTCGGGt                                                                                                                                                                                                                         >  2:134505/1‑250 (MQ=255)
                    cccGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGa                                                                                                                                                                                                                 <  1:114333/251‑1 (MQ=255)
                                 tGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAGTGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGGCAGGGAAACAAGATCCGCTGGCGTAt                                                                                                                                                                                                    >  1:228635/1‑251 (MQ=255)
                                  gCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTAt                                                                                                                                                                                                    <  1:141140/250‑1 (MQ=255)
                                  gCAATCGCGGTTCGGCAATGAACGGGCGGGGAGTGAAGACTGGCGTGGCTTGTGAGCCTTTAGCCCGTAAAAGCCCTGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GTGCAGAGTTAAGAAGCCTGGCAATGGGCGAAAGAGCTGTTGCTAACCGGCAAGAGCCAGCTTAACGTCTGTGACTATTTAGTTATTGGGTGGCGTGCGTCAGGTAAACAAGATCCGCTTGAGTAtt                                                                                                                                                                                                   <  2:243837/251‑1 (MQ=255)
                                   cAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTAttt                                                                                                                                                                                                  <  2:115434/251‑1 (MQ=255)
                                                     gAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACGATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGc                                                                                                                                                                                <  1:12473/251‑1 (MQ=255)
                                                             gcggcgTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAg                                                                                                                                                                        <  1:181235/251‑1 (MQ=255)
                                                                 cgTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGGAATTACTACTATGCAAAATGGAACCCC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGGAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCCGTGACAAGTTATTTAGAGGCAGGCGTGCGGCAGGTAAACCAGATACGCTGGCGAATTTAGCGCGTATCCGACTGGCGGTGAAAGAggg                                                                                                                                                                    >  2:192684/1‑251 (MQ=255)
                                                                             cGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCt                                                                                                                                                        <  1:200962/251‑1 (MQ=255)
                                                                                    tGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAAc                                                                                                                                                  >  1:155125/1‑250 (MQ=255)
                                                                                             ttAGCCCGGAAAAGCACGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGTTAACAGTGCTgg                                                                                                                                          <  1:62728/249‑1 (MQ=255)
                                                                                               aGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTTCTggcagg                                                                                                                                      <  2:131715/251‑1 (MQ=255)
                                                                                                     gAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAt                                                                                                                                <  1:62177/251‑1 (MQ=255)
                                                                                                       aaaGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGAAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAtgc                                                                                                                              >  1:238221/1‑251 (MQ=255)
                                                                                                             cGGAACTACCGAAGCGCAATGTAATTACTCCTATGCAAAATGGAACACC‑GGGTAGAGTGAAGAAGCCTGCCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGTCTGTGACTAGTGATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTACCACGGGCCTGGTAAGCGTGCTGGCAAGCCAGATGTCTGCCg                                                                                                                        <  2:109958/251‑1 (MQ=255)
                                                                                                                            gcAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCg                                                                                                         >  1:5711/1‑251 (MQ=255)
                                                                                                                              aaTGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGcc                                                                                                       >  1:87816/1‑251 (MQ=255)
                                                                                                                                 gtcaTTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGTCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGTCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGTTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTc                                                                                                     <  1:119124/247‑1 (MQ=255)
                                                                                                                                 gTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTc                                                                                                     <  1:47773/250‑1 (MQ=255)
                                                                                                                                          actaTGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCa                                                                                             <  2:193259/249‑1 (MQ=255)
                                                                                                                                                        aaCACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGATCTGTGGCTAACCGGCAAGAGCCAGCCTACCGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCATGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACGATCGCCTCGGCAATGATTTCACTTGCGaacaa                                                                             <  1:214266/251‑1 (MQ=255)
                                                                                                                                                             c‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCt                                                                         >  1:91611/1‑250 (MQ=255)
                                                                                                                                                                  cAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACTCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGGCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGTACAACCCTAATAc                                                                    <  2:94322/251‑1 (MQ=255)
                                                                                                                                                                  cAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAc                                                                    <  1:61646/251‑1 (MQ=255)
                                                                                                                                                                        gaagaaGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACATCCCTAATACGGTAc                                                               <  2:171434/250‑1 (MQ=255)
                                                                                                                                                                         aagaagCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCAGTTCACTGGCGAACAACCCTAATACGGTACTg                                                             <  1:96552/251‑1 (MQ=255)
                                                                                                                                                                            aagCCTGGCAATGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGACTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGAcc                                                          >  1:218498/1‑251 (MQ=255)
                                                                                                                                                                               ccTGGCAAGGGGTGAAAGAGCGGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGGTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCGCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTc                                                       <  2:87816/251‑1 (MQ=255)
                                                                                                                                                                                        gggCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGACGCGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTCTCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTCCTGACACTCGCGCGTAc                                               >  2:100389/1‑250 (MQ=255)
                                                                                                                                                                                            gAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTg                                          >  1:68684/1‑251 (MQ=255)
                                                                                                                                                                                              aaGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGgc                                        >  1:139698/1‑251 (MQ=255)
                                                                                                                                                                                                                    gagCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCaa                   >  1:53573/1‑250 (MQ=255)
                                                                                                                                                                                                                         agccTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATggcg             >  1:10544/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tgtgACAAGTTATTTAGTGTCTGGCGTGCGTCGGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTCCCACTGGCGCGACCGATTTTACCCGTCACATGGCGGCGGTGGCGtt  >  2:209578/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         
GCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTT  >  NC_000913/4631001‑4631478

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