Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,239,943 T→A K136* (AAA→TAA)  mglB ← D‑galactose/methyl‑galactoside ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,239,9430TA98.0% 167.2 / ‑3.7 51K136* (AAA→TAA) mglBD‑galactose/methyl‑galactoside ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (27/23);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

CTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCATTCTGAGAATCATTCATCAGCAGCTGAACATC  >  NC_000913/2239712‑2240183
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
cttTCGCCTGAGCGGTGTCCCACATTGCGATATATAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCCCGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTCCGAACTGAATCTGACCGTCTTTGTTCAGATCCCAAACCTGATTCGCCGCCCAGTGTTTAGCAATCAAATAGACTTGAATAAGGCAGGACTATTAGGAGTCAgagcgaaagt                                                                                                                                                                                                                                 >  1:235270/1‑240 (MQ=255)
     gCCTGGGCGGTGTCCCACATTGCGTTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAAGTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTCTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGc                                                                                                                                                                                                                          <  2:173508/251‑1 (MQ=255)
              gtggccctcaTTGCTGTATCTAACTGTGACTGTTCAGTTTTGATGCCTTTATCGTTCAAGTCTTTAATCCCGTAAGTGTTACGTGCTTCTTCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTTAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTGAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGtt                                                                                                                                                                                                                 <  2:238183/243‑2 (MQ=255)
                   ccACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTAt                                                                                                                                                                                                            >  1:263521/1‑251 (MQ=255)
                      caTTGCGGTATCTAACTGTAATCGTTCAGTTTTGAGGCCTTTATCGTTCAATTCTTGAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACgt                                                                                                                                                                                                                                 <  2:275793/227‑1 (MQ=255)
                      caTTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACgt                                                                                                                                                                                                                                 >  1:275793/1‑227 (MQ=255)
                         tGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAgc                                                                                                                                                                                                       <  1:253044/250‑1 (MQ=255)
                                     ctgtaactgtTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCtt                                                                                                                                                                                                  >  1:46947/1‑243 (MQ=255)
                                     ctgtaactggTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCtg                                                                                                                                                                                                  <  2:46947/243‑2 (MQ=255)
                                         aactgtTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGtt                                                                                                                                                                                       <  1:276502/250‑1 (MQ=255)
                                             gtTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCtt                                                                                                                                                                                    <  1:41880/249‑1 (MQ=255)
                                                  gTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTga                                                                                                                                                                              >  1:184879/1‑250 (MQ=255)
                                                       gATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGaaaa                                                                                                                                                                        <  1:305882/251‑1 (MQ=255)
                                                        aTGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGACCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCt                                                                                                                                                                                                   >  1:22590/1‑223 (MQ=255)
                                                        aTGCCTTTAGCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCATATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAACCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCt                                                                                                                                                                                                   <  2:22590/223‑1 (MQ=255)
                                                                  cGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACg                                                                                                                                                             >  1:292559/1‑251 (MQ=255)
                                                                      cAATTCTTTAATCACGTAATTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAAGCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGtttt                                                                                                                                                         <  1:186294/251‑1 (MQ=255)
                                                                              tAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCAcgc                                                                                                                                                 >  1:53202/1‑251 (MQ=255)
                                                                                 tCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGAAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGAGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGAAGGTTCATTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCt                                                                                                                                              >  2:255194/1‑251 (MQ=255)
                                                                                    cGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTGTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTGACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTtc                                                                                                                                           <  1:148382/251‑1 (MQ=255)
                                                                                     gTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTtct                                                                                                                                          >  1:215299/1‑251 (MQ=255)
                                                                                                        gCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGc                                                                                                                       >  1:94016/1‑251 (MQ=255)
                                                                                                             cGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCg                                                                                                                  <  1:285008/251‑1 (MQ=255)
                                                                                                                           cACCTGTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTTTCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGtt                                                                                                     <  2:2137/250‑1 (MQ=255)
                                                                                                                            aCCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTg                                                                                                    >  1:76807/1‑250 (MQ=255)
                                                                                                                             ccTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACACTGATTCGCCGCCCAGTGTTTAGCAATCACATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGAACGTTTTGCCCACGCGCTTTCTCCATAACCGTACCCGCAGCTGCCGGGTCAACCAGGTTgtt                                                                                                  >  2:4341/1‑249 (MQ=255)
                                                                                                                                 tCAGCAGTAAGAACAGAATCTGACCGTATTTGTTCAGATCCCAACCCTGAATAGCCGCCCAGTGTTTCGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTCTCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGAACCCGCAGATGCCGGGTCACCGAGGTTGATTGcc                                                                                              >  2:311233/1‑251 (MQ=255)
                                                                                                                                           gAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTc                                                                                     >  1:281504/1‑250 (MQ=255)
                                                                                                                                            atcTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCAAGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGCTGATTGCACGTGCCTTCAc                                                                                   >  2:256576/3‑251 (MQ=255)
                                                                                                                                            aaCTGAATCTGACCGTATTTGTTCAGATCCCAGCCCTGATTCGACGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCAc                                                                                   <  1:26272/251‑1 (MQ=255)
                                                                                                                                                gAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGGGCATTCACCCCt                                                                               >  1:259315/1‑251 (MQ=255)
                                                                                                                                                                gTTCAGATCCCAACCCTGATTCGCCGACCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGGCGTAGCTATCCAGCGCCTTACGAGTCGGTTCTTTGTTGAAGTAAACCACCGGCACGTTTTGCCCAAGCGCTTTCTCAATCACCGTACCCGCAACTGCCTGGACAACCCGGATGATTACCCGTGCCTTCACCCCTTACGCCAGcgatccgt                                                               >  2:136314/1‑244 (MQ=255)
                                                                                                                                                                gTTCAGATCCCAACCCTGATTCCCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCCCGCGCTTTCTCAATCACCGCACCCGCAGCTGCCGGATCAACCAGGTTGATTGCCAGTTACTTCACCCCTCTCTCCAGCAATACGt                                                               >  2:143829/1‑251 (MQ=255)
                                                                                                                                                                    aGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGAt                                                           <  1:306492/251‑1 (MQ=255)
                                                                                                                                                                        cccAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGTCTTTGTCTTAGCTATCCAGCGCCTTCCGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTTCCTTCACCCCTTTCGCCAGCAATACGTCGATCTg                                                        <  1:7716/250‑1 (MQ=255)
                                                                                                                                                                          cAACCCTGATTCGCCGCCCAGAGTTTAGCAATCAAAGCGCCTTGAATAATGCCGGACTCTTAGGAGTCCGCGCCAACGTCATAAGCTAAGTAGACGATATCCAACGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCCCGTTTTGCCCACGCGCTTTATCACACACCGTACACGCAGCTGCCGGGTAAACCAGGATGATTGCCAGGTCCTTCACCCATTTACCCAGCactcagtcagtctcacc                                                     >  2:205354/1‑235 (MQ=255)
                                                                                                                                                                             accTGATTCGCCGCCCAGTGTTTAGCAATCAAGTCGCCTTGAATAATGCCGGACTATTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCACACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGTTCAACAAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTAGATCTGATCggt                                                  >  2:258818/2‑249 (MQ=255)
                                                                                                                                                                                tGATTCGCCGCCCAGTGCTTAGCAATCAAATCGCCTTGTATAATGACGGCCTCTTAGGAGTCAGTGCCAACGAAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATCCGTCGATCTGATCGTTCTg                                               >  1:16363/1‑251 (MQ=255)
                                                                                                                                                                                 gaggcgccgcCCAGTGTTTAGCTATCAAATCGCCGTGAATAATTCCGGACTCTTAGGAGTCATTGCCAACGTAGTAGTCTTTGTCGTAGCTATCCAGCGCTTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCGCGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCGTCACCCCGTTCGCCAGCAATACGTCGTTCTGATCGTTCTg                                               <  2:179313/246‑1 (MQ=255)
                                                                                                                                                                                        cgctcaatggttAGCAATCAAATCGCCTGGAATAATGCCGTACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGGAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCCCCGGCACGTTTTGCCCACTCGCTTTCTCAACCTCCGTACCCGCAGCTGCCGGGTCAACAAGGTTGATTGCCAGTTCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGAGCGTTCTGCTTGTACt                                       <  2:184640/241‑1 (MQ=255)
                                                                                                                                                                                        cgccCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACt                                       <  1:168869/251‑1 (MQ=255)
                                                                                                                                                                                         cacggcccacttcccatctaaaccccTTGAATAATGCCGGACTCTTAGGAGACAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTg                                      <  2:314564/227‑1 (MQ=255)
                                                                                                                                                                                              gtgtTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCCATACGTCGATCTGATCGTTCTGCTTGGACTGGTCAt                                 >  1:47780/1‑251 (MQ=255)
                                                                                                                                                                                               tgtTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCAtt                                >  1:107914/1‑251 (MQ=255)
                                                                                                                                                                                                 tttAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTCATCGTTCTGCTTGGACTGGTCATTCt                              >  1:109123/1‑251 (MQ=255)
                                                                                                                                                                                                          cAAATCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCATTCTGAGAAtcat                     <  1:284943/251‑1 (MQ=255)
                                                                                                                                                                                                            aaTCGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTg                                      >  1:314564/1‑232 (MQ=255)
                                                                                                                                                                                                               cGCCTTGAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGCAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCATTCTGAGAATCATtcatc                >  1:163621/1‑251 (MQ=255)
                                                                                                                                                                                                                     gAATAATGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCATTCTGAGAATCATTCATcagcag           <  1:288522/250‑1 (MQ=255)
                                                                                                                                                                                                                         aatGCCGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACGGGTCATTCTGAGAATCATTCATCAGCAGCTGa       <  1:256576/250‑1 (MQ=255)
                                                                                                                                                                                                                             ccGGACTCTTAGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCATTCTGAGAATCATTCATCAGCAGCTGAACATc  >  1:166445/1‑251 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
CTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGCATCCGGATGGCCCGGTTCACCTTTCAGCAGTACGAACTGAATCTGACCGTCTTTGTTCAGATCCCAACCCTGATTCGCCGCCCAGTGTTTAGCAATCAAATCGCCTTGAATAATGCCGGACTCTTTGGAGTCAGTGCCAACGTAGTAGGCTTTGTCGTAGCTATCCAGCGCCTTACGAGACGGTTCTTTGTTGAAGAAAACCACCGGCACGTTTTGCCCACGCGCTTTCTCAATCACCGTACCCGCAGCTGCCGGGTCAACCAGGTTGATTGCCAGTGCCTTCACCCCTTTCGCCAGCAATACGTCGATCTGATCGTTCTGCTTGGACTGGTCATTCTGAGAATCATTCATCAGCAGCTGAACATC  >  NC_000913/2239712‑2240183

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