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.97.7% 185.1 / ‑3.7 43coding (198/663 nt)fsaBfructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (21/21);  total (22/21)
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.96e-01

GAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGAT  >  NC_000913/4139270‑4139755
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
gAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACCGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTcgc                                                                                                                                                                                                                                            >  2:123098/1‑251 (MQ=255)
gAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTcgc                                                                                                                                                                                                                                            >  2:324780/1‑251 (MQ=255)
         gAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAATATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTCTATACAGAAGTGCCAAGTGTAGAAATACCCTCTTTTTTCAGTATTTTAATTGCTGCTAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCCGCTCCTgg                                                                                                                                                                                                                                   >  2:113773/1‑251 (MQ=255)
               cgtttcCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTg                                                                                                                                                                                                                               >  1:149327/5‑249 (MQ=255)
                     cGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGGGCCAAGTGAAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGCGCAAAc                                                                                                                                                                                                                        >  1:382994/1‑250 (MQ=255)
                         cTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAAt                                                                                                                                                                                                                   >  1:124956/1‑251 (MQ=255)
                              cATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCt                                                                                                                                                                                                              <  1:68723/251‑1 (MQ=255)
                                   gcgcgGTTAACTTTCGGAGCAACGTATTTGGCCCCTGCCAGATCGGCTAATATCCCTTGTGCTGCGCGAGATACAGCAGTGCCAAGTGGAGTAATAGCCGCTTTGTTCAGTATTGTAATCGCTGCCAGACCTGCGGAAGTCACTTGGATTTTCACCACAATACCCGGAATAGCGTCGCGGAGGCGCTTCGCGTCTTCCACCATCCC‑TGCGCGTCGCGGCGCATGATCTGAGCAAACATAATGCCCtcatc                                                                                                                                                                                                          <  2:137500/249‑1 (MQ=255)
                                    cgcgGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCAc                                                                                                                                                                                                        >  1:186990/1‑251 (MQ=255)
                                        gTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGCGCAAACAGAATGCCCTCATCACCAAt                                                                                                                                                                                                    >  2:431610/1‑251 (MQ=255)
                                        gTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAAt                                                                                                                                                                                                    >  1:300304/1‑251 (MQ=255)
                                          tAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATc                                                                                                                                                                                                   <  1:440064/250‑1 (MQ=255)
                                             cATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCtt                                                                                                                                                                                               <  2:271232/251‑1 (MQ=255)
                                              aTACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCttt                                                                                                                                                                                              <  2:7280/251‑1 (MQ=255)
                                                     gCAACGTATTTTGCCCCTGCCAGTGCGGCTCATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTGTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATGCCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGa                                                                                                                                                                                       <  1:288617/251‑1 (MQ=255)
                                                           tATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGgc                                                                                                                                                                                 >  1:281425/1‑251 (MQ=255)
                                                                     cTGCCAGTGAGGCTAATTACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCCCTTccc                                                                                                                                                                       >  1:81053/1‑251 (MQ=255)
                                                                       tccAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTccc                                                                                                                                                                       <  1:256696/248‑1 (MQ=255)
                                                                                               gcggcgCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACCGCATGCCCTCATCACCAATCGCTTTTTGCAGCCGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGAt                                                                                                                                              >  2:354425/1‑250 (MQ=255)
                                                                                                cggcgCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATaa                                                                                                                                            <  1:73125/251‑1 (MQ=255)
                                                                                                 ggcgCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAAt                                                                                                                                           >  1:55782/1‑251 (MQ=255)
                                                                                                  gcgcTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGAATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATCTGGACTCCTTGCTGGCAGCGATAATg                                                                                                                                          >  2:358066/1‑251 (MQ=255)
                                                                                                       atataCAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTc                                                                                                                                       <  1:230533/249‑1 (MQ=255)
                                                                                                            cagcagTGCCAAGTGTAGTTATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAg                                                                                                                                <  2:186990/251‑1 (MQ=255)
                                                                                                                     cAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGGCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCGCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCAGATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCAccc                                                                                                                       <  2:281425/251‑1 (MQ=255)
                                                                                                                        gtgtAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCgg                                                                                                                    >  1:257924/1‑251 (MQ=255)
                                                                                                                           tagtaATACCCTCTTTTTTCAGGATTTTAATTGCTGCCAGACCTTCGGAAGTCACCTGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCGTCTTCCACCATCCC‑TGCGCGTCGCAGCTCATGGTCTGCGCAAACAGAATGCCCTCATCCCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGACCGTGTTAGCTGTCACACCGGCat                                                                                                                 >  2:259728/1‑250 (MQ=255)
                                                                                                                                  tcccTCTTTTTTCAGGATTTTAATTGCTGCCAGACCTTCGGAAGTCACCTGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGa                                                                                                           <  2:419702/249‑1 (MQ=255)
                                                                                                                                    ccTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAata                                                                                                        <  1:55909/251‑1 (MQ=255)
                                                                                                                                              cAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGcc                                                                                                >  1:342501/1‑249 (MQ=255)
                                                                                                                                                     ttAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTc                                                                                        <  1:447953/250‑1 (MQ=255)
                                                                                                                                                       aaTTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGAc                                                                                     >  2:439456/1‑251 (MQ=255)
                                                                                                                                                       aaTTGCTGCCAGACCTTCGGAAGGCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGGTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGAc                                                                                     <  1:302984/251‑1 (MQ=255)
                                                                                                                                                              gcCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGCAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCC‑GCCCGTCGAGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCACTTTTTGCCGACGCGGCAGCACTTCCCATATGGACTCATTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCGCTGGGGAATATGTGTGCCAGACGTTCGAATTCTGCg                                                                              >  1:331736/1‑251 (MQ=255)
                                                                                                                                                                cAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGAc                                                                            <  1:123098/251‑1 (MQ=255)
                                                                                                                                                                                aCCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTcc                                                               >  1:66536/1‑248 (MQ=255)
                                                                                                                                                                                 ccGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTcc                                                               >  1:172594/1‑247 (MQ=255)
                                                                                                                                                                                     gATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAg                                                       <  1:222045/251‑1 (MQ=255)
                                                                                                                                                                                        tttCACCACAATACCCGGAATAGCGGCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTTTCCAGATACAGtt                                                     <  2:149327/250‑1 (MQ=255)
                                                                                                                                                                                               acaATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCATACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatgatg                                             <  1:235389/251‑1 (MQ=255)
                                                                                                                                                                                                    aCCCGGAAAAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTTGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCGTTTTGCAGACGCGGCAGCACTTCCCATAAGGACTCCTTGCTTGCCGCGATAATGCTCGGGTTATTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGGCTTCTGCGACGTTAGCGGTGGCCAGATACAGTTCCATGATGTTTcc                                        <  2:244403/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cttcCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCt        >  1:22769/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tcCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCt        <  2:335579/249‑1 (MQ=255)
                                                                                                                                                                                                                                            catcccTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGAt  >  2:318214/4‑250 (MQ=255)
                                                                                                                                                                                                                                             atcacTGCGCGTCGCGGCTCATGGTCTGAGCCAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCCCTTCCCATATGGACTCCTTGCTGGCAGCGATCATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTAGACTTATGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTGTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGAt  >  2:60334/5‑249 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
GAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGAT  >  NC_000913/4139270‑4139755

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