Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,945 C→T R71R (CGC→CGT slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,9450CT95.1% 164.8 / ‑0.2 62R71R (CGC→CGTsltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (31/27);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

ATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCG  >  NC_000913/4630702‑4631188
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
aTTGATTTTTTTCCACTTTGAGTATGCGCTGGTGGAAAAAGCTAAAAAAGTTACCTGGCGGCTGGTGGCTGCCGGTGTCTGTCTGTTTACGTGCAGCATCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATAAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTAGCCCTACCTTGAAGACCGTCAGATCa                                                                                                                                                                                                                                                    <  2:477726/251‑1 (MQ=255)
aTTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTATGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAgctc                                                                                                                                                                                                                                                     >  1:627631/1‑247 (MQ=255)
aTTGATGTATTTACACTTAGAGGAGGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTTCCTGGAATACCGTCAGATCa                                                                                                                                                                                                                                                    <  1:464466/251‑1 (MQ=255)
   gATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACACACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCg                                                                                                                                                                                                                                                 >  2:71199/1‑251 (MQ=255)
                 tAGAGGATGCGCTTGGGGAAAAAGCCAAACAAGTTACCTGGCGGGTGTTGGCTGTCGGTGTCTGTCTGCTGACGGTCAGGAGCGTGGCGCGAGCCGACTCCCTGGATGAGCAGCGTAGCCGTTACGCGAAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATCCCGTCAGATCACCGATGACCtgatga                                                                                                                                                                                                                                     <  1:126402/249‑1 (MQ=255)
                  agagGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATc                                                                                                                                                                                                                                  >  1:559521/1‑251 (MQ=255)
                       aTGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACcg                                                                                                                                                                                                                             <  1:496349/251‑1 (MQ=255)
                              tgtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGACTGCTGACGGTCAGCAGCGTGGCGCGAGCCGCCTCACTGGATGAGCAGCGTAGCAGTTACGCGCAAATCATGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAAAGATGCCTGGACTGAAAGATTATCCACATTATCCCTACCTGGCATACCGTCAGATCACAGATGACCTGGTGAAACAAAcggcag                                                                                                                                                                                                                         >  2:366454/1‑246 (MQ=255)
                               gtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGACGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGGTCACCGATGACCTGATGAATCAACCGGCGGTg                                                                                                                                                                                                                       >  2:82426/1‑249 (MQ=255)
                               ggggaaaaccccaacaaaGTTACCTGGCGGGTGTGGGCGTCCGGTGTCGGTCTTCTGACGGTCAGCAGCTTGGAGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCACATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAATATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGAc                                                                                                                                                                                                                    <  2:589543/235‑1 (MQ=255)
                                 ggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGTCTGCCGGTGCCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGActg                                                                                                                                                                                                                  <  1:352323/251‑3 (MQ=255)
                                 ggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGActg                                                                                                                                                                                                                  <  1:461966/251‑3 (MQ=255)
                                              cAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGTCCCTAAATTTGt                                                                                                                                                                                                    >  2:20449/1‑251 (MQ=255)
                                               aaGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAACTCAAGCAGGCCTGGGATAATCGACAAATGTATGTGGTCGAACAAATGCTACCTGGACTGCAAGATTAACCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGTTGAATCAACCGGCGGTAAATGGTACTAACTTTGtt                                                                                                                                                                                                   >  2:353660/1‑251 (MQ=255)
                                                aGTTACCGGGCGGCTGTTGGCTGCCGGTGTCTGTTGGCTGACGGTCAGGAGCGTGTCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCTGCGGTTACTGTCACTAAGTTTGTTc                                                                                                                                                                                                  <  2:598542/251‑1 (MQ=255)
                                                       tGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCCCTAACTTTGTTCGCGCTa                                                                                                                                                                                            >  1:353216/1‑250 (MQ=255)
                                                         gcggcTGTTGGCTGCCGGTTTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgatga                                                                                                                                                                                                                                     <  2:364474/209‑1 (MQ=255)
                                                         gcggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgatga                                                                                                                                                                                                                                     >  1:364474/1‑209 (MQ=255)
                                                            gcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCa                                                                                                                                                                                      <  1:20449/251‑1 (MQ=255)
                                                              tGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGAACGTCACTAACTTTGTTGGCGCTTACCCCAcc                                                                                                                                                                                    >  2:373474/1‑250 (MQ=255)
                                                                 tGGCTGCGGGTGTCTGTATGCTGACGGTCCGCAGCGTGGCGCGAGACGCCTACCTGGATGAGCAGCGTAGCCGTTAAGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACACATGATGCCTGGACTGACAGGTTATCAAATTTATCCCTACCTGGAATACCGTCAGATCACCGATGCCCTGCTGAATACACCGGCGGTGTCTGTGCATAACTTTGTTCGCTCTAAACCCACGcat                                                                                                                                                                                 >  2:335787/1‑249 (MQ=255)
                                                                      gccgggggcggcctgctgACGGTCAGCAGCGTGGCGCGAGCCGACTCACTTGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGTGGACTGTCTCTAACTTTGCTCGCGCTAACCCCACGCTtccg                                                                                                                                                                              <  2:248744/237‑2 (MQ=255)
                                                                            gtctgtctGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACTCTTCCTCCCGCTCg                                                                                                                                                                      <  1:560023/251‑1 (MQ=255)
                                                                             tctgtctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGc                                                                                                                                                                     >  1:618920/1‑251 (MQ=255)
                                                                                     ctgACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGAGGTGACTGTCACTAACTTTGTTATCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGcc                                                                                                                                                             >  2:627777/1‑250 (MQ=255)
                                                                                             cagcagCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGCAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGCTGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGtt                                                                                                                                                     >  1:279981/1‑251 (MQ=255)
                                                                                                 agcGTGGCGCGAACCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACATGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTACCGCTCGCACGCTGCCATCTCGTTTcat                                                                                                                                                 >  1:318062/1‑249 (MQ=255)
                                                                                                        cgcgAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGCCTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGGCACTAACTTTGTTCACGCTAACCCCACGATTCCTCCCGCTCGAAAGCTGCAATCTCGTTTCGTCAATGa                                                                                                                                           >  1:1980/1‑250 (MQ=255)
                                                                                                            aGCCGACTCACTGGATGAGCCGCGTAGCCGTTACGCGCAAATCAAGCAGTCCTTGGATAAACGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTTAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCGCGTTTCGTCAATGAACTgg                                                                                                                                      <  2:87369/251‑1 (MQ=255)
                                                                                                               cGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTAACCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGgcgc                                                                                                                                   <  1:353660/251‑1 (MQ=255)
                                                                                                                     aCTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTg                                                                                                                             <  1:510282/251‑1 (MQ=255)
                                                                                                                            gAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCATGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTg                                                                                                                      >  1:223674/1‑251 (MQ=255)
                                                                                                                               cagcGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCg                                                                                                                   <  1:373474/251‑1 (MQ=255)
                                                                                                                                      gccGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTg                                                                                                            >  1:412412/3‑251 (MQ=255)
                                                                                                                                       ccGTTACGCTCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTTGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCTCCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGt                                                                                                           <  2:321099/250‑1 (MQ=255)
                                                                                                                                          ttACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGGGGCtttttct                                                                                                        >  2:533931/1‑246 (MQ=255)
                                                                                                                                                   aaTCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGcc                                                                                             <  1:660107/251‑1 (MQ=255)
                                                                                                                                                     ccAAGCAGGCCTGGGATAATCGACAAATGGATGTGGGCGAACAAATGATGCCTGGACTGAAAGATTATCCCCTTTATCCCTACCTGGAATACCGTCAGATCACCGATTACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGccag                                                                                           <  2:538987/250‑3 (MQ=255)
                                                                                                                                                        aGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGaa                                                                                       <  1:533939/251‑1 (MQ=255)
                                                                                                                                                              ccTGGGATAAGCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCGACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCGCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGccc                                                                                 <  2:193017/251‑1 (MQ=255)
                                                                                                                                                                        acgtcAAATGGATGTGGTCGAACAAATGATGCCTGGACTGGAAGATTATCCGCTTTATCCCTACCTGGAATACCGTCAGACCGCCGACGACCTGATGATTCAGCCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCGCGCTTCCGCCCGCTCGCACGCTGCAATCTCGGTTCGTGAATGCACCGGCGTTTCGTGAGCAGTGGCGTGTATTGTTAGCGTTTAGCCCGGAAAGGCCGGGatctgccg                                                                       <  2:318266/247‑8 (MQ=255)
                                                                                                                                                                         cGACAAATGGATGTGGACGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACCAACTTTTGTCCCGCTAACCCCAAGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGACAATGAACTGGGGCGGCGTGAGGGAGGGCGTGGCTTTTTAGGCCTTTTACCCGGaaaaaccgggacaccatg                                                                       >  2:178438/1‑237 (MQ=255)
                                                                                                                                                                           aCAAAGGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAAGCCCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAg                                                                    <  2:353216/251‑1 (MQ=255)
                                                                                                                                                                            cAAATTGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAgc                                                                   <  2:286030/251‑1 (MQ=255)
                                                                                                                                                                                       gggcGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAAtt                                                       <  1:294110/248‑1 (MQ=255)
                                                                                                                                                                                        gTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGACCTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAACAGACCGGAACTACCGAAGCGCACAGTAAtt                                                       >  1:51851/1‑250 (MQ=255)
                                                                                                                                                                                          cGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTTTTCTCCTTTATCACCGGAAAGACCCGGACAACCCAAGCGCCATGTAAAtact                                                    >  2:287449/1‑251 (MQ=255)
                                                                                                                                                                                          cGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTTTTAGCCTTTTGCCCGGAAAAACCCGGGACCACCAAAGCGCattgaacttact                                                     >  2:447467/1‑240 (MQ=255)
                                                                                                                                                                                                 atgatgCCTGGACTGAAGGATTATCCGCTTTATCCCTCCCTGGAATACCGCCAGATCACCTATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATg                                                                                                                                            <  2:240342/160‑1 (MQ=255)
                                                                                                                                                                                                 atgatgCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATg                                                                                                                                            >  1:240342/1‑160 (MQ=255)
                                                                                                                                                                                                   gatgCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCCATGTAATTACTACTATGca                                            >  1:358856/1‑249 (MQ=255)
                                                                                                                                                                                                    atgCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTTGCCCGGGAAAACCCCGGACTAACGGAGCGCCATGGAATTActactaagcaa                                           >  1:2578/1‑245 (MQ=255)
                                                                                                                                                                                                          ggACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCCAAATggg                                     >  1:603647/1‑249 (MQ=255)
                                                                                                                                                                                                             cTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTCACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGCACCGGCGCGGCGTGAAGACTGGCGTTGCTTGTTTACCTTTAGCCCGGGCAAGCCCGGGACCACCGCAGAGCAATGTACTTACTAATATTCAAAACGGAacac                                 >  2:651609/1‑251 (MQ=255)
                                                                                                                                                                                                                gaaagATTATCCACTTTATCCCTACCTGGAATCCCGTCAGATCACCGATGACCTTATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCTGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACAcc                                <  1:471997/246‑1 (MQ=255)
                                                                                                                                                                                                                aaGGATTATCCCCCTTATCCCTGCCTGGAATTCCGTCAGATCCCCGATGATCTGATGAAGCAACCGGCGGTGACTGCCACTAACTTTGTTCGCTCTAACCCTACGCTTCCTCCCGCTCGCACGCTGCAATCTCGCTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGTCTTTAGCCCGGAAAAGCCAGGAACTGCCGAAGCGCAATGTAATTACTACTATGCCAAATggtccaccgt                              <  2:1980/251‑9 (MQ=255)
                                                                                                                                                                                                                   agagtctCCGCTGTATCTCTACCTGGAATACCGACATATCACCGATGACCTGATGAATCATTCGGCGGTGACTGGCACTAACCTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCATTCTCGTTTCGTCAAGGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGTCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGc                            <  2:181659/245‑1 (MQ=255)
                                                                                                                                                                                                                    gATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTACCTTTAGACCGGAAAAGACCGGGACTAACGAAGAGCACTTGTATTAATACTATGCAAAAAGGAaccacgggc                            >  1:410801/1‑243 (MQ=255)
                                                                                                                                                                                                                         tccacTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGGAGTACAGAAGCGCAATTTAATTACTACTATGCCAAATTGAACAACCGGCAGGGTg                      >  1:488034/5‑251 (MQ=255)
                                                                                                                                                                                                                                        ccTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATATAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAAAAACCTGGCa         >  1:145881/1‑249 (MQ=255)
                                                                                                                                                                                                                                           ggAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAgggg    <  1:673734/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aTACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCg  >  1:165437/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
ATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCG  >  NC_000913/4630702‑4631188

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