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.

ACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTT  >  NC_000913/4307057‑4307570
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                
         cccAGGACCACCGCATCGGGATTGAACAGATTAATGCTGGTGGCAATGGCCCTTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGTCTGTTTTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCa                                                                                                                                                                                                                                        <  1:831/275‑1 (MQ=255)
               aTCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTAAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATACCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgcgc                                                                                                                                                                                                                                  <  2:283342/275‑1 (MQ=255)
                                 aaCAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAGTGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCCTTCCAGAGCAATTGGTTTCCAGGAACCCAGGATTGCCACACGCGCAGTGTTGGGCCATATCTCCCAGGGGGCTGTGACCCCGTTCGCCTGCCACACCGTGTCCACCCGTCCAAGGCGCACCGTTCATCcacccct                                                                                                                                                                                                                 >  1:64439/1‑270 (MQ=255)
                                                gtggCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTACCACACGCCCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCAAGTTGGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTAATCCCCACTGCGAACCCCAACCCCg                                                                                                                                                                                                  >  1:4157/1‑274 (MQ=255)
                                                 tggCAATGGCCCGTCCCCGTTTGTCAAGTCGGCCTTGGACGAAAGGGGCGTTTGCCGCATGGGGCAAAAAAGCGCGACAGGGGTAATTTCGTGGCTGTTGTTCGTACCAGCGTCTGAGCGCCATTCCGGAGCAATTGGTGTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTa                                                                                                                                                                                                <  4:368479/275‑1 (MQ=255)
                                                                     ttttCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCCCACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                                                                                                                                                                         >  2:473333/1‑230 (MQ=255)
                                                                                     gggACGAAAGGGGCGTTTCCGCATGGACGAAAGGATCGCGCAATGGGTAATTTCGGGGCTTTTGTTCGTACCAGCTCCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAgg                                                                                                                                                            <  2:430957/274‑1 (MQ=255)
                                                                                            aaGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac                                                                                                                                                                                                                     >  4:158685/1‑212 (MQ=255)
                                                                                                ggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTg                                                                                                                                                                                                                 <  2:158686/212‑1 (MQ=255)
                                                                                                   gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctacta                                                                                                                                               <  2:417239/275‑1 (MQ=255)
                                                                                                     tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACTCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATg                                                                                                                                                                                                                                                                            <  6:1028455/148‑1 (MQ=255)
                                                                                                      ttCCGCATAGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                                                                                                                                                           <  6:969776/148‑1 (MQ=255)
                                                                                                                aCGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGc                                                                                                                                                                                                                                                                 <  6:746517/148‑1 (MQ=255)
                                                                                                                   aaaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                                                                                                                                                              >  5:27812/1‑148 (MQ=255)
                                                                                                                                    gTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                                <  6:1073161/113‑1 (MQ=255)
                                                                                                                                    gTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                                >  5:1073161/1‑113 (MQ=255)
                                                                                                                                      aaTTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGt                                                                                                                                                                                                                                           >  6:908900/1‑148 (MQ=255)
                                                                                                                                             gggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                                                                                                                                                   >  2:185340/1‑165 (MQ=255)
                                                                                                                                                   ttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCg                                                                                                                                                                                                                              >  7:177610/1‑148 (MQ=255)
                                                                                                                                                      ttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                                                                                                                                                           <  8:486633/148‑1 (MQ=255)
                                                                                                                                                             cAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                                                                                                                                                 >  5:1020314/1‑87 (MQ=255)
                                                                                                                                                              aGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                                                                                                                                                          >  5:473432/1‑141 (MQ=255)
                                                                                                                                                              aGCGTCTTAGCGCCATGCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                                                                                                                                                          <  6:473432/141‑1 (MQ=255)
                                                                                                                                                               gCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                                                                                                                                                          <  8:745142/92‑1 (MQ=255)
                                                                                                                                                               gCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                                                                                                                                                          >  7:745142/1‑92 (MQ=255)
                                                                                                                                                                  tCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACg                                                                                                                                                                                                                                      >  8:708041/1‑125 (MQ=255)
                                                                                                                                                                  tCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACg                                                                                                                                                                                                                                      <  7:708041/125‑1 (MQ=255)
                                                                                                                                                                     tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGa                                                                                                                                                                                                                                                                           >  7:437624/1‑85 (MQ=255)
                                                                                                                                                                     tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGa                                                                                                                                                                                                                                                                           <  8:437624/85‑1 (MQ=255)
                                                                                                                                                                     tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAc                                                                                                                                                                                                            >  5:562660/1‑148 (MQ=255)
                                                                                                                                                                         gcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCg                                                                                                                                                                                                                                                      <  7:448999/102‑1 (MQ=255)
                                                                                                                                                                          ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                                >  7:557774/1‑75 (MQ=255)
                                                                                                                                                                          ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                                                                                                                                                           >  7:909748/1‑80 (MQ=38)
                                                                                                                                                                            aTTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                                                                                                                                                    <  5:913263/101‑1 (MQ=255)
                                                                                                                                                                             ttCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTCCACCCGTCCACGGCGCACCGTTCATCCAAACTGCGAACCCCTTCCCCGTACCGAGATAGGCCACCAGGACCAGTTGTTGCGTAAGGCGGTTTCACAATGCGTCCCAGAGCAGTTACCGGCAAATCGCGCGGGAAAGCACCACCGGACTCTTCAGGCTATTTCCAAGCTTacct                                                                     >  4:260099/1‑272 (MQ=255)
                                                                                                                                                                                  gagCAATTGGTTTCCAGGCACCCAGGATTGCCAAACGCGAAGTGTTGGGTCATATCTCCAAGGGGGATATGGCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCgag                                                                <  2:7458/275‑1 (MQ=255)
                                                                                                                                                                                   agCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAcc                                                                                                                                                                                              >  8:1054041/1‑148 (MQ=255)
                                                                                                                                                                                       aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                                                                                                                                                          <  7:1054041/148‑1 (MQ=255)
                                                                                                                                                                                         tGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCa                                                                                                                                                                                                                                                >  8:882908/1‑92 (MQ=255)
                                                                                                                                                                                         tGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCg                                                                                                                                                                                             >  8:1106480/1‑143 (MQ=255)
                                                                                                                                                                                          ggTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                                                                                                                                                 >  7:1014645/1‑58 (MQ=255)
                                                                                                                                                                                            tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                                                                                                                                                    <  5:124431/85‑1 (MQ=255)
                                                                                                                                                                                            tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                                                                                                                                                    >  6:124431/1‑85 (MQ=255)
                                                                                                                                                                                            tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgc                                                                                                                                                                                     <  8:413636/148‑1 (MQ=255)
                                                                                                                                                                                              tCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgt                                                                                                                                                                                                                                                     <  8:1042931/82‑1 (MQ=255)
                                                                                                                                                                                                 agggccccAGTATTGCACACCGCGCCGTGTTGGGTCATATCCCCCAGGGGGATAGACCCCAGTCCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATccg                                                  <  4:129307/269‑1 (MQ=255)
                                                                                                                                                                                                      cccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGa                                                                                                                                                                              <  7:186065/145‑1 (MQ=255)
                                                                                                                                                                                                       ccAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc                                                                                                                                                                                                                                                                          <  7:989088/52‑1 (MQ=18)
                                                                                                                                                                                                        cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                                                                                                                                                                <  8:525671/61‑1 (MQ=255)
                                                                                                                                                                                                        cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                                                                                                                                                                >  7:525671/1‑61 (MQ=255)
                                                                                                                                                                                                        cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAg                                                                                                                                                                         >  7:875706/1‑148 (MQ=255)
                                                                                                                                                                                                         aGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATcc                                                                                                                                                                                                                       <  6:270207/101‑1 (MQ=255)
                                                                                                                                                                                                          ggATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTTCCAGGAGAGTTGCAGGTTACAGTCGCGGGAAAACTCAACCGGACAATTTAGCGTATTTTCGCGCTTATCGGCGCGACCATATAAATCCGCCGCTGTTaa                                        >  2:394992/1‑275 (MQ=255)
                                                                                                                                                                                                           gATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttg                                                                                                                                                                      <  6:562660/148‑1 (MQ=255)
                                                                                                                                                                                                            aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgt                                                                                                                                                                     <  5:908900/148‑1 (MQ=255)
                                                                                                                                                                                                            aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgt                                                                                                                                                                     <  6:802576/148‑1 (MQ=255)
                                                                                                                                                                                                             ttGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                                                                                                                                                    <  8:177610/148‑1 (MQ=255)
                                                                                                                                                                                                                 cacaCGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                                                                                                                                                                                                                                     >  8:17300/1‑47 (MQ=255)
                                                                                                                                                                                                                   cacGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGaa                                                                                                                                                                             >  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)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                
ACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTT  >  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: