Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,139,511 Δ1 bp coding (198/663 nt) fsaB ← fructose‑6‑phosphate aldolase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,139,5110C.96.7% 258.6 / ‑1.5 60coding (198/663 nt)fsaBfructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (34/24);  total (35/25)
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

AACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGAT  >  NC_000913/4139271‑4139755
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
aaCCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTcgcg                                                                                                                                                                                                                                           <  2:21497/251‑1 (MQ=255)
 aCCGTACTAATGCCGTTTGGGCCCTGGGCATCTACGCGGTTAACATACGGAACAACGTATTTTGCCCGTGCCAGTGCGGCTCATAAACCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGATGCCAGACATTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCgg                                                                                                                                                                                                                                          >  2:352407/1‑251 (MQ=255)
  ccGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGc                                                                                                                                                                                                                                         <  1:90076/251‑1 (MQ=255)
       cGAATGACGTTTCGGCCCTGGGGCTCTACGAGGTTCACATACGGAGAAACGCATTTTGCCCCTGCCAATGAGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGCCCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCAGCGCGTCGCGGCTCATg                                                                                                                                                                                                                                    >  2:85678/1‑251 (MQ=255)
            gCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATACCCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCAGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCCCATGGTCTg                                                                                                                                                                                                                               >  2:61594/1‑251 (MQ=255)
            gCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCt                                                                                                                                                                                                                                <  1:338484/250‑1 (MQ=255)
              cgtttcCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAg                                                                                                                                                                                                                             <  1:199560/247‑1 (MQ=255)
                        cTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAAt                                                                                                                                                                                                                   >  2:19753/1‑251 (MQ=255)
                          gggCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAAGTGCTGCCATACCGTCGTAAGTCGCCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTGCGCTTCTTCCACCAGCCCTGCGCGTCGCGGCTCATGGTCTGATCAGACAGAAt                                                                                                                                                                                                                   <  2:475670/249‑1 (MQ=255)
                                 tACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCtcatc                                                                                                                                                                                                          >  1:210485/1‑251 (MQ=255)
                                 tACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCtcatc                                                                                                                                                                                                          >  2:153459/1‑251 (MQ=255)
                                                      aaCGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGAcg                                                                                                                                                                                     >  2:411897/1‑251 (MQ=255)
                                                           aTTTTGCCCCTGCCAGTGCGTCTAATAACCCTTGTGCGGCCCTATATACAGCAGTGCCAAGTGTAGTCATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCCGAAGTCACCGGGTTTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCGCGCGTCGCGGCTCATGGTCCGAGCGAACCGCAGGCCCTCATCTCCCATCGCGTTTTGCCGACGCGgcc                                                                                                                                                                                >  2:9488/1‑250 (MQ=255)
                                                                gCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTCATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCCCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACt                                                                                                                                                                           >  2:483542/1‑251 (MQ=255)
                                                                     tGCCAGTGCGGATAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTccc                                                                                                                                                                       >  1:444297/1‑250 (MQ=255)
                                                                             cGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCACGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGCAGTCCCCGGGATTTTCCCCCCAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCTCTGCGCGGCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACt                                                                                                                                                              >  1:210882/1‑251 (MQ=255)
                                                                              ggCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGATGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCTGCACTTCCCATATGGACTc                                                                                                                                                             >  2:526125/1‑251 (MQ=255)
                                                                                  aataaCCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTg                                                                                                                                                         >  1:377520/1‑251 (MQ=255)
                                                                                  aataaCCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGAGGCAGCACTTCCCATATGGACTCCTTg                                                                                                                                                         <  1:297275/251‑1 (MQ=255)
                                                                                          ttGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGc                                                                                                                                                 >  1:296133/1‑251 (MQ=255)
                                                                                            gtgCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCg                                                                                                                                                >  1:497529/1‑250 (MQ=255)
                                                                                                ggcgCTATATGCAGCAGTGCCAAGTGTAGTAATACCCTCTTGTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAAt                                                                                                                                           <  2:296133/251‑1 (MQ=255)
                                                                                                    cTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTc                                                                                                                                       >  1:478240/1‑251 (MQ=255)
                                                                                                      atataCAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCAGCACCAATCGCTGTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTc                                                                                                                                       <  2:497529/249‑1 (MQ=255)
                                                                                                             gcagTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGtt                                                                                                                              >  1:221402/1‑251 (MQ=255)
                                                                                                                gTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTc                                                                                                                           >  1:357304/1‑251 (MQ=255)
                                                                                                                  gCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCa                                                                                                                          <  1:93939/250‑1 (MQ=255)
                                                                                                                     aaGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCAcccc                                                                                                                      <  1:153459/251‑1 (MQ=255)
                                                                                                                          tagtaATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCa                                                                                                                  >  2:401169/1‑250 (MQ=255)
                                                                                                                                tACCCTCTTTTTTCAGTATTTTAATTGATGCTAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGATAGTTGTCACCCCGGCAATgggg                                                                                                            >  2:362321/1‑250 (MQ=255)
                                                                                                                                 aCCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATgg                                                                                                                                                                 >  1:140308/1‑196 (MQ=255)
                                                                                                                                 aCCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATgg                                                                                                                                                                 <  2:140308/196‑1 (MQ=255)
                                                                                                                                  cccTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGGCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTGCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTGGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAat                                                                                                         <  1:352407/251‑1 (MQ=255)
                                                                                                                                      ctTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTGCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAatat                                                                                                       <  2:42083/249‑1 (MQ=255)
                                                                                                                                             cAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCTGCAATGGGGAATATTCGTGCCAg                                                                                              <  1:401169/251‑1 (MQ=255)
                                                                                                                                             cAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCa                                                                                               <  2:120748/250‑1 (MQ=255)
                                                                                                                                                    ttAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGCCGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCg                                                                                       >  2:440146/1‑251 (MQ=255)
                                                                                                                                                     tAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTAAACACCACACCCGGAATAGCGTCGCGCCGGCCCTTCGCTTATTACACCATCACTGCGCGTCGCGGCTAATGGTCTGCGCAAACAGAATGCCCTCATCACCAATCGCTTTCTGCAGACGCGGCAGCACTACCCATATGGAATCCTTGCTGGCAGCGATAATGCTCGGGTTCGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCg                                                                                       >  1:470179/1‑250 (MQ=255)
                                                                                                                                                      aaTTGCTGCCAGACCGTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGa                                                                                      <  2:419879/250‑1 (MQ=255)
                                                                                                                                                         tgctgcCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACtt                                                                                   <  1:483542/250‑1 (MQ=255)
                                                                                                                                                              ccAGACCTTCGGAAGTCACCGGGATTTTCACCACAACACCAAGACTAGCGTCGTGCAGGCGATTCGCTTCTTTAACAATCCCTGAGCGTAGCGGATCAAGACATGAGCAAACCGAATGCCCTAATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCACGGCAATGGGGCCTATGCGTGCCAGACGTTCGACTTCTGCg                                                                              >  1:503528/1‑250 (MQ=255)
                                                                                                                                                                      gcGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCATCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGc                                                                      <  2:201030/249‑1 (MQ=255)
                                                                                                                                                                        ggAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                                                                   >  1:361799/1‑251 (MQ=255)
                                                                                                                                                                             tCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGAGCAGGAGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAACCAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATCTGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTcc                                                               >  1:333427/1‑250 (MQ=255)
                                                                                                                                                                               aCCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTcc                                                               >  1:180217/1‑248 (MQ=255)
                                                                                                                                                                                    gATTTTCACCACAATACCCGGAATAGCGGCGCGCAGGCTCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAg                                                       <  2:444297/251‑1 (MQ=255)
                                                                                                                                                                                        ttCACCACAATCCCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATTGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTGTTGCAGACGCGGCAGCACTTCCCATATGGACTCCGTGCTGGCAGCGATAGTGCTCGGGTTAGTTGTCTCCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACGTCTGCGACGTTAGCGGTGTCCAGACACAGTTc                                                    <  2:198389/250‑1 (MQ=255)
                                                                                                                                                                                        ttCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTcc                                                   >  1:347103/1‑251 (MQ=255)
                                                                                                                                                                                                      cGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                                                               >  1:363898/1‑209 (MQ=255)
                                                                                                                                                                                                      cGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTGCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCGCCCCGGCAATGGGGAATATGCGTTCCAGACGTTCGACTTCTGc                                                                               <  2:363898/209‑1 (MQ=255)
                                                                                                                                                                                                      cGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCCCATGGTCTGAGCAAACAGAATGCCCTCATCCCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATCTGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCttt                                     >  1:343832/1‑251 (MQ=255)
                                                                                                                                                                                                        gAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCtttt                                    >  1:373570/1‑250 (MQ=255)
                                                                                                                                                                                                               gTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCCTGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAaagttcc                            >  2:448140/1‑251 (MQ=255)
                                                                                                                                                                                                               gTCGCGCAGGCGCTTCACTTCTTCCACCATCCCTGCGCGTCTCGGCTCATGGTCTGAGCAAACATAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGTCAGCGATAATTCTCGGGTTAGTTGGCGCCCCGGCGATGGTTAATATGCGTGCTAGACGGGCGACTTCTGCGACGGCAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagttcc                            <  1:484348/251‑1 (MQ=255)
                                                                                                                                                                                                                  gcgcAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGCTACAGTTCCATGATGTTTCCTTTTAcagttccagt                         >  2:135124/1‑251 (MQ=255)
                                                                                                                                                                                                                         cgcttcgcttCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGtt                  <  2:478240/251‑1 (MQ=255)
                                                                                                                                                                                                                         cgcttcgcttCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGtt                  <  1:411897/251‑1 (MQ=255)
                                                                                                                                                                                                                               gcttcTTCCACAATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGACTGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTgcagcc            >  2:465011/1‑250 (MQ=255)
                                                                                                                                                                                                                                     tcCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCt        <  1:61594/249‑1 (MQ=255)
                                                                                                                                                                                                                                      ccaccaTCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGACTGCCATCATCACCAATAGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATACTGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGc     >  2:445212/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ccatcccTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGAt  >  1:449122/5‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
AACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGAT  >  NC_000913/4139271‑4139755

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: