Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 25,685 G→A R160K (AGA→AAA)  lspA → lipoprotein signal peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091325,6850GA97.8% 133.6 / ‑4.4 45R160K (AGA→AAA) lspAlipoprotein signal peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (25/19);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

TATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGC  >  NC_000913/25440‑25931
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
tATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑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‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGga‑tg                                                                                                                                                                                                                        <  1:155693/247‑1 (MQ=255)
                                   cGAAGGCCACGCAGAAGCTAAACAATA‑TAGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACATTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACATTAATAAACCCTGCCggctgagctgctg                                                                                                                                                                                                                 >  1:313479/1‑240 (MQ=255)
                                                gAAGCGAAAAAAAA‑TCGCTGACGCGGTGATT‑TTCGGCGGCGGGTTGGGCAGCCTGTAGGACCGCCTGTGGCACGGCTTCGTTGTCGATATGAGCGACTTCTACGTCGGGGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATcc                                                                                                                                                                                                   <  1:367503/251‑1 (MQ=255)
                                                    cTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGcc                                                                                                                                                                                               >  1:356349/1‑251 (MQ=255)
                                                    cTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGAACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGGTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGCAGCATCTTATCCGGcc                                                                                                                                                                                               >  1:70507/1‑251 (MQ=255)
                                                            atagcgtttACGCGTTGATT‑ATGGGCTGGGCGCGTGGCAACCTTTTCAGCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGAt                                                                                                                                                                                        <  1:94457/244‑1 (MQ=255)
                                                                                                 gCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGCAATCGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                    >  1:346049/1‑251 (MQ=255)
                                                                                                   atcgtgctcgcccGCCTGTGGCTCGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                                  <  2:363819/240‑1 (MQ=255)
                                                                                                      cTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGCTGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                               >  2:500353/1‑251 (MQ=255)
                                                                                                        gTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                             <  1:205020/251‑1 (MQ=255)
                                                                                                                 cctggtgccccggttCGTTGGCGATTTGATTGTCTTCTGCGTCGGTGGCTGGCTATGCGCCACCTTCAACCTTGCCGTTACTGCCATCTGTGTCGGTGCGGGAATGATTGTGCTGGAAGGTTTTTTTCCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAAGCGTAGTCCGTATAAGGCGTTTACGCCGCATCCGGCAAAAAt                                                                                                                                    <  2:490973/238‑1 (MQ=255)
                                                                                                                  cTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCAGGATAAGGCGTTTACACCGCATCCGGCAAAAATc                                                                                                                                   >  1:237425/1‑251 (MQ=255)
                                                                                                                    gtgGCACGGCTTCGTTGTCGATATGAGCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATc                                                                                                                                   <  1:406612/249‑1 (MQ=255)
                                                                                                                                cGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAagag                                                                                                                     <  2:437068/251‑1 (MQ=255)
                                                                                                                                    gTCGATATTATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCaaa                                                                                                                 <  1:378299/251‑1 (MQ=255)
                                                                                                                                       gttaggagcGACTTCTACGGCGGCGGCTGGCACTTCGCCACCTTCAACCTGGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTAGCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAACCCt                                                                                                              <  2:99220/243‑1 (MQ=255)
                                                                                                                                         tatGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGCAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCATTAAATATAAGAGAAAACCTGc                                                                                                            >  2:429861/1‑251 (MQ=255)
                                                                                                                                          atGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGc                                                                                                            <  1:234391/250‑1 (MQ=255)
                                                                                                                                          atGATCGAATTCTAAGACGGCGGCTGGAACTTCGCCACCTTCAACCTTGCCGGTACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAACAACAATAATAAACCCTGCCGGA‑AGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCAGTCACACATCCTTAAATATAAGAGCAAACCTGcc                                                                                                           >  2:336976/1‑250 (MQ=255)
                                                                                                                                          aggatcgctttCTACGTCGGCGGCTTGCACTTCGCCCCCTTCACCCTTGCCGATACTGCCATCTGTGTCGGTGCTGCACTGATTGTGCTGGAAGGTGTTTTCCCTTCTAAAGCGAAAAAACAATAATAAACCCGGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGc                                                                                                            <  2:277015/241‑1 (MQ=255)
                                                                                                                                                           cggcggCTGGCACTTCGCCACCTTCAACCTTGCCGATAATGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAGACAATAATAAACCCTGCCGGC‑TTCGATGCTGACGCATCTTATCAGGCCTACCGATTGCTGCGAACTCGTAGGCCAGATAAGGCGTTTACGCCGCCTCCTGCACACAGCCTTCAATATAACCGCCAACCTGCATGTCTGAATCAGTAc                                                                                            >  2:178241/1‑249 (MQ=255)
                                                                                                                                                                   ggCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATa                                                                                   <  1:121967/250‑1 (MQ=255)
                                                                                                                                                                     cacgtcgccccatttatcctggtcGATACTGCCCTCTGTGTCGGTGCGGCACGGATATTGTGGGAAGTTTTTTTACCTTCTAAAGCGAAAAAACAATAATAACCCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCTGCCTATAGATTGCTGCGATATCGTAGGCCTGATAAGGCGTTTACGCTTCATCCGGCAAAAATCCATAAATATAAGAGCAAACCTGCATGTCTGATTCTGTACAGAGCAATAgcg                                                                                <  2:434684/228‑1 (MQ=255)
                                                                                                                                                                       cTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTATGAATCTGTACAGAGCACTAGCGcc                                                                              >  2:22733/1‑251 (MQ=255)
                                                                                                                                                                           gCCAACTTCAACCTTGCCGATACTGCCATCTGTGTCGGGGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAAAGCAAACCTGCATGTCTTAATCTGTACAGAGCAATCGCGCCGTcc                                                                          >  2:264450/1‑251 (MQ=255)
                                                                                                                                                                                     aCCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAACACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGACTCTGTACAGAGCAATAGCGCCGTCCTGGTGCACtt                                                                >  2:399330/1‑251 (MQ=255)
                                                                                                                                                                                              aTACTGCCATCTGTTTCGGTGCGGCACTGATTGTTCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTaaa                                                       <  2:266671/251‑1 (MQ=255)
                                                                                                                                                                                                    ccATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCCTCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAACCTcg                                                  >  2:157466/1‑250 (MQ=255)
                                                                                                                                                                                                      ctctgtgtcggtaagtaaCTGATTGTGCTGTAAGGCTTTTTACCTTACAAAGCGAACAAACCATAATAAACCCTGCCGGA‑TGCGATGCTGAAGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGGTAAGGCGTTTAAGCCGCATCCGGCAAAAATCCTTAGATATAAGTGAAAACCTGCCTGTCTAAATCTGTACAAAGCCCTATCGCCGTACTGGGGCTTTTAACGCTCAAACTcgact                                               >  2:467314/18‑250 (MQ=255)
                                                                                                                                                                                                              cggagCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGCACAGAGCACAATCGCCGTCCTGGTGAACTGCAAACTACAACTCGACGATGGaacc                                       >  2:263186/5‑251 (MQ=255)
                                                                                                                                                                                                                  gCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTGGACGCTaaaaagagaagatggaaccacc                                    <  1:222826/250‑19 (MQ=255)
                                                                                                                                                                                                                  gCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAAGCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGGCCTGGTGCACTTCACTCTAAAACTCGACGATGGCACCAccg                                   <  2:295408/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGAAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCAccgccg                                >  2:224434/1‑251 (MQ=255)
                                                                                                                                                                                                                       aCTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCAccgccg                                >  1:58867/1‑249 (MQ=255)
                                                                                                                                                                                                                           aTTGTGTTGGAAGGTTTTTGCCCTTCTAAAGCAAAAAAACAATAATAACCCCGGCCGGAGTG‑GATGCTGACGCATCTTACCCGGCCTACAGATTGCTGCGAAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCAAGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCAAAAACTCGACGAGGGCACCACCGCCGAGTCTa                          <  1:277569/251‑1 (MQ=255)
                                                                                                                                                                                                                             tgtgCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGACACCACCGCCGAGTCTAcc                        >  1:440145/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ttttttACCTTCTAAAGCGAAAAAAAAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATACGGCATACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTCAAACTCGACGATTGCACAACCGCCGAGGCTACCAGCAACAACGGt            >  2:342154/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ttCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAACCCGgcg   >  2:307023/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   tCTAAAGCGAAAAAACAATACTAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTCGGCCGGAGAAGGCATTTACGCCGCATCCGGCAACAATCCTTAAATATAAGAACAAACCTGCATGTCTGAATCTGTACCGAGCAATAGCTCCGACCTGGAGCACTTCACGCTAAAACTCCACGGTGGCACCACCTCCGAGTTTACCCGCAACAACGGTAAACCGgcgc  >  2:315914/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
TATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGC  >  NC_000913/25440‑25931

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