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.98.0% 204.9 / ‑2.6 49coding (198/663 nt)fsaBfructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (21/27);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCG  >  NC_000913/4139267‑4139753
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ccTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGt                                                                                                                                                                                                                                             >  2:429589/1‑251 (MQ=255)
             aaTGCCGTTTCCCCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCCTCGCGGCTCATGGt                                                                                                                                                                                                                                <  2:314174/251‑1 (MQ=255)
              aTGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTc                                                                                                                                                                                                                               <  2:106014/251‑1 (MQ=255)
                 ccGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCATGCTCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGa                                                                                                                                                                                                                            <  1:444694/251‑1 (MQ=255)
                         gCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGACCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCAGCGCGTCGCGGCTCATGGTCTGAGCAAACAg                                                                                                                                                                                                                    <  2:364504/251‑1 (MQ=255)
                                      aCGCGGTTCACATACGGAGCAACGAATTTTAACCCTGCCAGTGCGGCTCATACCCCTTGAGCGGCGCTATATACACCAGTACCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACATTCGGAAGTCACCGGGATTTTCACCCCAATACCCGGAATAACGTCGCGCCGGCGCTTCGCTTCTTCCACAATCCCTGCGAGTCGCGGCCCATGGTCTGAGCAAACCGAATGCCCtcatcc                                                                                                                                                                                                       >  2:297189/1‑250 (MQ=255)
                                        gcgGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCAcc                                                                                                                                                                                                     <  1:296580/251‑1 (MQ=255)
                                             tAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATc                                                                                                                                                                                                 <  2:152672/250‑1 (MQ=255)
                                                                ttttGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGACGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGgc                                                                                                                                                                               >  2:5689/1‑249 (MQ=255)
                                                                                gCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCACGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGAc                                                                                                                                                             >  1:319424/1‑251 (MQ=255)
                                                                                gCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGa                                                                                                                                                              <  2:72066/250‑1 (MQ=255)
                                                                                    cTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCt                                                                                                                                                         >  1:246336/1‑251 (MQ=255)
                                                                                          aCCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCTCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTGGTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTgg                                                                                                                                                   <  1:322122/251‑1 (MQ=255)
                                                                                           cccTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGTCAGCACt                                                                                                                                                                         <  2:64461/228‑1 (MQ=255)
                                                                                           cccTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACt                                                                                                                                                                         >  1:64461/1‑228 (MQ=255)
                                                                                            ccTTGTGCGGCGCTAGATACAGCAGTGCCAAGTGTAGTAATACCCTCGTTTTTCAGTATTTTAATTGCTGCCAGACCGTCGGAAGTCACCGGGATTTTCACCACAATTCCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGGGGCTCATGGTCTGAGCAAACAGAATGCCCTCAGCACCAATCGCTTTTTGCAGTCGCGGCAGCACTTCCCATAGAGACACCTGGCTGGCa                                                                                                                                                 <  2:38746/251‑1 (MQ=255)
                                                                                                  gcggcgCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGAt                                                                                                                                            <  1:151332/250‑1 (MQ=255)
                                                                                                    ggcgCTATATACAGCAGTTCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCTCCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGTCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAAt                                                                                                                                         <  1:306162/251‑1 (MQ=255)
                                                                                                    ggcgCTATATAAAGCAGTGCCAAGTGTAGTAATACCCTATTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCAAAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCACTGCGCGTCGCGGCTCATGGTCTGAGCCAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAAt                                                                                                                                         >  1:77701/1‑251 (MQ=255)
                                                                                                          atataCAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCgg                                                                                                                                   <  2:385087/251‑1 (MQ=255)
                                                                                                                  cagTGCCAAGTGTTGTACTACCCTCTTTTTTCAGTATTTTCTTTGCTGCCAGACCTTCGGAAGTCCCCGGGCTTTTCACACCAATACCCGGACTAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCCTCCCTGCGCGTCGCGGCTAATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGCTAATGCTCGGGTTAGTTg                                                                                                                           >  2:468042/1‑251 (MQ=255)
                                                                                                                        cAAGTGTAGTAATACCCTCTTTTTTCAGTCTTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGACTGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCACATATGGACTCCTTGCTGGCAGAGATAATGCTCGGTTTAGTTGTCAccc                                                                                                                     >  2:214602/1‑251 (MQ=255)
                                                                                                                                  aattacctcctgtgtcactatcttAATTGCTGCCAGGCCTGCGGAAGTCACCGGGGTTTTCTCCACAATACCCGGCATAGCGTCGCGCAGGCGCTTCGCGTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACGTCCCATATGGACTCCTTGCTGGCAGCTATACGTCTCGGGTTAATTGTCGCCCCGGCAATgg                                                                                                            <  2:210251/228‑1 (MQ=255)
                                                                                                                                  aaTACCCTATTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATggg                                                                                                           >  2:442973/1‑251 (MQ=255)
                                                                                                                                            tttttCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGt                                                                                                 >  1:468109/1‑251 (MQ=255)
                                                                                                                                                tCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGTTTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCACACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGcc                                                                                              >  1:65905/1‑250 (MQ=255)
                                                                                                                                                         tAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCg                                                                                     >  1:452830/1‑250 (MQ=255)
                                                                                                                                                              gctgcCAGACCTTAGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCt                                                                               >  1:342490/1‑251 (MQ=255)
                                                                                                                                                               ctgcCAGACCTTCGGCAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTg                                                                              >  1:106730/1‑251 (MQ=255)
                                                                                                                                                                   cagaccttcgcccgaccccGGGCTTTTAACCACAATACCAGGCATAGCGTCGCGCAGGCGCTTCCCTTCTACCACCATCCCTGCCCGCCGCGGCTCATGGACTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCCTATGGACTCCTTGCTGGCCGCGATAATGCTAGGGTTTGTTGTCACCCCGGCCATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGAc                                                                          >  2:295972/18‑251 (MQ=255)
                                                                                                                                                                    agcccTTCGGAAGTCACCGGGATTTTCGCCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCCCCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCGTCTCCAATCGCGTTTTGCAGACGCGGCAGCACTTCCCATATGGACCCCTTGCTGGCAGTGATAATGCTCGGGTTAGTGTTCACCCCGGCAATGTGGAATATGCGTGCTAGACGTTCGGCTTCTGCGACg                                                                         <  2:319424/248‑1 (MQ=255)
                                                                                                                                                                     gACCTTCGGAAGTCACCGGGATTTTCACCACAATACCAGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGt                                                                        >  1:307632/1‑251 (MQ=255)
                                                                                                                                                                        cTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCGGGCGCTTAGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCCGCACTTCCCATATGGACTCCTTGCTGGCAGCGCTAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGtt                                                                       >  2:187604/1‑249 (MQ=255)
                                                                                                                                                                          tCGGAAGTTACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAAGCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCCATAATGCTCGGGTTAGTTGTCACCCCTGCAAGGGGTAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGc                                                                    <  1:318129/250‑1 (MQ=255)
                                                                                                                                                                            ggAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGTGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                                                                 <  1:227633/251‑1 (MQ=255)
                                                                                                                                                                               aGTCACCGGGATTTTCACCACACTACACGGCATCGCGAAGCGCAAGCGCTTCGCTTCTCCCACCATCCATGCGCGTCGCGGCTCCTGGTCTGCGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTc                                                              >  2:437664/1‑251 (MQ=255)
                                                                                                                                                                                gTCACCGAGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTTCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGGTGTCACCCCGGCAATTGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCTGTGTcc                                                             <  1:231158/251‑1 (MQ=255)
                                                                                                                                                                                  cACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATACTGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTcc                                                             >  1:271542/1‑249 (MQ=255)
                                                                                                                                                                                            ttCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCATCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGGTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGGTAGCGGTGTCCAGATACAGTTcc                                                 >  2:10211/1‑251 (MQ=255)
                                                                                                                                                                                                  acaATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTGGCTGGCAGCGATAATGCTCGGGTGAGTTGTCACCCCGGCAATGTGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatgat                                            <  1:225656/251‑1 (MQ=255)
                                                                                                                                                                                                         ccGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCGTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTGCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCtt                                    <  2:307632/251‑1 (MQ=255)
                                                                                                                                                                                                                aGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTGTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCGTTTAcagt                             <  2:395269/251‑1 (MQ=255)
                                                                                                                                                                                                                aGCGTCGCGCAGGCGCTTCGCGTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCACTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCAGATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcag                              <  2:248880/250‑1 (MQ=255)
                                                                                                                                                                                                                            gcgcTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGt                 <  2:246336/251‑1 (MQ=255)
                                                                                                                                                                                                                             cgcttcgcttCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTGCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTACCATGATGTTTCCTTTTACAGTTCCAGTTCATGtt                <  2:92444/251‑1 (MQ=255)
                                                                                                                                                                                                                             cgcttcgcttCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTGTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCGCCCCGGCAATGGGGAATATGCGTGCCGGACGTTCTACTTCTGCGACGTTAGCGTTTTCCAGATACAGGTCCATGAGGTTTCCGTTTACAGTTCCAGTTCATGtt                <  2:216243/251‑1 (MQ=255)
                                                                                                                                                                                                                             cgcttcgcttCTTCCACCATCCCTGCGCGGCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGtt                <  1:466414/251‑1 (MQ=255)
                                                                                                                                                                                                                                 tcgcttcTTCCACCATCCCTGCGCGTCGCCGCTCATGGTCTGAGCAAACATAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGTCTTCTGCGACGTTAGCGGTGTCCTGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTgc              <  2:139772/249‑1 (MQ=255)
                                                                                                                                                                                                                                           caccaTCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACGTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCGGCGACGTTAGCGGTGTCCAGATACAGTTCCAGTATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCg  <  1:468042/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
CCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCG  >  NC_000913/4139267‑4139753

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