Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 5,442 T→C L70P (CTA→CCA)  yaaX → DUF2502 domain‑containing protein YaaX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,4420TC94.7% 92.6 / ‑1.9 38L70P (CTA→CCA) yaaXDUF2502 domain‑containing protein YaaX
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (17/19);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

CGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCA  >  NC_000913/5199‑5679
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
cGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACgg                                                                                                                                                                                                                                             >  1:326569/1‑250 (MQ=255)
        aaCGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCgccgcc                                                                                                                                                                                                                                    >  1:142019/1‑251 (MQ=255)
           gggcaaataaaaaaGGAGTAACCTATGAAAAAGATGCAATCTATCGTACCCGCACTTTCCCTGTTTCTGGGCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGGCCCGGCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGGTAGTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCAcc                                                                                                                                                                                                                                 <  2:61838/241‑1 (MQ=255)
             gCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGACAATGGCGCGCCCACGGCTGGTGGAAACAACCTTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCCCCAccgc                                                                                                                                                                                                                               >  2:384557/1‑251 (MQ=255)
                              aaCCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCCATCGCTGGCACCCACACGGACCGCCGCCCCCGCCGCGCCCCCATA‑AGaa                                                                                                                                                                                                              >  1:366909/1‑250 (MQ=255)
                                 ctatgAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCCGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCAGA‑AGAAAGcc                                                                                                                                                                                                          >  2:16910/4‑250 (MQ=255)
                                   atgAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGGAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGctcc                                                                                                                                                                                                        >  1:4796/2‑251 (MQ=255)
                                     gAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCAGGGCAATCGCGGGCCCCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGctcctc                                                                                                                                                                                                      >  2:455782/1‑251 (MQ=255)
                                         aaGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGAGAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGa                                                                                                                                                                                                  <  1:437126/251‑1 (MQ=255)
                                           gATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGACGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGCGGCCATCGCTGGCACCCACCCGGACCGCCGCCACCCGCGCGCCACCATACAGAA‑GCTCCCCATGatc                                                                                                                                                                                                >  2:317597/1‑251 (MQ=255)
                                            aTGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAACACAGGCTGCGGAAGTTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTAACTGACGCGCCCACGGCTGGTGGAAAAAACATTATGAATGGCGAAGCAATCGCTGGCAACCCCACGGACCGCCGCCACCGCCGCACCACCCTA‑AAAAAGCTCCTCaatatc                                                                                                                                                                                                >  2:334205/1‑245 (MQ=255)
                                                   cTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCAcggc                                                                                                                                                                                        >  2:173704/1‑251 (MQ=255)
                                                     agcGTACTCGCACTTTCCCTGGTTCTGGTCGCACCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCAcgg                                                                                                                                                                                         <  1:41479/246‑1 (MQ=255)
                                                            ccGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGt                                                                                                                                                                               <  1:381450/250‑1 (MQ=255)
                                                             cgctcttccccgggTTCGGGTCGCTCCCATGGCAGCACAGGCTGGGGAAATTACGTTAGTCCTGACAGTAAAATTACAGATAGGCGATCGTGATAATCGGTGCTATTACGGGGATGGAGGTCACTGGCGCGCCCCCGGCTGGTGAAAACAACATTATGAATGGCGAGGCAATCGCTGGCCCCCACACGGCCCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATTACGGCGGTCATGGTc                                                                                                                                                                              <  2:316470/239‑1 (MQ=255)
                                                             cGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAACCTCCTCATGATCATCAcggcggcccaggtc                                                                                                                                                                              >  1:77136/1‑243 (MQ=255)
                                                              gCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTcc                                                                                                                                                                             <  1:128342/251‑1 (MQ=255)
                                                                  tttCCCTTGTTCTGGTCCCTACCATGGCAGCACAGGCCGGGGAAATTCCGTTAGTCCCGGCAGTAAAATGTCAGATAGGCGATCGTGATAACCGTGGCTATGACTGGGATGGAGGCCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCCCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAgg                                                                                                                                                                          <  2:77136/250‑1 (MQ=255)
                                                                   tccccTGTTTCTGGCCGCTCCCATGGCAGCACAGGCTGGGGAAATTACGTTAGCCCCGTCAGTAAAATTACAGATAGGGGATCGTGATAATCGGGGCTATGACTGGGTGGGAGTTCACGGGCGCGACCACGGCTGGGGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCCCCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGc                                                                                                                                                                         <  2:20503/248‑1 (MQ=255)
                                                                               ggTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCTTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGTTCCCGGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCg                                                                                                                                                             <  1:471396/249‑1 (MQ=255)
                                                                                     ccccATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGa                                                                                                                                                     <  1:16910/250‑1 (MQ=255)
                                                                                         aTGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGGCATGGTTAAGGCAAACC‑ACCCCGCTAAATGaaaa                                                                                                                                                  >  1:346252/1‑247 (MQ=255)
                                                                                            gcagcaCAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGCAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCTCCGCTAAATGACAAA‑TGc                                                                                                                                             <  1:90595/251‑1 (MQ=255)
                                                                                                             aaTTACGGTAGTCCGGTAAGTAAAATTACAGATAGGCTATCGGGATAATCGTGGCTATTACTGGGATGGAGGTCACGGGCGCGACCACGGCTGGGGGAAACAACATTATGAAGGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCa                                                                                                                            <  2:50172/251‑1 (MQ=255)
                                                                                                                                     aTTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAAGCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAGCGCCTGATGCGACGCTGgc                                                                                                   <  1:302124/251‑1 (MQ=255)
                                                                                                                                                            aTCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCCTGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAAACCGGCATT‑CAGCGCCCGATGCGACGCTGGGGCGTCTTATCAGGCCTACGtga                                                                             >  1:308771/1‑248 (MQ=255)
                                                                                                                                                               gTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGCCCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTc                                                                         <  1:173704/251‑1 (MQ=255)
                                                                                                                                                                 ggCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGGCATGGTCCCAGCCAACACTC‑CCGCCAAATGACAAAAGGCCGGG‑AACCATCCGGCATTACA‑CGCCCGATGACACGCTGGCGCGGCTTATccagcccacggggattctt                                                                       >  2:78853/1‑233 (MQ=255)
                                                                                                                                                                 ggCTACTACTGGGATGGAGGTCACTGTCGCGCCCACGGCTGGTTGAAACAACATTATGAATGGCGAGCCAATCGCTGGCACCCACACGGACCGCCGCCTCCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCATGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTg                                                                       <  1:78303/251‑1 (MQ=255)
                                                                                                                                                                            ggATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCCTGGTCCAAGCAAACA‑ACAACGCTAACAGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAACGCCTGCTGCGGCGCCGGGGCGTCTTTTCAAGCCGAAGTGAATTCTGACATATAttt                                                             >  1:239448/1‑249 (MQ=255)
                                                                                                                                                                               tGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACCAA‑TGCCGGGGAACAATCCGGCAAT‑CAGCGCCTGATGCGACCCTGGCGCGGCTTTTCAGGCCTACGGTAATTCTGCCATATATTggatt                                                         >  1:268476/1‑247 (MQ=255)
                                                                                                                                                                                        cTGGCGCGACCACGGCGGGTGGAAACAACATTAGGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAgg                                                                                                                                                                          <  1:300988/132‑1 (MQ=255)
                                                                                                                                                                                          ggCGCGACCACGGCTGGGGGAACAAACATGATGAAGGGCGAGGCAACCGCTGGCACCCCCACGGACCGCCGCCAACGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCCCCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCtttt                                              <  1:317597/251‑1 (MQ=255)
                                                                                                                                                                                                     ggCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGAt                                   <  1:459910/251‑1 (MQ=255)
                                                                                                                                                                                                            ggAAACAACATTAGGAAGGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGt                            <  1:363565/251‑1 (MQ=255)
                                                                                                                                                                                                              aaaCAACATTATGAATGGCGAGGCAATCGCTTGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTc                          <  1:311122/251‑1 (MQ=255)
                                                                                                                                                                                                                    caTTATGAAGGGCGAGGCAACCGCGGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCAcgccg                    <  1:103156/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCGCTGGCACCCACACGGACCGCCGCCACCGCCGCACAACCGTA‑AGAAAGCTCCTCATGATCATCACGGCGGGCATGGGACAGGCAAACA‑ACACCGCCAAAAGACAAA‑AGCCGGGTAAAAAACCGGCAGT‑CAGCGCCTGATGCGACGCCGGCGCGTCATATCAAGCCTACGTGAATTATGAAATATATTGCATCTGAATGCTTTTGTAAGACGGATAAGGAGTTCAcgcccaatcagacaatggcagca  >  1:443041/1‑231 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
CGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATA‑AGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACA‑TCACCGCTAAATGACAAA‑TGCCGGGTAACAATCCGGCATT‑CAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCA  >  NC_000913/5199‑5679

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