Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,631,146 G→A A138A (GCG→GCA)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,1460GA97.5% 118.2 / NA 40A138A (GCG→GCA) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/22);  minor base C (1/0);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCA  >  NC_000913/4630913‑4631387
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
gATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGcc                                                                                                                                                                                                                                  >  1:131715/1‑250 (MQ=255)
      ccacTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAAGTACTAGTATGCAAAATGGATCACC‑GTGCAGAGTg                                                                                                                                                                                                                             <  2:27585/247‑1 (MQ=255)
       cccTTTATCCCTACCTGGAATACCGTCAGATTACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAAGGGAACACC‑GGGCAGAGTgaa                                                                                                                                                                                                                           <  2:26910/249‑1 (MQ=255)
               cccTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCt                                                                                                                                                                                                                    <  1:72570/250‑1 (MQ=255)
               cccTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTg                                                                                                                                                                                                                   <  2:215121/251‑1 (MQ=255)
                  tACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGCCTGTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCTCCCACTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑GCCGGAAA‑AGCCCGGAACTACCGGAGCCCAAACTAATTACTACTTAGCCAAATGGAAAAAC‑GGGCAAGGGGAAGAAGCCTGGcc                                                                                                                                                                                                                >  1:115941/1‑250 (MQ=255)
                     cTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTAACGAAGCGCAATGTAATTTATACTATGCAAAAAGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAgg                                                                                                                                                                                                             >  1:115434/1‑251 (MQ=255)
                     cTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGGACTTCCGAAGCCCCATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGCAGAAGACTgggcagg                                                                                                                                                                                                             >  2:96552/1‑246 (MQ=255)
                         aaTACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGAGTGATGACTGACGTGGCTTGTTAGGCTTTAGGCCCGGAAA‑‑GCCCCGAACTCCCGAAGCGCAAAGCAATTACCACTCTGCAAAATGGAACACA‑GGGAAGAGTTAAGAAGCCTGGCCAGGGGCg                                                                                                                                                                                                         >  1:23524/1‑251 (MQ=255)
                                 cAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAATGGGCGAAAGAGCt                                                                                                                                                                                                 <  2:164754/251‑1 (MQ=255)
                                                          aaCCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCt                                                                                                                                                                        <  1:202591/251‑1 (MQ=255)
                                                              ggcggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAAACGGCAAGAGCCAGCCTAACg                                                                                                                                                                    >  1:193259/1‑251 (MQ=255)
                                                              ggcggTGACTGTCACTAACTTTGTTCGCAATAACCCCACGCTTCCTCCCGCTCGCAAGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAACAACCCG‑ACCTACCGAAGCGCAATGTAATTACTCCTATGCAACATGGAACACC‑GGGCACAGTGAAGAAGACTGGCAAGGGGGGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACg                                                                                                                                                                    >  2:12473/1‑251 (MQ=255)
                                                                  gTGACTGTCACTAACGTTGTTCGCGCTAACCCCACGCTTCCTCTCGCTCGCACGCTGCAATCTCGTTTCGTCAATGACCTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTAG‑CCCGGAAA‑AGCCCGGATCTACCGAAGCGCAATGTAATGCCTACTATGCAAAATGTAACACC‑GGGCAGAGGGAAGAAGCCTGGCAAGGGGCGAAAGCGCTGTGGCTAACCGGCAAGAGACAGTCTAACGCCtg                                                                                                                                                                <  2:207464/251‑1 (MQ=255)
                                                                     actgTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCtgtgg                                                                                                                                                             >  2:61646/4‑250 (MQ=255)
                                                                         tCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCACTGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGACTTTAT‑CCCGGAAA‑AGCCCGGAACTACCGAAACGCAATGTGATTACTACTATGCAAAATGGAACACCCGGGC‑GAGCGAAGAAGCCTGGCAAGGGGCGAAAGGGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCtgtgcca                                                                                                                                                           >  1:1679/1‑246 (MQ=255)
                                                                         tCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACa                                                                                                                                                           >  1:92138/1‑249 (MQ=255)
                                                                                        gcgcTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCTTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAGGCGCAATTTAATTACTACTAGGCAAAATGGAACACC‑GGGAAGGGGGAAGAAGCCTGGCAAGGGGCCAAAGAGCTGTGGGTAACCGGCAAGAGCCAGCCTAAAGGCTGTGACAAGGTATTTAGGGTCTgg                                                                                                                                          >  2:200962/1‑251 (MQ=255)
                                                                                          gcTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGc                                                                                                                                                                                    <  1:5987/207‑1 (MQ=255)
                                                                                                     gCTTCCTCCCGCTCGAACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGCAGCCTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGCACTACCGAAGCGCAATGTTATTACTACTATGCAAAATGGA‑CACCCGGGCAGCGTTGAGAAGCCTGGCCAGGGGCGAAAGAGCTGTGGCTAACCGGCACGAGCCAGCCTAACGCCTGTGACCAGTTATTTAGTGGCTGGCGTGCGTCGGGt                                                                                                                              >  2:134505/1‑250 (MQ=255)
                                                                                                            cccGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGa                                                                                                                      <  1:114333/251‑1 (MQ=255)
                                                                                                                         tGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAGTGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGGCAGGGAAACAAGATCCGCTGGCGTAt                                                                                                         >  1:228635/1‑251 (MQ=255)
                                                                                                                          gCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTAt                                                                                                         <  1:141140/250‑1 (MQ=255)
                                                                                                                          gCAATCGCGGTTCGGCAATGAACGGGCGGGGAGTGAAGACTGGCGTGGCTTGTGAGCCTTTAG‑CCCGTAAA‑AGCCCTGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GTGCAGAGTTAAGAAGCCTGGCAATGGGCGAAAGAGCTGTTGCTAACCGGCAAGAGCCAGCTTAACGTCTGTGACTATTTAGTTATTGGGTGGCGTGCGTCAGGTAAACAAGATCCGCTTGAGTAtt                                                                                                        <  2:243837/251‑1 (MQ=255)
                                                                                                                           cAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTAttt                                                                                                       <  2:115434/251‑1 (MQ=255)
                                                                                                                                             gAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACGATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGc                                                                                     <  1:12473/251‑1 (MQ=255)
                                                                                                                                                     gcggcgTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAg                                                                             <  1:181235/251‑1 (MQ=255)
                                                                                                                                                         cgTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGGAATTACTACTATGCAAAATGGAACCCC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGGAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCCGTGACAAGTTATTTAGAGGCAGGCGTGCGGCAGGTAAACCAGATACGCTGGCGAATTTAGCGCGTATCCGACTGGCGGTGAAAGAggg                                                                         >  2:192684/1‑251 (MQ=255)
                                                                                                                                                                     cGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCt                                                             <  1:200962/251‑1 (MQ=255)
                                                                                                                                                                            tGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAAc                                                       >  1:155125/1‑250 (MQ=255)
                                                                                                                                                                                     ttAG‑CCCGGAAA‑AGCACGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGTTAACAGTGCTgg                                               <  1:62728/249‑1 (MQ=255)
                                                                                                                                                                                       aG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTTCTggcagg                                           <  2:131715/251‑1 (MQ=255)
                                                                                                                                                                                              gAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAt                                     <  1:62177/251‑1 (MQ=255)
                                                                                                                                                                                                aa‑aGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGAAAAATGGAACACC‑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)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
GATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG‑CCCGGAAA‑AGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACC‑GGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCA  >  NC_000913/4630913‑4631387

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