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

CATTATGCGCGTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTAC  >  NC_000913/25351‑25841
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
aaTTAGGCGCGTAACTATGGCGCGGCTGTTAGTTCCCTTGCCGATAGCGGGGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATGGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTg                                                                                                                                                                                                                                                  <  1:444583/250‑1 (MQ=255)
        gcgTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCACTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGCCTTCTACGTCGGCGGCTGGCACTTcc                                                                                                                                                                                                                                          >  1:163732/1‑250 (MQ=255)
         cgTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGGTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTAGTTTTCGATATGATCGACTTCTCCGTCGGCTGCCGGCACTTCGc                                                                                                                                                                                                                                         >  2:234391/1‑251 (MQ=255)
             aCTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGGATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACTCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGAGCGCTGGGCAACCTGTTCGGCCTCCTGTGGCACGGCTTCGTTTGCTGTTTTATCGACTTCTACGTTGGCGGCTGGCACTTCGCCAcc                                                                                                                                                                                                                                     >  2:396011/1‑251 (MQ=255)
                atggcgctgtttttAGTTTCCGTGCCGATCGCGGCGGCTGGCAGGGTTGGTTCTTTGCCGGTATTGCGTTTGTTATTAGCGTGTACCTGGCAGTGATGATGTATCGCTCAAGGGCCACTCAGAAGCTAAACATTATCGCTTACGCGCTGTTTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCAcc                                                                                                                                                                                                                                     <  1:354754/237‑1 (MQ=255)
                            ttAGTTTCCTGGCCGATAGCGGCGGCTGGCAGCGTTGGTTTTTTGCCGGTATTGCGTTTGGTATTAGCGTGACCCTGGCAGTGAGGATGTACCCCTCGAAGCCCACGCAGAAGCAAAACAATATCGTTTACGCGCTGATTATTGGCGGCGCCCGGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCTGCTGGCACTTCGCCACCTTCAACCTTGCCGa                                                                                                                                                                                                                       <  1:115823/250‑1 (MQ=255)
                                          gATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGTTACTGCCACCtgttt                                                                                                                                                                                                        >  1:295408/1‑249 (MQ=255)
                                               gggcggcTGGCAGCGTTGGTTCTTGGCCGGTATTGCGTTGGTTATTAGCGTGTCCCTGGCAGTGATTATGTACCGCTCGAAGGCCACGCAGAAGCTAAACAATTTCGCTTACGCGCTGATTATTGGCGGCGCGCTGGTCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCCCCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTg                                                                                                                                                                                                   <  1:450990/250‑1 (MQ=255)
                                                 gcggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCg                                                                                                                                                                                                 >  1:118285/1‑251 (MQ=255)
                                                       ggcggCGTTGGTTTTTTGCTGGTATTGCGATTGGTATTAGCGGGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCAAAACAATATCGCTTCCGCGCTGATTATGGGCGGCGCGCGGGGCACCCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTg                                                                                                                                                                                           <  1:169452/247‑1 (MQ=255)
                                                         gaGCGTTGTTTCTTTGCCTTTATTGCTATTGGTATTAGCGTGATCCTGGAAGTGATGATGTATCGCTCGAAGTCCACGCAGAAGCTAAACATTTTCGCTTACGCGCGGATTATTGGCGGCGCGCGGGGCAACCTGTTCGTCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGAt                                                                                                                                                                                         <  1:41218/250‑1 (MQ=255)
                                                               ggtttCTTTGTGGGTATTGCGATTGGTATTAGCGGGATCCTGGCAGTGTGGATGTTTCGTACGAAGGCCACGCAGAAGCAAAACAATATCGCTTGCGCTCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCt                                                                                                                                                                                   <  1:296191/248‑1 (MQ=255)
                                                                 gTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGTCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTACACCTTGCCGCTACTGCCATCTGTGTCGGTTCTGCACCGATTGGGCTgg                                                                                                                                                                                 >  2:140210/1‑251 (MQ=255)
                                                                   tgttggCCGGTATGGCGATTGGTATGAGGTTTACCCTGCCAGTGTTGAGTTATCGCTCGAGGGCCCCGCAGAAGTAAAACAATATCGCTTACGCGCTGATTATTGGCGTCGCGCGGGGCAACCTTTTCGTCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGa                                                                                                                                                                                <  1:40085/245‑1 (MQ=255)
                                                                          cGGTATTGCGATTGGTATTCGCGTGATCCTGGCAGTGATGCTGTATCGATCGAGGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCTGTATGATCGTCTTCTACGTCGGCGGCTGGCTCTTCGCCACCTTCAACCTTGCCGATATAGCCATCTTTTTCGGTATGGCCCTGATTGTGCTGGACGGgtttt                                                                                                                                                                        >  2:193490/1‑251 (MQ=255)
                                                                                     ttGGTATTAGCGTGATCCTGGCAGTGCTGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGAGCTAGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCAGGCACTTCGCCACCTTCAACCTTTCCGATACTGCCATCGTTGTCGGTGCGGCAATGATTGTGCTTGTAGGTTTTTTACCttttaa                                                                                                                                                              >  2:90449/1‑246 (MQ=255)
                                                                                         tATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGGTATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGGTTTTTACCTTCTAAAGCgg                                                                                                                                                         >  2:235443/1‑250 (MQ=255)
                                                                                                      ccTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAGGCGAAAAAACaataat                                                                                                                                             >  1:277015/1‑251 (MQ=255)
                                                                                                        tGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAataata                                                                                                                                            >  1:41282/1‑250 (MQ=255)
                                                                                                          gCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAAc                                                                                                                                         >  1:464882/1‑251 (MQ=255)
                                                                                                             gggttgtgggatctcgcagaGGCCACGCAGAAGCAAACAAAAGTCGCTTCGGGGTGATTTATGGGCGGCGCGCGGGACAACCTGTCCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACGTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTTCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCt                                                                                                                                      <  2:464882/232‑1 (MQ=255)
                                                                                                             gTGAGGATGTTTCGCTCGAAGGCCCGGAAGAAGTTAAACAATAGCGTTTACGCGCTAATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCt                                                                                                                                      <  1:396011/251‑1 (MQ=255)
                                                                                                                     gtaccGCTCGAGGGCACCGCAGAAGCAAAACAATATCTCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGgatg                                                                                                                              <  1:155693/247‑1 (MQ=255)
                                                                                                                             cGAAGGCCACGCAGAAGCTAAACAATATAGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACATTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACATTAATAAACCCTGCCggctgagctgctg                                                                                                                      >  1:313479/1‑240 (MQ=255)
                                                                                                                                          gAAGCGAAAAAAAATCGCTGACGCGGTGATTTTCGGCGGCGGGTTGGGCAGCCTGTAGGACCGCCTGTGGCACGGCTTCGTTGTCGATATGAGCGACTTCTACGTCGGGGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATcc                                                                                                         <  1:367503/251‑1 (MQ=255)
                                                                                                                                              cTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGcc                                                                                                     >  1:356349/1‑251 (MQ=255)
                                                                                                                                              cTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGAACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGGTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGCAGCATCTTATCCGGcc                                                                                                     >  1:70507/1‑251 (MQ=255)
                                                                                                                                                      atagcgtttACGCGTTGATTATGGGCTGGGCGCGTGGCAACCTTTTCAGCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGAt                                                                                              <  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)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
CATTATGCGCGTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTAC  >  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: