Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 789,580 C→A E134* (GAA→TAA)  galK ← galactokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913789,5800CA95.1% 179.6 / ‑2.3 62E134* (GAA→TAA) galKgalactokinase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (30/28);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATC  >  NC_000913/789335‑789821
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
tttGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTg                                                                                                                                                                                                                                               >  1:458099/1‑251 (MQ=255)
 ttGAAGTTACTGTTGATGATGACGACGGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTg                                                                                                                                                                                                                                               <  2:323393/250‑1 (MQ=255)
 ttGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGCTATGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGCCCGCGCCTTCCAGTg                                                                                                                                                                                                                                               >  2:17821/1‑250 (MQ=255)
  tgcaGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGAGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTTC‑GCAGTCTCGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGAGACTTACAGTGa                                                                                                                                                                                                                                              >  2:143204/4‑250 (MQ=255)
       ttACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCgg                                                                                                                                                                                                                                         >  2:117332/1‑250 (MQ=255)
            gTTGATGATGACGACAGCCACTCCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAAc                                                                                                                                                                                                                                   <  1:378363/251‑1 (MQ=255)
               gatgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGCTTTGTGC‑GCCGTCCAGCGGCAGCTGCTAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACtt                                                                                                                                                                                                                                 >  2:103199/1‑250 (MQ=255)
                          cAGCCACACCTTTGGGCATGGAAACTGCTTTGGGCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCAccc                                                                                                                                                                                                                     <  2:87098/251‑1 (MQ=255)
                           aGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCt                                                                                                                                                                                                                    <  2:84984/251‑1 (MQ=255)
                            gCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTg                                                                                                                                                                                                                   <  1:422535/251‑1 (MQ=255)
                                         gCATGGAAACTGCGTTGGGCCCCAGTGGGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCg                                                                                                                                                                                                      <  2:277415/251‑1 (MQ=255)
                                           aTGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCt                                                                                                                                                                                                    <  1:277170/251‑1 (MQ=255)
                                                    gCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAGGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCAt                                                                                                                                                                                           >  1:12494/1‑251 (MQ=255)
                                                                 gTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGccgccg                                                                                                                                                                             >  1:209529/1‑251 (MQ=255)
                                                                    aGCGGCAATCGATCAGCAATGCAGGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCTTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCGGCAATACGGGTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTATCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAg                                                                                                                                                                          <  2:250177/251‑1 (MQ=255)
                                                                    aGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGa                                                                                                                                                                            >  2:396894/1‑249 (MQ=255)
                                                                       ggCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGt                                                                                                                                                                                         <  1:162614/234‑1 (MQ=255)
                                                                       ggCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGt                                                                                                                                                                                         >  2:162614/1‑234 (MQ=255)
                                                                        gCAATCGATCAGCATGGAATTATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTTCTTCCTGACCGTGAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTCACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCtgt                                                                                                                                                                      <  2:290735/251‑1 (MQ=255)
                                                                             cGATCAGCAAGTCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGCTCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAATCGGAAGAACTTAACCCGGCACCCTGCGGCCCATTGCCGCTGATCACCATGT‑CCACGCCGCCAAAGCTgttgtt                                                                                                                                                                  >  2:231029/1‑250 (MQ=255)
                                                                              gATCGGCAAGGCATGATCTTTCTTCCCGAGCGCGGCAATTAGCTGATCCATGATCCCGCATTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGTCAGCTGATAAAGCTGCTGCAATACGGTTCCGCCCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGTACCCCTGCGGCACCTTGCCGTTGTTCACCATGT‑CCACGCCACCGACCCTgttgttc                                                                                                                                                                 >  2:253366/1‑249 (MQ=255)
                                                                                 cAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCACGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTCGTTaagc                                                                                                                                                              >  1:281228/1‑247 (MQ=255)
                                                                                  aGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGc                                                                                                                                                              >  2:56654/1‑249 (MQ=255)
                                                                                  aGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAg                                                                                                                                                            <  1:91501/251‑1 (MQ=255)
                                                                                             aTCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGTTCACCCTGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGtt                                                                                                                                                 >  2:388408/1‑251 (MQ=255)
                                                                                                  tcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTcac                                                                                                                                             <  1:244128/250‑1 (MQ=255)
                                                                                                        cGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgccc                                                                                                                                      >  1:76770/1‑250 (MQ=255)
                                                                                                          ggcgcGGAAAGTAGCTTATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCCCCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgccacg                                                                                                                                    <  2:190726/250‑1 (MQ=255)
                                                                                                             gcgGAAATTAGCTGAGCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCGGCTACCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTACCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCTCGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAAc                                                                                                                                 <  1:65862/251‑1 (MQ=255)
                                                                                                              cgGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACg                                                                                                                                >  1:458811/1‑251 (MQ=255)
                                                                                                                     tAGCTGACCCCTGATCCCGCAGTTACAGCCTACAAACGGGTTTTCTGCTTCCTGACCGTTAAGCGAGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTa                                                                                                                          <  2:9841/250‑1 (MQ=255)
                                                                                                                           aTCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCAtt                                                                                                                   <  1:19044/251‑1 (MQ=255)
                                                                                                                               aTGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGACCGAAGTTAGCCCATTGAt                                                                                                                >  1:485598/1‑250 (MQ=255)
                                                                                                                                 gATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGaa                                                                                                                                                                           >  2:185799/1‑189 (MQ=255)
                                                                                                                                 gATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGaa                                                                                                                                                                           <  1:185799/189‑1 (MQ=255)
                                                                                                                                              cAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGAc                                                                                                >  2:182753/1‑251 (MQ=255)
                                                                                                                                              cAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGAAGTTTCACCACGCCACGACCGAAGTTAGCCCATTGATAGTTTTCATGTGCgcc                                                                                                >  2:218709/1‑249 (MQ=255)
                                                                                                                                                 ccTACAAACTGGTTTTCTGCTTCCTGCTCGTTAAGCGCGCTTTTGGC‑GTCGTCCAGCGGCAGATGATAAAGCTGCTTCAATCCGGTTCCGCCCGCGACTTTCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGTTGTTGTTACGCAGTTGCAGATGTTTCCCCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCTACAAt                                                                                             <  2:331270/251‑1 (MQ=255)
                                                                                                                                                 ccTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAAt                                                                                             >  2:59286/1‑251 (MQ=255)
                                                                                                                                                       aaCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGgcgc                                                                                       <  1:56654/251‑1 (MQ=255)
                                                                                                                                                           ggtgttCTGCTTCCTGACCGTTAAGCGCGATT‑GTGCCGCCGCCCAGCGGCAGATGATAAAGATGCTGCAATACGGTTCCGTCCGCGACTTACAGTGAAGCGGAAGAACTTAGCCCGGCACCCTGTGGCACATTGCCGCTGATCCCCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCCCCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATc                                                                                    <  2:383490/246‑1 (MQ=255)
                                                                                                                                                           ggTTTTCTGCTTCATGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGGTGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCAGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGCACGTAGTTAGCCCATTGATAGTTTTCATTGGCGTCAATGGGCGAATCg                                                                                   >  2:230619/1‑251 (MQ=255)
                                                                                                                                                            gTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCg                                                                                   >  2:476162/1‑250 (MQ=255)
                                                                                                                                                                tCTGCTTCCTGACCTTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAggg                                                                               <  1:304317/250‑1 (MQ=255)
                                                                                                                                                                     ttCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCAGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGCCAATGGGCGCATCGAGGGACAACt                                                                         >  1:415949/1‑251 (MQ=255)
                                                                                                                                                                         tGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgtc                                                                     >  1:377405/1‑251 (MQ=255)
                                                                                                                                                                                   gcgcGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGACAACTCGTCGAGCTGAttt                                                           >  2:191097/1‑251 (MQ=255)
                                                                                                                                                                                        aTTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCAGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGa                                                                              <  1:333007/227‑1 (MQ=255)
                                                                                                                                                                                        aTTTGTGC‑GACGTCCAGCGGCAGATGATAAAGCTGCAGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCAACCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGCACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGATGGa                                                                              >  2:333007/1‑227 (MQ=255)
                                                                                                                                                                                            gtgC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCg                                                  <  1:463265/251‑1 (MQ=255)
                                                                                                                                                                                               c‑gcCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCt                                               <  1:103199/251‑1 (MQ=255)
                                                                                                                                                                                                 gcCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGCGCTGATTTTCATAATCGGCTg                                              >  1:268187/1‑251 (MQ=255)
                                                                                                                                                                                                  ccGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGc                                             >  2:79379/1‑251 (MQ=255)
                                                                                                                                                                                                   cGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGcc                                            <  1:226940/251‑1 (MQ=255)
                                                                                                                                                                                                    gTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTAAGCAGTTGCAGATGTTTCACCACGCCACGAACGTAATTAGCCCATTGATAGTTTTCATGTGCGACACTGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATcagccccc                                            >  2:232793/1‑243 (MQ=255)
                                                                                                                                                                                                       cAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCa                                        <  1:476162/251‑1 (MQ=255)
                                                                                                                                                                                                        aGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCAc                                       <  1:117332/251‑1 (MQ=255)
                                                                                                                                                                                                         gCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACCTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGAACACCATTT‑CCACGCCGCCGAAGCTTTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCAcg                                      <  1:17821/251‑1 (MQ=255)
                                                                                                                                                                                                            ccAGATGATAACGCTGCTGCCATACGGTTCCGACCGCGACTTACCGTGAAGCGGACGACCTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCCCGACGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGCCATTGGGCGCATCGAGGGAAAACTCGCCGAGCTGATTTTCATAATAGGCTGCCATCCcgcg                                    >  2:40206/2‑250 (MQ=255)
                                                                                                                                                                                                                                  tACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCgtggtg             >  2:437828/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cGGTTCCGACCGCGACCTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGTGGCACATTTCCGCTGATCACCA‑GTGCCACGCCGCCGAAGCTGGTGTTCCGCAGTTTCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTGTCATGTGCGACAATGGGTGCATCGAGGTAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGca           <  2:253987/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ccGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATc  <  1:437828/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
TTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGC‑GCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGT‑CCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATC  >  NC_000913/789335‑789821

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