Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,113,831 C→T 6.3% R166H (CGT→CAT)  wcaN ← putative regulatory subunit for GalU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,113,8310CT 6.3% 205.3 / 4.5 64R166H (CGT→CAT) wcaNputative regulatory subunit for GalU
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (2/2);  total (30/34)
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.03e-01

GTAACTGTCGCCGGTCATCAGCATTGCATCAACGGATTGTTTTTTCGCCAGCTCGGCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAG  >  NC_000913/2113545‑2114119
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                 
gTAACTGTCGCCGGTCATCAGCATTGCATCAACGGATTGTTTTTTCGCCAGCTCGGCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGt                                                                                                                                                                                                                                                                                                 <  1:432127/289‑1 (MQ=255)
        cGCCGGTCATCAGCATTGCATCAACGGATTGTTTTTTCGCCAGCTCGGCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCa                                                                                                                                                                                                                                                                                          >  2:568352/1‑288 (MQ=255)
            ggTCATCAGCATTGCATCAACGGATTGTTTTTTCGCCAGCTCGGCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCAcc                                                                                                                                                                                                                                                                                     <  1:576841/289‑1 (MQ=255)
                          cATCAACGGATTGTTTTTTCGCCAGCTCGGCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGt                                                                                                                                                                                                                                                                       >  2:676955/1‑289 (MQ=255)
                           aTCAACGGATTGTTTTTTCGCCAGCTCGGCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGtt                                                                                                                                                                                                                                                                      <  1:568352/289‑1 (MQ=255)
                            tCAACGGATTGTTTTTTCGCCAGCCCGGCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCACCAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCCCTTTAGTCTGGATGACGGAGTATTAAGAGAGGTCACCCGGCATACGTTTCGCCAGCACCTGGCTCCGTCCCGtt                                                                                                                                                                                                                                                                      >  1:258490/1‑288 (MQ=255)
                                       ttttttCGCCAGCTCGGCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGt                                                                                                                                                                                                                                                         <  2:394981/289‑1 (MQ=255)
                                                    tCGGCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCGCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATgg                                                                                                                                                                                                                                               <  1:50487/286‑1 (MQ=255)
                                                       gCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACTGTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGgcagca                                                                                                                                                                                                                                          <  1:699440/288‑1 (MQ=255)
                                                           tAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGtt                                                                                                                                                                                                                                     <  2:241684/289‑1 (MQ=255)
                                                                 cagtcagCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGCCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTACGTCGGCATGACGGAGTATTCAGGAGGCCCCCCGGCAATCCTTTTTCCACTCCCCGGTCTCGGCCCGGTTCCGTGAAACGTTGCAAtacggccagcagtttcgtacct                                                                                                                                                                                                                               >  1:148938/1‑268 (MQ=255)
                                                                             tACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAg                                                                                                                                                                                                                            <  2:551792/280‑1 (MQ=255)
                                                                                 tCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGGCTGAGTCCAGCGTCTGCCGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTg                                                                                                                                                                                                               <  1:77869/289‑1 (MQ=255)
                                                                                     cATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCa                                                                                                                                                                                                            >  2:9085/1‑288 (MQ=255)
                                                                                                          tCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCa                                                                                                                                                                                      >  2:703473/1‑289 (MQ=255)
                                                                                                          tCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGCCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCGGGCATCGTCGCTCACAACGTCTGGCa                                                                                                                                                                                      >  1:28265/1‑289 (MQ=255)
                                                                                                             aGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGt                                                                                                                                                                                    >  2:151573/1‑288 (MQ=255)
                                                                                                                tCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCa                                                                                                                                                                                <  1:676955/289‑1 (MQ=255)
                                                                                                                    gCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGgcagc                                                                                                                                                                                                                                           <  2:448339/226‑1 (MQ=255)
                                                                                                                    gCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGgcagc                                                                                                                                                                                                                                           >  1:448339/1‑226 (MQ=255)
                                                                                                                    gCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCa                                                                                                                                                                                >  2:16775/1‑285 (MQ=255)
                                                                                                                    gCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCa                                                                                                                                                                                <  1:16775/285‑1 (MQ=255)
                                                                                                                    gCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGAc                                                                                                                                                                            >  1:64634/1‑289 (MQ=255)
                                                                                                                    gCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGAc                                                                                                                                                                            >  1:457232/1‑289 (MQ=255)
                                                                                                                               gCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCGGGCATCGTCGATGACAACGTCTGGCAGTACCAAGACAAATGGGTTGt                                                                                                                                                                 >  1:678399/1‑289 (MQ=255)
                                                                                                                                                 ccTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTcgc                                                                                                                                                >  2:658713/1‑288 (MQ=255)
                                                                                                                                                  cTACGGCCATGATGTCTGGGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCgcgc                                                                                                                                              <  1:1528/289‑1 (MQ=255)
                                                                                                                                                  cTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCgcgc                                                                                                                                              <  1:122345/289‑1 (MQ=255)
                                                                                                                                                            gatgTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCa                                                                                                                                                                                >  2:253844/1‑245 (MQ=255)
                                                                                                                                                            gatgTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCa                                                                                                                                                                                <  1:253844/245‑1 (MQ=255)
                                                                                                                                                               gTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGa                                                                                                                                  <  1:151573/288‑1 (MQ=255)
                                                                                                                                                               gTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAg                                                                                                                                 <  2:258490/289‑1 (MQ=255)
                                                                                                                                                                 cggggTCCACGGTGTGCGGCTGATCCGGGTTTTGGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGt                                                                                                                                <  1:13905/284‑1 (MQ=255)
                                                                                                                                                                  tGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCCCGACAAATCGGTTGTCACCAATGGCAGGTCGCGCCCACAAAATGGAGTg                                                                                                                               >  1:30240/1‑288 (MQ=255)
                                                                                                                                                                  tGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTgg                                                                                                                              >  2:560819/1‑289 (MQ=255)
                                                                                                                                                                      tCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGccta                                                                                                                          <  1:703473/289‑1 (MQ=255)
                                                                                                                                                                        cAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCCATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACATATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGcctaa                                                                                                                         >  2:39903/1‑288 (MQ=255)
                                                                                                                                                                                                      aaaTTGAACAATGCGGCTGACTGTACGCTACCGGTCAAACGGCTTGTTACGCAGGATGACGTAGTAGTCAAAGAGGTAAACTGGCATACGTTTTGACAGCACCTGTCTGCGGCCCGGTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCGGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCGGACGAACGTTCATa                                                                                          <  1:661825/289‑1 (MQ=255)
                                                                                                                                                                                                                 tGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATATGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGcac                                                                                                                                            >  2:438436/1‑228 (MQ=255)
                                                                                                                                                                                                                 tGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATATGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGcac                                                                                                                                            >  2:438428/1‑228 (MQ=255)
                                                                                                                                                                                                                 tGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATATGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGcac                                                                                                                                            <  1:438428/228‑1 (MQ=255)
                                                                                                                                                                                                                 tGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATATGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGcac                                                                                                                                            <  1:438436/228‑1 (MQ=255)
                                                                                                                                                                                                                 tGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGGAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCCATCACAACGTCTGGCAGTAGCACGACAAATGGGTTGTCACCAATGGCCGGTCGCCCACACAAAACGGAATAGCCTAAACCTACAGGTTCCACCTGACGCACGTTACTAATGGtacagc                                                                                >  1:635645/1‑283 (MQ=255)
                                                                                                                                                                                                                 tGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCAACAATGGCAGGTCGCGCACACCAAATGGAGTGGCCTGAACCTAACGGTTCGCACTGACGCACGGTCATAATTGTCActc                                                                                >  1:569457/1‑286 (MQ=255)
                                                                                                                                                                                                                           tttACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAg                                                                     <  2:457232/289‑1 (MQ=255)
                                                                                                                                                                                                                              aCCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAGATg                                                                  <  2:30240/289‑1 (MQ=255)
                                                                                                                                                                                                                                    aCGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTg                                                                                                                                                                  >  1:122823/1‑187 (MQ=255)
                                                                                                                                                                                                                                    aCGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTg                                                                                                                                                                  <  2:122823/187‑1 (MQ=255)
                                                                                                                                                                                                                                    aCGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAGATGGa                                                                >  2:347191/1‑285 (MQ=255)
                                                                                                                                                                                                                                    aCGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAGATGGa                                                                <  1:347191/285‑1 (MQ=255)
                                                                                                                                                                                                                                     cGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCTCTGACGCACGTTCATAATGGTACCGACCGGCGGACCGATGGAcgga                                                            >  1:253288/1‑285 (MQ=255)
                                                                                                                                                                                                                                     cGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTACCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAg                                                                     <  2:323874/279‑1 (MQ=255)
                                                                                                                                                                                                                                     cGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTACCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAg                                                                     >  1:323874/1‑279 (MQ=255)
                                                                                                                                                                                                                                              ggCTCGTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGcc                                                   <  2:678399/288‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGa                                                                                                                                  <  2:560455/197‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGa                                                                                                                                  >  1:560455/1‑197 (MQ=255)
                                                                                                                                                                                                                                                          tGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGt                                                                                               <  2:32119/232‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGt                                                                                               >  1:32119/1‑232 (MQ=255)
                                                                                                                                                                                                                                                            gATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATaa                                                                                         >  2:261520/1‑236 (MQ=255)
                                                                                                                                                                                                                                                            gATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATaa                                                                                         <  1:261520/236‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ggTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGGCCGGCGGa                                                                        >  1:325263/1‑233 (MQ=255)
                                                                                                                                                                                                                                                                                ggTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGGCCGGCGGa                                                                        <  2:325263/233‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      ccGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTaa          <  1:658713/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           aTACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCAt     <  2:182905/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              cGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAag  <  2:148938/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              cGTTTGGCCAGCACCTGGCGGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCAcgc                                <  2:255287/259‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                 
GTAACTGTCGCCGGTCATCAGCATTGCATCAACGGATTGTTTTTTCGCCAGCTCGGCAATAGCATCAGTCAGCTGAATACGTCCCCATGCACCAGGCTGAGTACGTTCCAGTTCCGGCCAAATATCGGCAGAAAGCACATAGCGACCTACGGCCATGATGTCTGAGTCCAGCGTCTGCGGCTGATCCGGTTTTTCGATAAATTCAACAATGCGGCTGACTTTACCCTCACGGTCCAGCGGCTCTTTAGTCTGGATGACGGAGTATTCAGAGAGGTCACCCGGCATACGTTTTGCCAGCACCTGGCTGCGGCCCGTTTCGTTGAAACGTGCAATCATGGCAGCAAGGTTGTAACGTAGCGGGTCGGCGCTGGCATCGTCGATCACAACGTCTGGCAGTACCACGACAAATGGGTTGTCACCAATGGCAGGTCGCGCACACAAAATGGAGTGGCCTAAACCTAAAGGTTCGCCCTGACGCACGTTCATAATGGTCACGCCCGGCGGACAGATGGACTGTACTTCCGCCAGCAGTTGACGCTTCACGCGCTGCTCAAGGAGTGATTCTAACTCATAAG  >  NC_000913/2113545‑2114119

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