Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 25,598 A→G D131G (GAC→GGC)  lspA → lipoprotein signal peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091325,5980AG97.6% 116.4 / ‑4.2 42D131G (GAC→GGC) lspAlipoprotein signal peptidase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (20/21);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

CATTATGCGCGTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGT‑ATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTAC  >  NC_000913/25351‑25841
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
aaTTAGGCGCGTAACTATGGCGCGGCTGTTAGTTCCCTTGCCGATAGCGGGGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGT‑ATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATGGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTg                                                                                                                                                                                                                                                  <  1:444583/250‑1 (MQ=255)
        gcgTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGT‑ATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCACTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGCCTTCTACGTCGGCGGCTGGCACTTcc                                                                                                                                                                                                                                          >  1:163732/1‑250 (MQ=255)
         cgTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGGTGGTTCTTTGCCGGTATTGCGATTGGT‑ATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTAGTTTTCGATATGATCGACTTCTCCGTCGGCTGCCGGCACTTCGc                                                                                                                                                                                                                                         >  2:234391/1‑251 (MQ=255)
             aCTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGG‑ATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACTCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGAGCGCTGGGCAACCTGTTCGGCCTCCTGTGGCACGGCTTCGTTTGCTGTTTTATCGACTTCTACGTTGGCGGCTGGCACTTCGCCAcc                                                                                                                                                                                                                                     >  2:396011/1‑251 (MQ=255)
                atggcgctgtttttAGTTTCCGTGCCGATCGCGGCGGCTGGCAGGGTTGGTTCTTTGCCGGTATTGCGTTTGTT‑ATTAGCGTGTACCTGGCAGTGATGATGTATCGCTCAAGGGCCACTCAGAAGCTAAACATTA‑TCGCTTACGCGCTGTTT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCAcc                                                                                                                                                                                                                                     <  1:354754/237‑1 (MQ=255)
                            ttAGTTTCCTGGCCGATAGCGGCGGCTGGCAGCGTTGGTTTTTTGCCGGTATTGCGTTTGGT‑ATTAGCGTGACCCTGGCAGTGAGGATGTACCCCTCGAAGCCCACGCAGAAGCAAAACAATA‑TCGTTTACGCGCTGATT‑ATTGGCGGCGCCCGGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCTGCTGGCACTTCGCCACCTTCAACCTTGCCGa                                                                                                                                                                                                                       <  1:115823/250‑1 (MQ=255)
                                          gATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGT‑ATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGTTACTGCCACCtgttt                                                                                                                                                                                                        >  1:295408/1‑249 (MQ=255)
                                               gggcggcTGGCAGCGTTGGTTCTTGGCCGGTATTGCG‑TTGGTTATTAGCGTGTCCCTGGCAGTGATTATGTACCGCTCGAAGGCCACGCAGAAGCTAAACAATT‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGTCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCCCCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTg                                                                                                                                                                                                   <  1:450990/250‑1 (MQ=255)
                                                 gcggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGT‑ATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCg                                                                                                                                                                                                 >  1:118285/1‑251 (MQ=255)
                                                       ggcggCGTTGGTTTTTTGCTGGTATTGCGATTGGT‑ATTAGCGGGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCAAAACAATA‑TCGCTTCCGCGCTGATT‑ATGGGCGGCGCGCGGGGCACCCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTg                                                                                                                                                                                           <  1:169452/247‑1 (MQ=255)
                                                         gaGCGTTGTTTCTTTGCCTTTATTGCTATTGGT‑ATTAGCGTGATCCTGGAAGTGATGATGTATCGCTCGAAGTCCACGCAGAAGCTAAACATTT‑TCGCTTACGCGCGGATT‑ATTGGCGGCGCGCGGGGCAACCTGTTCGTCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGAt                                                                                                                                                                                         <  1:41218/250‑1 (MQ=255)
                                                               ggtttCTTTGTGGGTATTGCGATTGGT‑ATTAGCGGGATCCTGGCAGTGTGGATGTTTCGTACGAAGGCCACGCAGAAGCAAAACAATA‑TCGCTTGCGCTCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCt                                                                                                                                                                                   <  1:296191/248‑1 (MQ=255)
                                                                 gTTCTTTGCCGGTATTGCGATTGGT‑ATTAGCGTGATCCTGTCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTACACCTTGCCGCTACTGCCATCTGTGTCGGTTCTGCACCGATTGGGCTgg                                                                                                                                                                                 >  2:140210/1‑251 (MQ=255)
                                                                   tgttggCCGGTATGGCGATTGGT‑ATGAGGTTTACCCTGCCAGTGTTGAGTTATCGCTCGAGGGCCCCGCAGAAGTAAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGTCGCGCGGGGCAACCTTTTCGTCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGa                                                                                                                                                                                <  1:40085/245‑1 (MQ=255)
                                                                          cGGTATTGCGATTGGT‑ATTCGCGTGATCCTGGCAGTGATGCTGTATCGATCGAGGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCTGTATGATCGTCTTCTACGTCGGCGGCTGGCTCTTCGCCACCTTCAACCTTGCCGATATAGCCATCTTTTTCGGTATGGCCCTGATTGTGCTGGACGGgtttt                                                                                                                                                                        >  2:193490/1‑251 (MQ=255)
                                                                                     ttGGT‑ATTAGCGTGATCCTGGCAGTGCTGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGAGCTAGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCAGGCACTTCGCCACCTTCAACCTTTCCGATACTGCCATCGTTGTCGGTGCGGCAATGATTGTGCTTGTAGGTTTTTTACCttttaa                                                                                                                                                              >  2:90449/1‑246 (MQ=255)
                                                                                         t‑ATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGGTATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGGTTTTTACCTTCTAAAGCgg                                                                                                                                                         >  2:235443/1‑250 (MQ=255)
                                                                                                      ccTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAGGCGAAAAAACaataat                                                                                                                                             >  1:277015/1‑251 (MQ=255)
                                                                                                        tGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAataata                                                                                                                                            >  1:41282/1‑250 (MQ=255)
                                                                                                          gCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAAc                                                                                                                                         >  1:464882/1‑251 (MQ=255)
                                                                                                             gggttgtgggatctcgcagaGGCCACGCAGAAGC‑AAACAAAAGTCGCTT‑CGGGGTGATTTATGGGCGGCGCGCGGGACAACCTGTCCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACGTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTTCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCt                                                                                                                                      <  2:464882/232‑1 (MQ=255)
                                                                                                             gTGAGGATGTTTCGCTCGAAGGCCCGGAAGAAGTTAAACAATA‑GCGTTTACGCGCTAATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCt                                                                                                                                      <  1:396011/251‑1 (MQ=255)
                                                                                                                     gtaccGCTCGAGGGCACCGCAGAAGCAAAACAATA‑TCTCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGgatg                                                                                                                              <  1:155693/247‑1 (MQ=255)
                                                                                                                             cGAAGGCCACGCAGAAGCTAAACAATA‑TAGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACATTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACATTAATAAACCCTGCCggctgagctgctg                                                                                                                      >  1:313479/1‑240 (MQ=255)
                                                                                                                                          gAAGCGAAAAAAAA‑TCGCTGACGCGGTGATT‑TTCGGCGGCGGGTTGGGCAGCCTGTAGGACCGCCTGTGGCACGGCTTCGTTGTCGATATGAGCGACTTCTACGTCGGGGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATcc                                                                                                         <  1:367503/251‑1 (MQ=255)
                                                                                                                                              cTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGcc                                                                                                     >  1:356349/1‑251 (MQ=255)
                                                                                                                                              cTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGAACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGGTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGCAGCATCTTATCCGGcc                                                                                                     >  1:70507/1‑251 (MQ=255)
                                                                                                                                                      atagcgtttACGCGTTGATT‑ATGGGCTGGGCGCGTGGCAACCTTTTCAGCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGAt                                                                                              <  1:94457/244‑1 (MQ=255)
                                                                                                                                                                                           gCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGCAATCGTAGGCCGGATAAGGCGTTTAcgc                                                          >  1:346049/1‑251 (MQ=255)
                                                                                                                                                                                             atcgtgctcgcccGCCTGTGGCTCGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTAcgccg                                                        <  2:363819/240‑1 (MQ=255)
                                                                                                                                                                                                cTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGCTGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                     >  2:500353/1‑251 (MQ=255)
                                                                                                                                                                                                  gTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                   <  1:205020/251‑1 (MQ=255)
                                                                                                                                                                                                           cctggtgccccggttCGTTGGCGATTTGATTGTCTTCTGCGTCGGTGGCTGGCTATGCGCCACCTTCAACCTTGCCGTTACTGCCATCTGTGTCGGTGCGGGAATGATTGTGCTGGAAGGTTTTTTTCCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAAGCGTAGTCCGTATAAGGCGTTTACGCCGCATCCGGCAAAAAt                                          <  2:490973/238‑1 (MQ=255)
                                                                                                                                                                                                            cTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCAGGATAAGGCGTTTACACCGCATCCGGCAAAAATc                                         >  1:237425/1‑251 (MQ=255)
                                                                                                                                                                                                              gtgGCACGGCTTCGTTGTCGATATGAGCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATc                                         <  1:406612/249‑1 (MQ=255)
                                                                                                                                                                                                                          cGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAagag                           <  2:437068/251‑1 (MQ=255)
                                                                                                                                                                                                                              gTCGATATTATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCaaa                       <  1:378299/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gttaggagcGACTTCTACGGCGGCGGCTGGCACTTCGCCACCTTCAACCTGGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTAGCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAACCCt                    <  2:99220/243‑1 (MQ=255)
                                                                                                                                                                                                                                   tatGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGCAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCATTAAATATAAGAGAAAACCTGc                  >  2:429861/1‑251 (MQ=255)
                                                                                                                                                                                                                                    atGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGc                  <  1:234391/250‑1 (MQ=255)
                                                                                                                                                                                                                                    atGATCGAATTCTAAGACGGCGGCTGGAACTTCGCCACCTTCAACCTTGCCGGTACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAACAACAATAATAAACCCTGCCGGAAGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCAGTCACACATCCTTAAATATAAGAGCAAACCTGcc                 >  2:336976/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aggatcgctttCTACGTCGGCGGCTTGCACTTCGCCCCCTTCACCCTTGCCGATACTGCCATCTGTGTCGGTGCTGCACTGATTGTGCTGGAAGGTGTTTTCCCTTCTAAAGCGAAAAAACAATAATAAACCCGGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGc                  <  2:277015/241‑1 (MQ=255)
                                                                                                                                                                                                                                                     cggcggCTGGCACTTCGCCACCTTCAACCTTGCCGATAATGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAGACAATAATAAACCCTGCCGGCTTCGATGCTGACGCATCTTATCAGGCCTACCGATTGCTGCGAACTCGTAGGCCAGATAAGGCGTTTACGCCGCCTCCTGCACACAGCCTTCAATATAACCGCCAACCTGCATGTCTGAATCAGTAc  >  2:178241/1‑249 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
CATTATGCGCGTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGT‑ATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTAC  >  NC_000913/25351‑25841

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