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,9840GT91.4% 155.1 / ‑2.3 59T84T (ACG→ACT) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base T (29/24);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCT  >  NC_000913/4630748‑4631228
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                           
cAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGA‑CTGT‑CCCTAAATTTGt                                                                                                                                                                                                                                              >  2:20449/1‑251 (MQ=255)
 aaGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAACTCAAGCAGGCCTGGGATAATCGACAAATGTATGTGGTCGAACAAATGCTACCTGGACTGC‑AAGATTAACCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGTTGAATCAACCGGCGGTAAATGGT‑‑ACTAACTTTGtt                                                                                                                                                                                                                                             >  2:353660/1‑251 (MQ=255)
  aGTTACCGGGCGGCTGTTGGCTGCCGGTGTCTGTTGGCTGACGGTCAGGAGCGTGTCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCTGCGGTTA‑CTGT‑CACTAAGTTTGTTc                                                                                                                                                                                                                                            <  2:598542/251‑1 (MQ=255)
         tGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CCCTAACTTTGTTCGCGCTa                                                                                                                                                                                                                                      >  1:353216/1‑250 (MQ=255)
              gcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCa                                                                                                                                                                                                                                <  1:20449/251‑1 (MQ=255)
                tGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGA‑ACGT‑CACTAACTTTGTTGGCGCTTACCCCAcc                                                                                                                                                                                                                              >  2:373474/1‑250 (MQ=255)
                   tGGCTGCGGGTGTCTGTATGCTGACGGTCCGCAGCGTGGCGCGAGACGCCTACCTGGATGAGCAGCGTAGCCGTTAAGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACACATGATGCCTGGACTGACAGG‑TTATCAAATTTATCCCTACCTGGAATACCGTCAGATCACCGATGCCCTGCTGAATACACCGGCGGTGT‑CTGTGCA‑TAACTTTGTTCGCTCTAAACCCACGcat                                                                                                                                                                                                                           >  2:335787/1‑249 (MQ=255)
                        gccgggggcggcctgctgACGGTCAGCAGCGTGGCGCGAGCCGACTCACTTGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGTGGA‑CTGT‑CTCTAACTTTGCTCGCGCTAACCCCACGCTtccg                                                                                                                                                                                                                        <  2:248744/237‑2 (MQ=255)
                              gtctgtctGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACTCTTCCTCCCGCTCg                                                                                                                                                                                                                <  1:560023/251‑1 (MQ=255)
                               tctgtctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGc                                                                                                                                                                                                               >  1:618920/1‑251 (MQ=255)
                                       ctgACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGAGGTGA‑CTGT‑CACTAACTTTGTTATCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGcc                                                                                                                                                                                                       >  2:627777/1‑250 (MQ=255)
                                               cagcagCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGCAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGCTGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGtt                                                                                                                                                                                               >  1:279981/1‑251 (MQ=255)
                                                   agcGTGGCGCGAACCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACATGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTACCGCTCGCACGCTGCCATCTCGTTTcat                                                                                                                                                                                           >  1:318062/1‑249 (MQ=255)
                                                          cgcgAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGCCTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGG‑CACTAACTTTGTTCACGCTAACCCCACGATTCCTCCCGCTCGAAAGCTGCAATCTCGTTTCGTCAATGa                                                                                                                                                                                     >  1:1980/1‑250 (MQ=255)
                                                              aGCCGACTCACTGGATGAGCCGCGTAGCCGTTACGCGCAAATCAAGCAGTCCTTGGATAAACGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTTA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCGCGTTTCGTCAATGAACTgg                                                                                                                                                                                <  2:87369/251‑1 (MQ=255)
                                                                 cGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTAACCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGgcgc                                                                                                                                                                             <  1:353660/251‑1 (MQ=255)
                                                                       aCTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTg                                                                                                                                                                       <  1:510282/251‑1 (MQ=255)
                                                                              gAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCATGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTg                                                                                                                                                                >  1:223674/1‑251 (MQ=255)
                                                                                 cagcGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCg                                                                                                                                                             <  1:373474/251‑1 (MQ=255)
                                                                                        gccGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTg                                                                                                                                                      >  1:412412/3‑251 (MQ=255)
                                                                                         ccGTTACGCTCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTTGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCTCCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGt                                                                                                                                                     <  2:321099/250‑1 (MQ=255)
                                                                                            ttACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGGGGCtttttct                                                                                                                                                  >  2:533931/1‑246 (MQ=255)
                                                                                                     aaTCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGc‑c                                                                                                                                       <  1:660107/251‑1 (MQ=255)
                                                                                                       ccAAGCAGGCCTGGGATAATCGACAAATGGATGTGGGCGAACAAATGATGCCTGGACTGA‑AAGATTATCCCCTTTATCCCTACCTGGAATACCGTCAGATCACCGATTACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGc‑cag                                                                                                                                     <  2:538987/250‑3 (MQ=255)
                                                                                                          aGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑aa                                                                                                                                 <  1:533939/251‑1 (MQ=255)
                                                                                                                ccTGGGATAAGCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATCCGCTTTATCCCGACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCGCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGccc                                                                                                                           <  2:193017/251‑1 (MQ=255)
                                                                                                                          acgtcAAATGGATGTGGTCGAACAAATGATGCCTGGACTGG‑AAGATTATCCGCTTTATCCCTACCTGGAATACCGTCAGACCGCCGACGACCTGATGATTCAGCCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCGCGCTTCCGCCCGCTCGCACGCTGCAATCTCGGTTCGTGAATGCACCGGCGTTTCGTGAGCAGTGGCGTGTATTGTTAG‑CGTTTAGC‑CCGG‑AAAGGCCGGGatctgccg                                                                                                                 <  2:318266/247‑8 (MQ=255)
                                                                                                                           cGACAAATGGATGTGGACGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACCAACTTTTGTCCCGCTAACCCCAAGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGACAATGAACTGGGGCGGCGTGAGGGAGGGCGTGGCTTTTTAGGCCTTTTAC‑CCGG‑aaaaaccgggacaccatg                                                                                                                 >  2:178438/1‑237 (MQ=255)
                                                                                                                             aCAAAGGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAAGCCCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAg                                                                                                              <  2:353216/251‑1 (MQ=255)
                                                                                                                              cAAATTGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAgc                                                                                                             <  2:286030/251‑1 (MQ=255)
                                                                                                                                         gggcGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAAtt                                                                                                <  1:294110/248‑1 (MQ=255)
                                                                                                                                          gTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGACCTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AACAGACCGGAACTACCGAAGCGCA‑CA‑GTAAtt                                                                                                >  1:51851/1‑250 (MQ=255)
                                                                                                                                            cGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTTTTCT‑CCTTTATCACCGG‑AAAGACCCGGA‑CAACCCAAGCGCC‑AT‑GTAAAtact                                                                                             >  2:287449/1‑251 (MQ=255)
                                                                                                                                            cGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTTTTAG‑CCTTTTGC‑CCGG‑AAAAACCCGGGACCACCAAAGCGCattgaacttact                                                                                               >  2:447467/1‑240 (MQ=255)
                                                                                                                                                   atgatgCCTGGACTGA‑AGGATTATCCGCTTTATCCCTCCCTGGAATACCGCCAGATCACCTATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATg                                                                                                                                                                                      <  2:240342/160‑1 (MQ=255)
                                                                                                                                                   atgatgCCTGGACTGA‑AGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATg                                                                                                                                                                                      >  1:240342/1‑160 (MQ=255)
                                                                                                                                                     gatgCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCC‑AT‑GTAATTACTACTATGca                                                                                     >  1:358856/1‑249 (MQ=255)
                                                                                                                                                      atgCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTG‑CCTTTTGC‑CCGGGAAAACCCCGGA‑CTAACGGAGCGCC‑AT‑GGAATTActactaagcaa                                                                                    >  1:2578/1‑245 (MQ=255)
                                                                                                                                                            ggACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCCAAATggg                                                                              >  1:603647/1‑249 (MQ=255)
                                                                                                                                                               cTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTCACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGCACCGGCGCGGCGTGAAGACTGGCGTTGCTTGTTTA‑CCTTTAGC‑CCGG‑GCAAGCCCGGGACCACCGCAGAGCA‑AT‑GTACTTACTAATATTCAAAACGGAacac                                                                          >  2:651609/1‑251 (MQ=255)
                                                                                                                                                                  gaaagATTATCCACTTTATCCCTACCTGGAATCCCGTCAGATCACCGATGACCTTATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CTGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACAcc                                                                         <  1:471997/246‑1 (MQ=255)
                                                                                                                                                                  a‑aGGATTATCCCCCTTATCCCTGCCTGGAATTCCGTCAGATCCCCGATGATCTGATGAAGCAACCGGCGGTGA‑CTGC‑CACTAACTTTGTTCGCTCTAACCCTACGCTTCCTCCCGCTCGCACGCTGCAATCTCGCTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTG‑TCTTTAGC‑CCGG‑AAAAGCCAGGAACTGCCGAAGCGCA‑AT‑GTAATTACTACTATGCCAAATggtccaccgt                                                                       <  2:1980/251‑9 (MQ=255)
                                                                                                                                                                     agagtctCCGCTGTATCTCTACCTGGAATACCGACATATCACCGATGACCTGATGAATCATTCGGCGGTGA‑CTGG‑CACTAACCTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCATTCTCGTTTCGTCAAGGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGTCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGc                                                                     <  2:181659/245‑1 (MQ=255)
                                                                                                                                                                      gATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTA‑CCTTTAGA‑CCGG‑AAAAGACCGGGACTAACGAAGAGCACTT‑GTA‑TTAATACTATGCAAAAAGGAaccacgggc                                                                     >  1:410801/1‑243 (MQ=255)
                                                                                                                                                                           tccacTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGGAGTACAGAAGCGCA‑AT‑TTAATTACTACTATGCCAAATTGAACAACCGGCAGGGTg                                                               >  1:488034/5‑251 (MQ=255)
                                                                                                                                                                                          ccTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑ATAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAAAAACCTGGCa                                                  >  1:145881/1‑249 (MQ=255)
                                                                                                                                                                                             ggAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑gggg                                            <  1:673734/251‑1 (MQ=255)
                                                                                                                                                                                                aTACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCg                                          >  1:165437/1‑250 (MQ=255)
                                                                                                                                                                                                          atcaacgaagaccTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTGTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCtgtg                               <  2:488307/239‑1 (MQ=255)
                                                                                                                                                                                                           tCACCGATGACCTGATGAATCAACCGGCGGTTA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGAGTGGCTTGTTTG‑TATTTAGC‑CCGG‑GAAAACCCGGGAATACCGAAACGCA‑ATTGCAAT‑AATAATATGCAAAAAGGAAAACCGGGCAGAGGGAAGAAACCCGGCAAAGGGGGGAAaggctgagg                               >  2:645864/1‑244 (MQ=255)
                                                                                                                                                                                                                 atgacCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAAcc                        <  1:82426/251‑1 (MQ=255)
                                                                                                                                                                                                                   gaccTGATGAATCAACCGGCTGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCGTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGGGCTGTGGCTAACCg                       <  2:627631/247‑1 (MQ=255)
                                                                                                                                                                                                                    accTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GGAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAAACggg                     >  1:566074/3‑250 (MQ=255)
                                                                                                                                                                                                                            aaTCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAA‑CCTTTAGC‑CCGG‑AAAAACCCGGGACTACCGAAGCGCA‑AT‑GTAATTACTACTATTCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCGGTGGCTAAACGGCAAGagccc             >  1:368735/1‑250 (MQ=255)
                                                                                                                                                                                                                               cAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGCAAGagccagcc          <  1:532176/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cggcggTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑GAAAACCCGGAACTACCGCAGCGCA‑AT‑GTAATTACTTATCTGCAAAATGGAACACCGGGCAGAGTGGAGAAGCCTGGGAA‑GGGGGGAAAGAGCTGTGGCTAACCGGCAAAAGGCAGCCTaa       >  1:391693/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cggcggTGA‑CGGT‑CCCTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTTAAGACTGGCGTGGCTTGTTAG‑CCGTTAGC‑CCGG‑AAAAGCCCGGACCTACCGAAGCGCA‑AT‑GTAATTACTACTATGCATAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAAc      <  1:267581/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCGGCTCGCACGCTGCAATCTCGATTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTTGC‑CCGG‑AGAAACCCGGAAAAAACGAAGCGCA‑AT‑GTGATTACTACCATGCAACAAGGAAAAACAGGCAGAGGGAAGACAACTTGCAA‑GGGGCGAAAGAGCTAGGGCCCACCGGCCAGagccaggccagagcac  >  2:122770/1‑241 (MQ=255)
                                                                                                                                                                                                                                       ggTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCCAAATGGAACACCGGGCAGAGTGGAGAAGCCTGGCAA‑GGGGCGGAAGAGCTGTGGCTTACCGGCAAGAGCCAGCCTAACGccc  >  1:237614/1‑250 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                           
CAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCT  >  NC_000913/4630748‑4631228

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