Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 789,579 2 bp→AA coding (400‑401/1149 nt) galK ← galactokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913789,5790TA100.0% 156.0 / NA 48E134L (GAA→TTA) galKgalactokinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (25/23);  total (25/23)
*NC_000913789,5800CA97.9% 145.2 / NA 48E134L (GAA→TTA) galKgalactokinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (24/23);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTG  >  NC_000913/789347‑789818
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                              
gTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCACCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAAc                                                                                                                                                                                                                               >  2:156986/1‑251 (MQ=255)
   gatgatgaCGACAGCCACACCTTTGGGCACGGAAACTGCTTTGGTCCCCGTCGAGCGGCAATCGAACAGCAAGGCATGATCTTTCTTGCCGACCGCGGAAATTAGCTGTTCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGTTTCCTGACCGTTAAGCGCAGTTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCATTACGGTTCCGACCGCGACTACAAGTCAAGTGGAAGAACtt                                                                                                                                                                                                                             >  2:228152/1‑250 (MQ=255)
         gacgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACGTAACCCg                                                                                                                                                                                                                       <  1:306732/250‑1 (MQ=255)
                      ccTTTGGGCATGGAAACTGCCTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGTGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCGGCTGCAATACGGTTCCGACCGCGGCTAACAGTGAAGCGGAAGAACGTAACCCGTCACCCTGCGGcac                                                                                                                                                                                                         <  2:154276/251‑1 (MQ=255)
                          tGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACAtt                                                                                                                                                                                                      <  1:227583/250‑1 (MQ=255)
                                                ccccAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCAcg                                                                                                                                                                                <  2:247763/250‑1 (MQ=255)
                                                   cAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCAcgccgc                                                                                                                                                                            >  1:332180/1‑251 (MQ=255)
                                                    aGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACgccgcc                                                                                                                                                                           <  1:284290/251‑1 (MQ=255)
                                                    aGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACgccgcc                                                                                                                                                                           >  1:290193/1‑251 (MQ=255)
                                                         agcggggaTCGCTCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTCCAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAACACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGa                                                                                                                                                                         >  2:231138/8‑249 (MQ=255)
                                                             gCAATCGATCAGCAAGGCCTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCtgt                                                                                                                                                                   <  1:215434/251‑1 (MQ=255)
                                                                aTCGATCCGCAAGCCCTGATCTATCTTCCCGAGCCCGGAAATTCGCTGAACCATCATCCCGCCGTCACAGCCTACAAACTGGTCTTCTGCCTCCTGACCGTTAAGAGCGATTTGCGCGACGTCCACAGGCACATGATAACGCTGCTGCAATCCGGTTCCGACCGCGGCTAACAGTGAAGCGGCAGAACTTAACCCGGCACCCTGCGGCCCATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCtgttgt                                                                                                                                                                >  2:148550/1‑251 (MQ=255)
                                                                        gCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGt                                                                                                                                                        >  1:117570/1‑251 (MQ=255)
                                                                            ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGACGACGCTGTTGTTACGCAGTTGc                                                                                                                                                     >  2:132145/1‑250 (MQ=255)
                                                                             gcatggtgtttcttGCGGAGCGCGGAAATTAGTTGATCCCTGATCCCGCTGTTACAGCCTAAAAACTGGTTTTATGATTCCTGAACGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGAGGATCCAACTGTTGCACTTCGGTTCCGAACGCGACTAACAGTGAAGCGGAAGAACTCAACCCGGCACCCAGCGGCTCAGTGCCGCTGCTCACCCGGTCCACGCCGCCGAAGCGGGTGGTTCGCAGTTGctg                                                                                                                                                   >  2:301204/9‑249 (MQ=255)
                                                                              cATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGc                                                                                                                                                     >  2:3535/1‑248 (MQ=255)
                                                                                ggagctttcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTGAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATg                                                                                                                                                <  2:311449/247‑1 (MQ=255)
                                                                                        cttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACCAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTCCGCAGTTGCAGATGTTTcacc                                                                                                                                         >  2:14322/1‑250 (MQ=255)
                                                                                           gCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgc                                                                                                                                     >  1:91626/1‑251 (MQ=255)
                                                                                            ccGAGCGCGGAAAGTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAGGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTACCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcc                                                                                                                                    <  2:124901/251‑1 (MQ=255)
                                                                                                 cgcgGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGa                                                                                                                                >  1:273742/1‑250 (MQ=255)
                                                                                                           aGCTGAACCATGATCCCGCAGTTACAGCCTACAAACTGGTGTTCTGCTTCCTGACCGTTAAGCGCGTTTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTATCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGGTTTTGTGCGGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAg                                                                                                                      <  2:364446/250‑1 (MQ=255)
                                                                                                                 tCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTg                                                                                                               >  1:47519/1‑251 (MQ=255)
                                                                                                                  ccATGATCCCGCAGTGACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTTCGGCACATTGCCGCTGATCGCCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGa                                                                                                              <  2:288576/251‑1 (MQ=255)
                                                                                                                        tCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGAGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATCCGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGttt                                                                                                        >  1:362149/1‑251 (MQ=255)
                                                                                                                              aGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCGGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCTCCCTGCGTCACATTGCCGCTGATCGCCATGTCCGCGCCGCCGAAGGTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAtgt                                                                                                  <  1:163466/251‑1 (MQ=255)
                                                                                                                                         aCAAACTGGTTTTCTGCTTCCTGACCTTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTGCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATggg                                                                                       <  1:262191/251‑1 (MQ=255)
                                                                                                                                          cAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGCAGCTGATGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATCGGTTTCATTTGCGACAATGGgc                                                                                      >  2:128377/1‑251 (MQ=255)
                                                                                                                                                ggTTTGCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATc                                                                                 <  2:230533/250‑1 (MQ=255)
                                                                                                                                                       tGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCAGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGaa                                                                          >  1:40069/1‑250 (MQ=255)
                                                                                                                                                               gACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgtc                                                                  <  1:167990/250‑1 (MQ=255)
                                                                                                                                                               gACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCGGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgtcg                                                                 <  2:283007/251‑1 (MQ=255)
                                                                                                                                                                     taatcgcgATTTGTCCGCCGTCCAGCGCGAGTTGATAAAGCTGCTGCAATACGGTTCCGACCGTGTCTAACGGTGAAGCGGAAGAACTTAACCCGGCACCCTGCTGCAGATTGCCTCTGAGTACCCTTTCTACGGCGCCGAAGCTGTTTTGTCGCAGGTGTAGATGTTTAACGACGCCACGTGCGTAGTTAGCCCATTGAGAGTTTTCATGGGCTGGAATGGGCGCATCGAGGGAAAATGCGTCGAGccg                                                            <  2:226715/246‑3 (MQ=255)
                                                                                                                                                                       aGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGATGAAATAAGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTACCCCGGCACCCTGCGGCACCTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCCCGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGAtt                                                         >  1:125160/1‑251 (MQ=255)
                                                                                                                                                                                  tgCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCgg                                              >  2:111127/1‑251 (MQ=255)
                                                                                                                                                                                   gcgcCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGc                                             <  1:35541/251‑1 (MQ=255)
                                                                                                                                                                                       cGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCATTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGcc                                         <  1:215103/251‑1 (MQ=255)
                                                                                                                                                                                       cGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGcc                                         >  1:79111/1‑251 (MQ=255)
                                                                                                                                                                                                     tgatAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACa                                                                                            >  2:128340/1‑186 (MQ=255)
                                                                                                                                                                                                     tgatAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACa                                                                                            <  1:128340/186‑1 (MQ=255)
                                                                                                                                                                                                      gatAAAGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGATCTTAACCCGGCACCCTGCGGCACATTGCCGCTTATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTa                          <  1:228152/251‑1 (MQ=255)
                                                                                                                                                                                                         aaaGCTGCTGCAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACgg                       >  1:314783/1‑251 (MQ=255)
                                                                                                                                                                                                              ggctgcAAGACGGTGCCGACCGCGCCTAACAGTGAAGCGGAAGAACTTAACCCGGCCCCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATc                  <  2:220428/250‑1 (MQ=255)
                                                                                                                                                                                                               gctgcAACACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGAGCACCATGTCCACGCCGCCGAAGCTGTGGTTACGCAGTTGCAGAGGTGTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTGTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATc                  <  2:107572/250‑1 (MQ=255)
                                                                                                                                                                                                                 tgcAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGGACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAgaccctt                                                                                                                >  1:31156/1‑148 (MQ=255)
                                                                                                                                                                                                                   cAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTTTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCgtg             <  1:314115/251‑1 (MQ=255)
                                                                                                                                                                                                                   cAATACGGTTCCGACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCAACATGTCCACGCCGCCGGAGCTGTTGTTACGCAGTTGCAGATGCTTCCCCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGACGAGCTGGTTTTCATAATAGGCTGCCATCACGCGAACTTTACGGTCATAGCgtg             >  2:11927/1‑251 (MQ=255)
                                                                                                                                                                                                                               gACCGCGACTAACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTg  >  2:306122/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                              
GTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTG  >  NC_000913/789347‑789818

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