Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,984 G→T T84T (ACG→ACT slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,9840GT94.4% 144.1 / ‑2.2 56T84T (ACG→ACTsltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (19/32);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

AAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCC  >  NC_000913/4630738‑4631227
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
aaaaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCa                                                                                                                                                                                                                                                    <  1:330521/251‑1 (MQ=255)
 aaaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCAc                                                                                                                                                                                                                                                   >  1:73368/1‑251 (MQ=255)
  aaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGCCCTGTTGAATCAACCGGCGGGGACTGTCAcc                                                                                                                                                                                                                                                  >  2:322175/1‑250 (MQ=255)
                  cTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCt                                                                                                                                                                                                                                   >  1:462772/1‑250 (MQ=255)
                    ggcggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAc                                                                                                                                                                                                                                >  1:269812/1‑251 (MQ=255)
                     ggggcTGTTGGCTGCCGGGGTCTGTCTGTTGACGGTCAGAAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAcc                                                                                                                                                                                                                               <  2:330387/249‑1 (MQ=255)
                       ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGAGGTGACAGTAACTAACTTTTTTCGCCCTAAcccc                                                                                                                                                                                                                             >  1:107477/1‑251 (MQ=255)
                           gTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACg                                                                                                                                                                                                                          <  2:133864/250‑1 (MQ=255)
                             tGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCCCTAACTTTGTTCTCGCTAACCCCACgatt                                                                                                                                                                                                                       >  2:137657/1‑248 (MQ=255)
                             gggCTGCCGGTGTCTGGCGGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCtt                                                                                                                                                                                                                       <  2:361529/250‑1 (MQ=255)
                                             tctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGc                                                                                                                                                                                                       >  1:235764/1‑251 (MQ=255)
                                                    aCGGTCAACAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGACGTGGTCGAACAAATGATGCCTGGACTGAAGGGTTATCCGCTTTATCCCTACCTGGAAAACCGCCAGATCACCGATGATCAGATGAATCAACCGGCGGTGACGGTCAATAACTTTGATCACGCTAACCCCACGCTTCCTCCCGCTCGCCCGCTGCAAtc                                                                                                                                                                                                >  1:111233/1‑251 (MQ=255)
                                                      ggTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                                                                                                                                                                              >  1:155982/1‑251 (MQ=255)
                                                       gTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCg                                                                                                                                                                                             >  1:211951/1‑251 (MQ=255)
                                                              gagTGGCGCGAGCCGAGCCACTGGAGGTGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTTAAAGATTATCCTCTTTACCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACATTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTc                                                                                                                                                                                      <  2:506566/249‑1 (MQ=255)
                                                                       gAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTg                                                                                                                                                                             >  2:21120/1‑251 (MQ=255)
                                                                            gACTCACTGGATGAGCAGCGTAGTCGTTATGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGTCGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTccc                                                                                                                                                                                                                 >  2:23962/1‑210 (MQ=255)
                                                                            gACTCACTGGATGAGCAGCGTAGTCGTTATGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTccc                                                                                                                                                                                                                 <  1:23962/210‑1 (MQ=255)
                                                                            gACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGAGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGAGACCGTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCTCCCGCCCGCCCGCTGCAATCCCGTTTCGTCAATggaacggcgcg                                                                                                                                                                        >  2:462864/1‑241 (MQ=255)
                                                                                         agcagcGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTgg                                                                                                                                                           <  1:199320/251‑1 (MQ=255)
                                                                                                  gccGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGGCTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTg                                                                                                                                                  <  2:374047/249‑1 (MQ=255)
                                                                                                           cgcAAATCAAGCAGGCCTGGGATAAGCGACAAATGGATGTGTTCGAACAAATGATGCCTGGACTGAAAGATTATCCCCTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAGCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCttt                                                                                                                                        <  2:281841/251‑1 (MQ=255)
                                                                                                            gcAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTa                                                                                                                                       <  1:166867/251‑1 (MQ=255)
                                                                                                                      cAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGaaaa                                                                                                                             <  1:153390/251‑1 (MQ=255)
                                                                                                                         gCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGcc                                                                                                                         <  1:420969/251‑1 (MQ=255)
                                                                                                                          ccTGGGATAATCGACAAATGGATGTGGTCGAACAAACGATGCCTGGACTGAAAGATTATCCACTTTAGCCCTCCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCTCGCTGCCTCCCGCTCGCACGCTGCAATCGCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGccc                                                                                                                        <  2:204057/251‑1 (MQ=255)
                                                                                                                            gtggagaaTCGACAAATGGATGTGGTCGAACAAAGGATGCCTGTACTGAAAGATTAGCCACTTTATCCCTACCTGTAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCgg                                                                                                                      <  2:258358/245‑1 (MQ=255)
                                                                                                                                aTAATCGACAAATGGATGTGGTCGAACAACTGATGCCTGGACTGAAAGATTATCCACTTTCTCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTTACTTTGTTCGAGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCCATGAACTGGCGCGGCGGGAAGGCTGGCGTGGCTTGTTTGGCCTTTTCCCGGGAAAGCCCGGagcc                                                                                                                  >  2:493050/1‑248 (MQ=255)
                                                                                                                                aTAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACt                                                                                                                  <  1:493050/251‑1 (MQ=255)
                                                                                                                                     cGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCg                                                                                                              <  1:137657/250‑1 (MQ=255)
                                                                                                                                         aaaTGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGGATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAgcg                                                                                                         <  2:172334/251‑1 (MQ=255)
                                                                                                                                                      cGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCGGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATtact                                                                                           <  2:53257/251‑1 (MQ=255)
                                                                                                                                                      cGAACAAAAGATGCCTGGACTGAAAGATTATCCCCTTTATCCCTTCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATtact                                                                                           <  2:267482/251‑1 (MQ=255)
                                                                                                                                                                 tgCCAGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTATCCCCACGCTTCCCCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGcaaa                                                                                 <  2:442337/250‑4 (MQ=255)
                                                                                                                                                                 ggCCTGGACTGAAAGATTAGCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCTCTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACGTGCGCGGCGTGAATACTGGCGTGGCTTGTTAGTCTTGAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCaaaa                                                                                <  1:62715/250‑1 (MQ=255)
                                                                                                                                                                   ccTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAACCCGGAACTACCGGAGCGCCAATTAATTACTACTATGCAAAATg                                                                              >  1:63547/1‑251 (MQ=255)
                                                                                                                                                                    cTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGAGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGGCTTTAGCCCGGAAAAGGCCAGCACTACCGGAGCGCACTGTAATTACTACTATGCCAAATgg                                                                             >  2:40670/1‑251 (MQ=255)
                                                                                                                                                                       gACTGAAAGATTATCTACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTACCCCCACGCTTCCGCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGTCTGGCGTGGCTTGTTAGCCTTTAGCTCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAac                                                                          <  2:455632/251‑1 (MQ=255)
                                                                                                                                                                       gACTGAAAGATTATCCCCTTTATCCCTACCTGGAATTCCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAAGGGAac                                                                          <  2:455641/251‑1 (MQ=255)
                                                                                                                                                                       gACTGAAAGATTATCCACTTTATCCCTACCTTGAATACCGTCAGATCCCCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACCTTGTTCGCGCTAACCCCACGCTTCCGCTCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAac                                                                          <  1:364007/251‑1 (MQ=255)
                                                                                                                                                                       gACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTTGCCCGGAAAAACCCCGGAATACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAac                                                                          >  2:177446/1‑251 (MQ=255)
                                                                                                                                                                          ggaaagATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGTCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACAcc                                                                       <  2:471795/246‑1 (MQ=255)
                                                                                                                                                                                aTTATCCACTTTATCCCTACCTGGCATACAGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGACCTGGCGCGGCGAGAAGACTGGCGTGGCTTGTTAGACTTTAGACCGGAAAAGCCCGGGACTACCGAAGCGGAAAGTAATTAATACTATGCAAAATGGAAAACCgggg                                                                   >  1:357877/1‑248 (MQ=255)
                                                                                                                                                                                          ttATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTgaagaa                                                        <  1:262683/250‑1 (MQ=255)
                                                                                                                                                                                           tATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGaagaag                                                       <  1:21120/250‑1 (MQ=255)
                                                                                                                                                                                             cccctccctGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACGACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCt                                                    <  2:452649/250‑1 (MQ=255)
                                                                                                                                                                                                    ccggaaacacttcAGATCACCTATGACCTGAGGAACCAACCGGCGGTGCCTGTCACTAACTTTGTTCGCGCTAACCCACCGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATTAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAgg                                             <  2:119721/239‑1 (MQ=255)
                                                                                                                                                                                                                cAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATg                                                                                                   >  1:477644/1‑185 (MQ=255)
                                                                                                                                                                                                                cAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATg                                                                                                   <  2:477644/185‑1 (MQ=255)
                                                                                                                                                                                                                 aGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTAGCCCGGGAAAACCCGGAACTACCGAAGCGCAATGGAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAAAAGCCTGGCCAGGGGCGAAagaggtg                                >  1:270689/1‑248 (MQ=255)
                                                                                                                                                                                                                           tgaccTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTAACGAAGCGCAAAGTAATTACTACTATGCAAAATGGAACACCGGGACGAGGGAAAAAGCCTGGCACGGGGCGAAAGAGCAGTGGCCAACCg                      >  1:494586/5‑251 (MQ=255)
                                                                                                                                                                                                                                tgatgaATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACACTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAAATTGCGCGGCGTGAAGACTGGCGTTGCTTGTTAGCCTTTAACCCGGAAAAGCCCGGAACTACCCAAGCGCAATGTAATTTCAACTTTGTAAAATGGAACACCGGAGAGAGTGAAAAAGCCCTGCAAGGGGGGaaaaacattgggcaaacctgcca                  >  1:191645/1‑230 (MQ=255)
                                                                                                                                                                                                                                 gatgaATCAACCGGCGGTGACTGTCACTCACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGGAAAGCCCGGAACTAACGAAGCGCAATGTAATTACTAATCAGCACAATGGAACACCGGGCAGAGTGAAGACACCTGGCAAGGGGCGAAAGAGCTGTGGGTAACCgggcaa                 >  1:196000/1‑246 (MQ=255)
                                                                                                                                                                                                                                  atgaATCAACCGGCGGTGACTGTCACTAACGTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAagag               <  1:32720/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cggcggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAAc     <  1:179769/251‑1 (MQ=255)
                                                                                                                                                                                                                                               cggTGACTGTCACTATCTTTGTTCGCTCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCTTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACGATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCATGGGGCGAAAGAGCTGTGGCTAACCGGCAAGATCCCGCCTAACGcc  <  2:499609/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
AAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCC  >  NC_000913/4630738‑4631227

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: