Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 790,832 Δ1 bp coding (198/1047 nt) galT ← galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913790,8320G.98.4% 281.2 / ‑4.3 62coding (198/1047 nt)galTgalactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (30/31);  total (31/31)
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

GCACCCATGGGTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTG  >  NC_000913/790585‑791077
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
gCACCCATGGGTACGGTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCGAAGTCATTAGTGAAAACGTAAGTCCCGGTGGAATCGGGtttt                                                                                                                                                                                                                                                    <  2:246297/251‑1 (MQ=255)
         ggTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTCTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAATTCCCGGTGTAAGCGGGTTTTTATCGCCtt                                                                                                                                                                                                                                          >  2:103599/1‑250 (MQ=255)
          gTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGt                                                                                                                                                                                                                                         <  1:443405/251‑1 (MQ=255)
            aCGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCa                                                                                                                                                                                                                                       <  1:179042/251‑1 (MQ=255)
            aCGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGATGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTATCGCCTGTc                                                                                                                                                                                                                                        >  1:27996/1‑250 (MQ=255)
                               ggTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTGTTACTGTGATCCGGTGAAAAGCAGATCGCCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGTCATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCa                                                                                                                                                                                                                    <  1:178137/251‑1 (MQ=255)
                                gTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTgcgc                                                                                                                                                                                                                     >  2:173987/1‑249 (MQ=255)
                                gTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCGCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGGGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCag                                                                                                                                                                                                                   <  2:240057/251‑1 (MQ=255)
                                                     aCGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTAGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTACCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACATGCGCAGAGGAAGCAATCTGGATCGTgc                                                                                                                                                                                              >  2:372807/1‑251 (MQ=255)
                                                     aCGATTCCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGGCACCCGCACATTCCCTGCGCAGAGGAAGCAATCGGGATCGTgc                                                                                                                                                                                              <  1:166737/251‑1 (MQ=255)
                                                       gATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTgcgc                                                                                                                                                                                            >  2:87003/1‑251 (MQ=255)
                                                          ttCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAg                                                                                                                                                                                          <  1:332055/250‑1 (MQ=255)
                                                              gTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAac                                                                                                                                                                                     >  1:103137/1‑251 (MQ=255)
                                                               tCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAac                                                                                                                                                                                     >  2:106876/1‑250 (MQ=255)
                                                                       gcAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTtg                                                                                                                                                                            >  1:453990/1‑251 (MQ=255)
                                                                               gAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCACGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATGACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCg                                                                                                                                                                     <  2:453990/250‑1 (MQ=255)
                                                                                 gCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGttt                                                                                                                                                                  <  1:263000/251‑1 (MQ=255)
                                                                                      ggCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCGGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGGTCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCTCCCGCACATTGCCTGCGCAGAGGAAGCAATCTGAATCGTGCGCAGGTAACACCTGTTTGGCTGGTTTTTCCt                                                                                                                                                               <  2:295180/249‑1 (MQ=255)
                                                                                          gCGTTTTACTGTGATCCGGTGAAAAGCAGAGCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGccc                                                                                                                                                         <  2:17488/251‑1 (MQ=255)
                                                                                               ttACTGTGATCCGGCGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGc                                                                                                                                                    >  1:380481/1‑251 (MQ=255)
                                                                                                      gATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAgggg                                                                                                                                              <  2:112442/250‑1 (MQ=255)
                                                                                                          cGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCtt                                                                                                                                         >  1:171779/1‑251 (MQ=255)
                                                                                                              gAAAATCAGATCACCCGGCTGGGGCCGCGCGCGCTCTGGCAACGCATCAGCTGCTCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGGAAGGCCCGGTGTAATCGGGTGTTTATCGCCTGTCACCCGCACATTACCTGCGCAGATGAAGCAAGCTGGATCGGGCGCAGGCATCACCTGTTTGGCTGGCGTTTCCTgc                                                                                                                                                             <  2:470280/227‑1 (MQ=255)
                                                                                                              gAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTgc                                                                                                                                                             >  1:470280/1‑227 (MQ=255)
                                                                                                               aaaaGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGcac                                                                                                                                                                                                                                >  1:262983/1‑159 (MQ=255)
                                                                                                               aaaaGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGcac                                                                                                                                                                                                                                <  2:262983/159‑1 (MQ=255)
                                                                                                                aaaGGAGATCAGCCGGATGGGGCCGCGCTCGCTCTGGCACCCAACCCGCGGAGCATAACTTTCTGGCGCATCGGGCGTGTCAGACATCAAAGCCGACAAGTCCTTAGTGACAACGTCAGTCCAGGTGTAATCGGGTTTTTCACGCATGTCACCCGCACCTTACATGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAGCACCTGTTTGGCTGGCGTGTCCTGCGCCCCCGGccgggcgcgattaacggg                                                                                                                                   >  2:93307/1‑235 (MQ=255)
                                                                                                                   gCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCATGGGCGCTTAGCgcggtg                                                                                                                                <  1:107261/251‑1 (MQ=255)
                                                                                                                   gCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTa                                                                                                                                        <  2:4615/243‑1 (MQ=255)
                                                                                                                   gCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTa                                                                                                                                        >  1:4615/1‑243 (MQ=255)
                                                                                                                    cAGATCACCCGTCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCgcggtgc                                                                                                                               <  2:181703/251‑1 (MQ=255)
                                                                                                                     aGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACCTTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCgcggtgcg                                                                                                                              >  2:326538/1‑251 (MQ=255)
                                                                                                                     aGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCCGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCgcggtgcg                                                                                                                              >  2:355924/1‑251 (MQ=255)
                                                                                                                         cACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGcggtgcggtg                                                                                                                           <  1:173987/250‑1 (MQ=255)
                                                                                                                              ggCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCt                                                                                                                                                                                 <  2:115823/191‑1 (MQ=255)
                                                                                                                              ggCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCt                                                                                                                                                                                 >  1:115823/1‑191 (MQ=255)
                                                                                                                               gCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAAcc                                                                                                                      >  1:164632/1‑249 (MQ=255)
                                                                                                                                              gcTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCggg                                                                                                     >  1:50531/1‑250 (MQ=255)
                                                                                                                                                   ggCAACGCATCCGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAg                                                                                                  <  2:138545/249‑1 (MQ=255)
                                                                                                                                                   ggCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAg                                                                                                  <  2:50531/249‑1 (MQ=255)
                                                                                                                                                     cAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACTCCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCtt                                                                                                                                         <  2:287771/208‑1 (MQ=255)
                                                                                                                                                     cAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCtt                                                                                                                                         >  1:287771/1‑208 (MQ=255)
                                                                                                                                                         gCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCTCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGt                                                                                          <  1:299359/251‑1 (MQ=255)
                                                                                                                                                              aGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAgcgg                                                                                     >  1:125723/1‑251 (MQ=255)
                                                                                                                                                                        gACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGgat                                                                           <  1:87003/251‑1 (MQ=255)
                                                                                                                                                                                      tCTGGAGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGAtt                                                              >  2:192796/1‑250 (MQ=255)
                                                                                                                                                                                         ggCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATTAACGGGATTaaa                                                           <  1:386114/250‑1 (MQ=255)
                                                                                                                                                                                            gtgtCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGc                                                       <  1:116196/251‑1 (MQ=255)
                                                                                                                                                                                              gtCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGt                                                     >  2:318321/1‑251 (MQ=255)
                                                                                                                                                                                                 aGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCAt                                                  <  2:103137/251‑1 (MQ=255)
                                                                                                                                                                                                      tCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCg                                             >  1:177288/1‑251 (MQ=255)
                                                                                                                                                                                                       cAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCg                                             <  1:355924/250‑1 (MQ=255)
                                                                                                                                                                                                         aaGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGt                                                                                             >  1:137184/1‑200 (MQ=255)
                                                                                                                                                                                                         aaGCCGCAAAGTCAGTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCCCCCGCACATTACCTGCGCAGAGGAAGCAATCGGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGt                                                                                             <  2:137184/200‑1 (MQ=255)
                                                                                                                                                                                                            ccGCAAAGTCATTAGTGACAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCtt                                       >  1:128259/1‑251 (MQ=255)
                                                                                                                                                                                                            ccGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGGTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCCTTGCCCGGTGAGCGGGTTGTAGCGGCGCTGTGGATGATCAAAGGGATTAAATTGCGTCATGGTCGTTCCtt                                       >  2:330143/1‑251 (MQ=255)
                                                                                                                                                                                                                  agtaaTTAGTGAAAACGTAAGTTCCGGTGTAAGCGGGTTTTTATAGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGGTTGGCTGGCGGTTCAGGCGCCACCTGCCAGGGGCGCTTATCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCgg                                 <  2:78619/247‑1 (MQ=255)
                                                                                                                                                                                                                         aGTGAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCCCCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCGGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATGAAATTGCGTCATGGTCGTTCCTTAATCGGGATATcc                          <  2:174491/251‑1 (MQ=255)
                                                                                                                                                                                                                            gAAAACGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACCTTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTCACAACTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCCGAATCCATTGCCCGGTGAGCGGGGTGTAGCGGCGATGTGGATGATCAACGGGGTTAACTTGCGTCATGGTCGTTCCTTAATCGTGATATCCCtg                       >  2:175544/1‑251 (MQ=255)
                                                                                                                                                                                                                              aaacgtcggtccagaagtAATCGGGTTTTTATCGCCTGTAACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTTGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGCGCGGGTTGTAGCGGCGATGTGGGTGATCAACGGGATTAAGTTGCGACATGGTCGTTCATTAATCGGGCTATCCCtgtg                     >  2:44218/17‑251 (MQ=255)
                                                                                                                                                                                                                                 cGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCAt                                                                                                             <  2:382050/160‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTAAGTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCAt                                                                                                             >  1:382050/1‑160 (MQ=255)
                                                                                                                                                                                                                                      gTCCCGGTGTAATCGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGGCGTTCCTTAATCGGGATATCCCTGTGGATGGCGt             >  1:340939/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  atcgggTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgccagtg  >  2:261652/4‑251 (MQ=255)
                                                                                                                                                                                                                                                  atcgggTTTTTATCGCATGTCACCCGCACATTACCTGCGCAGAGGAAGCATTCTGGTTCTTGCGCCGGTAACACCTGTTTGTCTGCCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGATTCCATTGCCCGGGGACGGGGTTGTAGCGGCGTTGTGGATTATCAACGGGCTTAAATTGCGTCTGGGTCGTTCCTTATTCGGGATTTCCCtgtgtttggaggtactgcacgtg  >  2:165485/4‑232 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
GCACCCATGGGTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTG  >  NC_000913/790585‑791077

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