Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,307,298 G→A 17.2% P168S (CCC→TCC)  alsK ← D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,307,2980GA17.2% 330.8 / 60.0 134P168S (CCC→TCC) alsKD‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (58/53);  new base A (12/11);  total (70/64)
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.54e-01

Only 100 of 149 total aligned reads displayed.

ACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTT‑TCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTT  >  NC_000913/4307057‑4307570
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                
         cccAGGACCACCGCATCGGGATTGAACAGATTAATGCTGGTGGCAATGGCCCTTGCCGCGTTT‑TCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGTCTGTTTTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCa                                                                                                                                                                                                                                        <  1:831/275‑1 (MQ=255)
               aTCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTT‑TAAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATACCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgcgc                                                                                                                                                                                                                                  <  2:283342/275‑1 (MQ=255)
                                 aaCAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTT‑TCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAGTGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCCTTCCAGAGCAATTGGTTTCCAGGAACCCAGGATTGCCACA‑CGCGCAGTGTTGGGCCATATCTCCCAGGGGGCTGTGACCCCGTTCGCCTGCCACACCGTGTCCACCCGTCCAAGGCGCACCGTTCATCcacccct                                                                                                                                                                                                                 >  1:64439/1‑270 (MQ=255)
                                                gtggCAATGGCCCGTGCCGCGTTT‑TCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTACCACA‑CGCCCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCAAGTTGGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTAATCCCCACTGCGAACCCCAACCCCg                                                                                                                                                                                                  >  1:4157/1‑274 (MQ=255)
                                                 tggCAATGGCCCGT‑CCCCGTTTGTCAAGTCGGCCTTGGACGAAAGGGGCGTTTGCCGCATGGGGCAAAAAAGCGCGACAGGGGTAATTTCGTGGCTGTTGTTCGTACCAGCGTCTGAGCGCCATTCCGGAGCAATTGGTGTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTa                                                                                                                                                                                                <  4:368479/275‑1 (MQ=255)
                                                                     ttt‑tCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCCCA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                                                                                                                                                                         >  2:473333/1‑230 (MQ=255)
                                                                                     gggACGAAAGGGGCG‑TTTCCGCATGGACGAAAGGATCGCGCAATGGGTAATTTCGGGGCTTTTGTTCGTACCAGCTCCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAgg                                                                                                                                                            <  2:430957/274‑1 (MQ=255)
                                                                                            aaGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac                                                                                                                                                                                                                     >  4:158685/1‑212 (MQ=255)
                                                                                                ggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTg                                                                                                                                                                                                                 <  2:158686/212‑1 (MQ=255)
                                                                                                   gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctacta                                                                                                                                               <  2:417239/275‑1 (MQ=255)
                                                                                                     tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACTCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATg                                                                                                                                                                                                                                                                            <  6:1028455/148‑1 (MQ=255)
                                                                                                      ttCCGCATAGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                                                                                                                                                           <  6:969776/148‑1 (MQ=255)
                                                                                                                aCGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGc                                                                                                                                                                                                                                                                 <  6:746517/148‑1 (MQ=255)
                                                                                                                   aaaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                                                                                                                                                              >  5:27812/1‑148 (MQ=255)
                                                                                                                                    gTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                                <  6:1073161/113‑1 (MQ=255)
                                                                                                                                    gTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                                >  5:1073161/1‑113 (MQ=255)
                                                                                                                                      aaTTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGt                                                                                                                                                                                                                                           >  6:908900/1‑148 (MQ=255)
                                                                                                                                             gggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                                                                                                                                                   >  2:185340/1‑165 (MQ=255)
                                                                                                                                                   ttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCg                                                                                                                                                                                                                              >  7:177610/1‑148 (MQ=255)
                                                                                                                                                      ttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                                                                                                                                                           <  8:486633/148‑1 (MQ=255)
                                                                                                                                                             cAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                                                                                                                                                 >  5:1020314/1‑87 (MQ=255)
                                                                                                                                                              aGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                                                                                                                                                          >  5:473432/1‑141 (MQ=255)
                                                                                                                                                              aGCGTCTTAGCGCCATGCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                                                                                                                                                          <  6:473432/141‑1 (MQ=255)
                                                                                                                                                               gCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                                                                                                                                                          <  8:745142/92‑1 (MQ=255)
                                                                                                                                                               gCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                                                                                                                                                          >  7:745142/1‑92 (MQ=255)
                                                                                                                                                                  tCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACg                                                                                                                                                                                                                                      >  8:708041/1‑125 (MQ=255)
                                                                                                                                                                  tCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACg                                                                                                                                                                                                                                      <  7:708041/125‑1 (MQ=255)
                                                                                                                                                                     tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGa                                                                                                                                                                                                                                                                           >  7:437624/1‑85 (MQ=255)
                                                                                                                                                                     tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGa                                                                                                                                                                                                                                                                           <  8:437624/85‑1 (MQ=255)
                                                                                                                                                                     tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAc                                                                                                                                                                                                            >  5:562660/1‑148 (MQ=255)
                                                                                                                                                                         gcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCg                                                                                                                                                                                                                                                      <  7:448999/102‑1 (MQ=255)
                                                                                                                                                                          ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                                >  7:557774/1‑75 (MQ=255)
                                                                                                                                                                          ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                                                                                                                                                           >  7:909748/1‑80 (MQ=38)
                                                                                                                                                                            aTTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                                                                                                                                                    <  5:913263/101‑1 (MQ=255)
                                                                                                                                                                             ttCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTCCACCCGTCCACGGCGCACCGTTCATCCAAACTGCGAACCCCTTCCCCGTACCGAGATAGGCCACCAGGACCAGTTGTTGCGTAAGGCGGTTTCACAATGCGTCCCAGAGCAGTTACCGGCAAATCGCGCGGGAAAGCACCACCGGACTCTTCAGGCTATTTCCAAGCTTacct                                                                     >  4:260099/1‑272 (MQ=255)
                                                                                                                                                                                  gagCAATTGGTTTCCAGGCACCCAGGATTGCCAAA‑CGCGAAGTGTTGGGTCATATCTCCAAGGGGGATATGGCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCgag                                                                <  2:7458/275‑1 (MQ=255)
                                                                                                                                                                                   agCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAcc                                                                                                                                                                                              >  8:1054041/1‑148 (MQ=255)
                                                                                                                                                                                       aTTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                                                                                                                                                          <  7:1054041/148‑1 (MQ=255)
                                                                                                                                                                                         tGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCa                                                                                                                                                                                                                                                >  8:882908/1‑92 (MQ=255)
                                                                                                                                                                                         tGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCg                                                                                                                                                                                             >  8:1106480/1‑143 (MQ=255)
                                                                                                                                                                                          ggTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                                                                                                                                                 >  7:1014645/1‑58 (MQ=255)
                                                                                                                                                                                            tttCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                                                                                                                                                    <  5:124431/85‑1 (MQ=255)
                                                                                                                                                                                            tttCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                                                                                                                                                    >  6:124431/1‑85 (MQ=255)
                                                                                                                                                                                            tttCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgc                                                                                                                                                                                     <  8:413636/148‑1 (MQ=255)
                                                                                                                                                                                              tCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgt                                                                                                                                                                                                                                                     <  8:1042931/82‑1 (MQ=255)
                                                                                                                                                                                                 agggccccAGTATTG‑CACACCGCGCCGTGTTGGGTCATATCCCCCAGGGGGATAGACCCCAGTCCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATccg                                                  <  4:129307/269‑1 (MQ=255)
                                                                                                                                                                                                      cccAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGa                                                                                                                                                                              <  7:186065/145‑1 (MQ=255)
                                                                                                                                                                                                       ccAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc                                                                                                                                                                                                                                                                          <  7:989088/52‑1 (MQ=18)
                                                                                                                                                                                                        cAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                                                                                                                                                                <  8:525671/61‑1 (MQ=255)
                                                                                                                                                                                                        cAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                                                                                                                                                                >  7:525671/1‑61 (MQ=255)
                                                                                                                                                                                                        cAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAg                                                                                                                                                                         >  7:875706/1‑148 (MQ=255)
                                                                                                                                                                                                         aGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATcc                                                                                                                                                                                                                       <  6:270207/101‑1 (MQ=255)
                                                                                                                                                                                                          ggATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTTCCAGGAGAGTTGCAGGTTACAGTCGCGGGAAAACTCAACCGGACAATTTAGCGTATTTTCGCGCTTATCGGCGCGACCATATAAATCCGCCGCTGTTaa                                        >  2:394992/1‑275 (MQ=255)
                                                                                                                                                                                                           gATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttg                                                                                                                                                                      <  6:562660/148‑1 (MQ=255)
                                                                                                                                                                                                            aTTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgt                                                                                                                                                                     <  5:908900/148‑1 (MQ=255)
                                                                                                                                                                                                            aTTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgt                                                                                                                                                                     <  6:802576/148‑1 (MQ=255)
                                                                                                                                                                                                             ttGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                                                                                                                                                    <  8:177610/148‑1 (MQ=255)
                                                                                                                                                                                                                 caca‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                                                                                                                                                                                                                                     >  8:17300/1‑47 (MQ=255)
                                                                                                                                                                                                                   ca‑cGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGaa                                                                                                                                                                             >  7:1076045/1‑133 (MQ=255)
                                                                                                                                                                                                                        cgcAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAgg                                                                                                                                                                                      >  6:28711/1‑120 (MQ=255)
                                                                                                                                                                                                                        cgcAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAgg                                                                                                                                                                                      <  5:28711/120‑1 (MQ=255)
                                                                                                                                                                                                                        cgcAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAg                                                                                                                                                                                       <  6:1020246/119‑1 (MQ=255)
                                                                                                                                                                                                                          cAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                                                                                                                                                          >  7:827621/1‑114 (MQ=255)
                                                                                                                                                                                                                          cAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgt                                                                                                                                                                     >  5:766839/1‑135 (MQ=255)
                                                                                                                                                                                                                            gtgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAg                                                                                                                                                                         >  6:823654/1‑129 (MQ=255)
                                                                                                                                                                                                                            gtgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAg                                                                                                                                                                         <  5:823654/129‑1 (MQ=255)
                                                                                                                                                                                                                            gtgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCAGCCGCTGTTAACGGCAGGTTAGGCGTAg                       >  2:66955/1‑275 (MQ=255)
                                                                                                                                                                                                                                   gTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccg                                                                                                                                                                                   >  8:95872/1‑112 (MQ=255)
                                                                                                                                                                                                                                    tCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc                                                                                                                                                                                                      >  6:1003412/1‑92 (MQ=255)
                                                                                                                                                                                                                                    tCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc                                                                                                                                                                                                      <  5:1003412/92‑1 (MQ=255)
                                                                                                                                                                                                                                    tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATa                                                                                                                                                                                        >  5:1100344/1‑106 (MQ=39)
                                                                                                                                                                                                                                    tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATa                                                                                                                                                                                        <  6:1100344/106‑1 (MQ=39)
                                                                                                                                                                                                                                     cATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATcc                                                                                                                                                                                                     <  8:684445/92‑1 (MQ=255)
                                                                                                                                                                                                                                     cATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTc                                                                                                                                                    <  7:144899/141‑1 (MQ=255)
                                                                                                                                                                                                                                     cATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg                                                                                                                                             >  6:583110/1‑148 (MQ=255)
                                                                                                                                                                                                                                       tatCTCCCAGGGGGATATGACCCAGTTCGCTTGCCACACCGTGTGTACCCCTCCCCGGCGCACCGTTCATCCACACTGCGCACCCCATCCCCGTACCTAGACAAGCCGCCTCAACCGGTTGTTTCGCAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCCGTTTAACGTCGCGGGAACACTCTTCAGGTCACTTCCCTGTACTTGCGAGCTTATCTGCGAGAGCACATAAGTCCGCCGCTGTTCACGGCAGCGTACGCGTaataaatatgtga            >  4:524449/1‑263 (MQ=255)
                                                                                                                                                                                                                                       tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATg               >  4:189713/1‑272 (MQ=255)
                                                                                                                                                                                                                                       tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGACCTTATCGGCGAGATCATATAAATCCGCCCATGTTAACGGCAGGTTAGGCGTAGAcatacttggt            >  4:460136/1‑269 (MQ=255)
                                                                                                                                                                                                                                       tatCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                                                                                                                                                          >  6:483764/1‑69 (MQ=38)
                                                                                                                                                                                                                                         tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttg                                                                                                                                                                      >  8:201585/1‑119 (MQ=255)
                                                                                                                                                                                                                                         tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc                                                                                                                                         >  5:629118/1‑148 (MQ=255)
                                                                                                                                                                                                                                         tcccccTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc                                                                                                                                         >  7:982698/4‑148 (MQ=255)
                                                                                                                                                                                                                                          ttcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCa                                                                                                                                        >  5:812042/2‑148 (MQ=255)
                                                                                                                                                                                                                                          ctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCa                                                                                                                                        >  8:29677/1‑148 (MQ=255)
                                                                                                                                                                                                                                          ctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCa                                                                                                                                        >  5:502906/1‑148 (MQ=255)
                                                                                                                                                                                                                                           ccccTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                                                                                                                                                                         >  7:282545/2‑66 (MQ=39)
                                                                                                                                                                                                                                            cccTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAgg                                                                                                                                      <  6:8767/148‑1 (MQ=255)
                                                                                                                                                                                                                                            cccAGGGGGATATGACCCAGTTCTCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAgg                                                                                                                                      <  8:793187/148‑1 (MQ=255)
                                                                                                                                                                                                                                            cccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAgg                                                                                                                                      >  5:119762/1‑148 (MQ=255)
                                                                                                                                                                                                                                             ccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATa                                                                                                                                                                                        >  7:392137/1‑97 (MQ=255)
                                                                                                                                                                                                                                              cAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGt                                                                                                                                                        >  6:455172/1‑128 (MQ=255)
                                                                                                                                                                                                                                              cAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgag                                                                                                                                    <  8:875706/148‑1 (MQ=255)
                                                                                                                                                                                                                                               tgggagATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac                                                                                                                                                                                                                     <  7:918255/65‑1 (MQ=255)
                                                                                                                                                                                                                                               tgggagATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac                                                                                                                                                                                                                     >  8:918255/2‑66 (MQ=255)
                                                                                                                                                                                                                                                gggggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGttt                                                                                                                                                      <  5:458049/128‑1 (MQ=255)
                                                                                                                                                                                                                                                gggggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGt                                                                                                                                                        <  6:1044221/126‑1 (MQ=255)
                                                                                                                                                                                                                                                gggggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGt                                                                                                                                                        >  5:1044221/1‑126 (MQ=255)
                                                                                                                                                                                                                                                gggggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTcc                                                                                                                                          <  8:1141413/140‑1 (MQ=255)
                                                                                                                                                                                                                                                gggggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTcc                                                                                                                                          >  7:1141413/1‑140 (MQ=255)
                                                                                                                                                                                                                                                 ggggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCt   >  2:330960/1‑274 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                
ACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTT‑TCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACA‑CGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTT  >  NC_000913/4307057‑4307570

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑