Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,497,385 4 bp→CAA coding (530‑533/1656 nt) opgD → glucans biosynthesis protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,497,3850TC91.7% 106.4 / 1.4 36F177S (TTT→TCT) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (21/12);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01
*NC_0009131,497,3860T.94.3% 134.4 / 5.0 36coding (531/1656 nt)opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (21/12);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01
*NC_0009131,497,3870CA94.4% 112.1 / ‑0.1 36L178N (CTC→AAC) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (21/13);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01
*NC_0009131,497,3880TA94.4% 112.3 / ‑0.1 36L178N (CTC→AAC) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (21/13);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

GTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATTGTGAGAAAAGTC  >  NC_000913/1497151‑1497615
                                                                                                                                                                                                                                          |||                                                                                                                                                                                                                                    
gTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTAtt                                                                                                                                                                                                                       >  1:662282/1‑251 (MQ=255)
        gACGCTCAGTTCTTCCATATGGGAATGTGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGcc                                                                                                                                                                                                               <  2:76708/251‑1 (MQ=255)
         aCGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGTGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCg                                                                                                                                                                                                              >  2:424982/1‑251 (MQ=255)
               aGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTtgatg                                                                                                                                                                                                        <  1:147554/251‑1 (MQ=255)
                  tcttcCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTgatgat                                                                                                                                                                                                      <  2:303694/250‑1 (MQ=255)
                       cATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACata                                                                                                                                                                                                 <  1:352044/251‑1 (MQ=255)
                          atgggaatgggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCaa                                                                                                                                                                                             <  1:637693/251‑1 (MQ=255)
                                   ggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTg                                                                                                                                                                                    >  1:340337/1‑251 (MQ=255)
                                        ccGTCGCCGCGTTCGTATGGTTTCTGTAGATCCAGCAACACAGCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGc                                                                                                                                                                               <  1:482767/251‑1 (MQ=255)
                                               cgagTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGCGTTATTTCCGCGCCGTTGCTGATACATATCAATACGGTTTGTCGTCCCGCGGc                                                                                                                                                                        >  1:250194/4‑251 (MQ=255)
                                                              tCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGGGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGCTCGACACt                                                                                                                                                         >  1:346968/1‑251 (MQ=255)
                                                                     aTCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGGGATGTAGTTTCATC‑AACGGCGCGAGTTATTTCCGCGCCGGTGATGATACATATCAATACGGTTTGGCGGCCCTCGGCCTGGCGATCGACACTTACACCg                                                                                                                                                  <  1:311538/251‑1 (MQ=255)
                                                                            aaCACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGCCACTTACACCGACAGTa                                                                                                                                            >  1:362768/1‑250 (MQ=255)
                                                                                caTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCACTACGGTTTGTAGGCCCGCGACCTTGCTATCGACCCTTACACCGAAAGTAaaga                                                                                                                                        >  2:786348/1‑250 (MQ=255)
                                                                                   cTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAaagaag                                                                                                                                      >  1:737291/1‑249 (MQ=255)
                                                                                             aaatgtccTTTCGCCCGGAGGTGTCCAAATACAACGATGCAGTGGTTGAATCAAAACAAGTAGAAGGGCAAAGCGATCTCGGCTTGTCCGGTTTCCGCGGGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATC‑CTCAACGGCTCGAGTTATCTCCGCGCCGTTGATGATACATATCAATACGGTTTGGCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCTACt                                                                                                                          <  2:479175/244‑1 (MQ=255)
                                                                                             aaaTTCACTGTCGCCCGGAGTTGTGCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCTTTTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACt                                                                                                                          <  2:382978/251‑1 (MQ=255)
                                                                                                                     taaaaaaaaaCGATGCAGGTGTTGATACAAAAAAATTAGAAGGGCAAAGCGAACTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACg                                                                                                   <  2:362768/242‑1 (MQ=255)
                                                                                                                           aCAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGCCAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAAc                                                                                            >  2:367628/1‑251 (MQ=255)
                                                                                                                           aCAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAAc                                                                                            <  1:779364/251‑1 (MQ=255)
                                                                                                                                gATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACTGTAAACCCgggg                                                                                       >  1:408868/1‑251 (MQ=255)
                                                                                                                                    cAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCa                                                                                    >  1:397477/1‑250 (MQ=255)
                                                                                                                                     aGGTGTTGATACAAAACAATTAGAAGGGCCAAGCGATCTAGGCTTTGCCGGTTTTCGCGGGTCTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAAc                                                                                  >  1:228578/1‑251 (MQ=255)
                                                                                                                                            gATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCttt                                                                           >  1:366113/1‑251 (MQ=255)
                                                                                                                                             aTACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCTGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTa                                                                          <  1:424982/251‑1 (MQ=255)
                                                                                                                                                      aaTTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATg                                                                 >  1:317620/1‑251 (MQ=255)
                                                                                                                                                           gAAGGGCAAAGCGATCTCGACTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACAGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGCTAGACACTTACACCGACCGTAAAGAATAGATCCCCGACTTTACCGCCTTCTGGCTTGATACGGTAAAACCGGGGGCAACTACCTTTAACGTTTATGccgtg                                                            >  2:281296/1‑247 (MQ=255)
                                                                                                                                                                   aaGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGCTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGTGCAACTACCTTTACCGTTTATGCGTTGCTCGAt                                                      >  1:345586/1‑250 (MQ=255)
                                                                                                                                                                             ggCTTGGCCGGTTTTCGGGTTTTTAAAGCCCCGGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGTCGCGAGTTATTTCCGCGCCGTTGATGATACAGATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCAtt                                          <  2:366113/251‑1 (MQ=255)
                                                                                                                                                                              gCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCAtt                                          >  1:164540/1‑250 (MQ=255)
                                                                                                                                                                                  tGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGCTACATCTCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGAGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCCGCATTACTgg                                     >  2:82084/1‑251 (MQ=255)
                                                                                                                                                                                   gCCGGTTTTCGCGTGTTTAAAGCCCACGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTggg                                    >  1:664695/1‑250 (MQ=255)
                                                                                                                                                                                       gatttCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATAATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTCATACGGGAACACCGGGAGCACCTACCTTTACCGTTTATGCGTTGCTCGCTAGCGCCAGAATTACTGGTGCCt                                >  2:720179/3‑251 (MQ=255)
                                                                                                                                                                                                   ttAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACt                     <  1:504705/250‑1 (MQ=255)
                                                                                                                                                                                                               aaCTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTAGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATtgtg           >  2:813221/1‑248 (MQ=255)
                                                                                                                                                                                                                       cgcCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATTGTGAGAAAAGTc  >  2:453590/1‑249 (MQ=255)
                                                                                                                                                                                                                                          |||                                                                                                                                                                                                                                    
GTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATTGTGAGAAAAGTC  >  NC_000913/1497151‑1497615

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