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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,8370TC97.2% 206.2 / ‑4.5 71S35S (AGT→AGC) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (35/34);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCG  >  NC_000913/4630592‑4631079
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
cTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGtt                                                                                                                                                                                                                                                   >  1:295135/1‑250 (MQ=255)
  aaTCCCAGTGGGGCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTATCGTGTCGTAAACGGCAATTACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCAATGATGTATTTACACTTAGAGGATTCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTGGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCTCTGGATGATCAGCGTAGCCGTTat                                                                                                                                                                                                                                                 <  2:605671/250‑2 (MQ=255)
               cGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCa                                                                                                                                                                                                                                   <  1:171294/251‑1 (MQ=255)
                gATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGc                                                                                                                                                                                                                                    >  1:181659/1‑249 (MQ=255)
                gATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAg                                                                                                                                                                                                                                  <  1:360874/251‑1 (MQ=255)
                gATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAg                                                                                                                                                                                                                                  <  1:283220/251‑1 (MQ=255)
                          ccAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATATATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCTCAAATCAAGCAGGCCTGGGATa                                                                                                                                                                                                                        <  2:55247/250‑1 (MQ=255)
                                aaagtaaaaAAGTGTCCGTAACTGGGCGTAAACGGCAATGACTGGTTAGCATAAAGCTATTTCGCGGCATGACGCTGCATTGATGTATTTACGCTTAGAGGATGCGCTGGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCGGCCGGTGTCTGGCTGCGGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGTATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAc                                                                                                                                                                                                                  <  2:88572/251‑1 (MQ=255)
                                 aagtaaaaAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAc                                                                                                                                                                                                                  >  1:54719/1‑250 (MQ=255)
                                 aagtaaaaAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAc                                                                                                                                                                                                                  >  1:87369/1‑250 (MQ=255)
                                 aaataaaaAAGTGTCCGTAACGTGGCGTAAACGGCAATGACCGGTTAGCATAAATCAATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTCCCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACa                                                                                                                                                                                                                 <  2:152116/248‑1 (MQ=255)
                                     aaaaaaGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATg                                                                                                                                                                                                             <  2:54719/251‑1 (MQ=255)
                                       aaaaGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATgg                                                                                                                                                                                                            >  1:114568/1‑250 (MQ=255)
                                        aaaGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCGTGGGCTAATCGACAAATGGAt                                                                                                                                                                                                          >  1:457147/1‑251 (MQ=255)
                                        aaaGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGAt                                                                                                                                                                                                          >  1:75936/1‑251 (MQ=255)
                                           gtgtCCGTAAAGTGGCGTAAACGGCAATGACTGGATAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGAAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCAGGGGTCTGTCTGCTGACGGTCAGCAGCGGGGCGAGAGCCGACTCACTGGATAAGCAGAGTAGCCGTTAAGCCCAAATCAAGCAGGCCCGGGGTAATCGACaacaggagatg                                                                                                                                                                                                       >  2:600263/1‑242 (MQ=255)
                                                   aaCGTGTCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACTCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGGGTCTGTCTGCTGCCGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCACAAATCAAGCAGGCCTGGGATAATCGACAAATGGATTTGGTCgcac                                                                                                                                                                                                >  2:126402/1‑247 (MQ=255)
                                                      ggggCGTAAACGGCAATGTCTGGTTATCATAAATCTATTCCGCGGCATGACGCTGAATGGATGTATTTCCACGTAGAGGATGCGCTTTTGGAAAAGGCCAAACAAGTTCCCTGGCGGCTGTTGGCTGCCGGGGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGTCTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTTGTCGAACAAat                                                                                                                                                                                            <  2:75936/249‑1 (MQ=255)
                                                      gTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAat                                                                                                                                                                                            >  1:538987/1‑251 (MQ=255)
                                                          ccgatcgggCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCGGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAatgatg                                                                                                                                                                                        <  2:123436/244‑1 (MQ=255)
                                                                 ggCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAatgatg                                                                                                                                                                                        >  1:123436/1‑244 (MQ=255)
                                                                             ttAGTATAAATCTAATCCGCGGCATGAAGCTGCATGGATGTATTTTCACTTAGAGGATGCGCTGGTGGAAAAAGCCAAAAAAGTTACCGGTCGGCTGTTGGCTGCAGGTGTCTGGCTGCGGACGGACAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGAGGTCGACCAATTGATGCCTGGACTGA‑GAGATTAt                                                                                                                                                                     <  2:426374/250‑1 (MQ=255)
                                                                               aGCATAAAGCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATcc                                                                                                                                                                   <  2:114568/250‑1 (MQ=255)
                                                                                     aaTCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCATCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGTAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTAt                                                                                                                                                            <  1:405626/251‑1 (MQ=255)
                                                                                     aaTCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTAt                                                                                                                                                            <  1:270830/251‑1 (MQ=255)
                                                                                                cgGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTgg                                                                                                                                                  >  1:286030/1‑250 (MQ=255)
                                                                                                cgGCATGACGCTGCATGGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTgg                                                                                                                                                  >  2:489045/1‑250 (MQ=255)
                                                                                                    aTGACGATTCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTCCCGGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGAAGCGTGGCGCGAGCCGACTCACGTGATGAGCAGAGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTTGAATAc                                                                                                                                             <  2:295135/251‑1 (MQ=255)
                                                                                                    aTGACACTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATAc                                                                                                                                             >  1:102375/1‑251 (MQ=255)
                                                                                                         gCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCGGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCgtc                                                                                                                                         >  2:443889/1‑248 (MQ=255)
                                                                                                              aTTGATTTTTTTCCACTTTGAGTATGCGCTGGTGGAAAAAGCTAAAAAAGTTACCTGGCGGCTGGTGGCTGCCGGTGTCTGTCTGTTTACGTGCAGCATCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATAAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTAGCCCTACCTTGAAGACCGTCAGATCa                                                                                                                                   <  2:477726/251‑1 (MQ=255)
                                                                                                              aTTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTATGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAgctc                                                                                                                                    >  1:627631/1‑247 (MQ=255)
                                                                                                              aTTGATGTATTTACACTTAGAGGAGGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTTCCTGGAATACCGTCAGATCa                                                                                                                                   <  1:464466/251‑1 (MQ=255)
                                                                                                                 gATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACACACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCg                                                                                                                                >  2:71199/1‑251 (MQ=255)
                                                                                                                               tAGAGGATGCGCTTGGGGAAAAAGCCAAACAAGTTACCTGGCGGGTGTTGGCTGTCGGTGTCTGTCTGCTGACGGTCAGGAGCGTGGCGCGAGCCGACTCCCTGGATGAGCAGCGTAGCCGTTACGCGAAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATCCCGTCAGATCACCGATGACCtgatga                                                                                                                    <  1:126402/249‑1 (MQ=255)
                                                                                                                                agagGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATc                                                                                                                 >  1:559521/1‑251 (MQ=255)
                                                                                                                                     aTGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACcg                                                                                                            <  1:496349/251‑1 (MQ=255)
                                                                                                                                            tgtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGACTGCTGACGGTCAGCAGCGTGGCGCGAGCCGCCTCACTGGATGAGCAGCGTAGCAGTTACGCGCAAATCATGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAAAGATGCCTGGACTGA‑AAGATTATCCACATTATCCCTACCTGGCATACCGTCAGATCACAGATGACCTGGTGAAACAAAcggcag                                                                                                        >  2:366454/1‑246 (MQ=255)
                                                                                                                                             gtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGACGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGGTCACCGATGACCTGATGAATCAACCGGCGGTg                                                                                                      >  2:82426/1‑249 (MQ=255)
                                                                                                                                             ggggaaaaccccaacaaaGTTACCTGGCGGGTGTGGGCGTCCGGTGTCGGTCTTCTGACGGTCAGCAGCTTGGAGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCACATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AATATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑c                                                                                                   <  2:589543/235‑1 (MQ=255)
                                                                                                                                               ggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGTCTGCCGGTGCCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑ctg                                                                                                 <  1:352323/251‑3 (MQ=255)
                                                                                                                                               ggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑ctg                                                                                                 <  1:461966/251‑3 (MQ=255)
                                                                                                                                                            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)
                                                                                                                                                                  aCCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGCCGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTg                                                                                                                                                   >  2:180581/1‑183 (MQ=255)
                                                                                                                                                                  aCCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTg                                                                                                                                                   <  1:180581/183‑1 (MQ=255)
                                                                                                                                                                     tGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CCCTAACTTTGTTCGCGCTa                                                                           >  1:353216/1‑250 (MQ=255)
                                                                                                                                                                       gcggcTGTTGGCTGCCGGTTTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgatga                                                                                                                    <  2:364474/209‑1 (MQ=255)
                                                                                                                                                                       gcggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgatga                                                                                                                    >  1:364474/1‑209 (MQ=255)
                                                                                                                                                                          gcTGTTGGCTGCCGGTGTGTGTCGGCAGACGGTCAGCAGCGTGGCGCGGGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCTACAAATGGATGTGGGCGAACAAATGATGCCTggtcta                                                                                                                                                                               <  2:431862/148‑5 (MQ=255)
                                                                                                                                                                          gcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTg                                                                                                                                                                               >  1:431862/1‑148 (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)
                                                                                                                                                                                                                        cgAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTa                                                                                                                                                             >  1:106730/1‑119 (MQ=255)
                                                                                                                                                                                                                        cgAGCCGACTCACGGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGTGATAATCGACAAATGGATGTGGTCGAACAAAGGATGCCTGGACTGA‑AAGATTATCCACTTTa                                                                                                                                                             <  2:106730/119‑1 (MQ=255)
                                                                                                                                                                                                                          aGCCGACTCACTGGATGAGCCGCGTAGCCGTTACGCGCAAATCAAGCAGTCCTTGGATAAACGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTTA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCGCGTTTCGTCAATGAACTgg                     <  2:87369/251‑1 (MQ=255)
                                                                                                                                                                                                                             cGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGTACTGA‑AAGATTAGCCACTTTa                                                                                                                                                             <  2:671252/114‑1 (MQ=255)
                                                                                                                                                                                                                             cGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTa                                                                                                                                                             >  1:671252/1‑114 (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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCG  >  NC_000913/4630592‑4631079

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