Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,771,979 G→T intergenic (‑235/‑302) yibI ← / → mtlA DUF3302 domain‑containing protein YibI/mannitol‑specific PTS enzyme IICBA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,771,9790GT92.7% 197.0 / 5.5 68intergenic (‑235/‑302)yibI/mtlADUF3302 domain‑containing protein YibI/mannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base T (32/31);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AAACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACAT  >  NC_000913/3771739‑3772221
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
aaccATATCTGCCGCGGATTCACGTTTATCAGTGTTTTTTTTGGGTTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTTTAGATCACATCCGGTGAATGCCGTAATAAAAGATCTTGTGATTCAGATCACAAATATTCAACaaa                                                                                                                                                                                                                                          <  1:544483/248‑1 (MQ=255)
 aaCAGATCTGCCTCGGAGTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAAc                                                                                                                                                                                                                                         <  1:506528/251‑1 (MQ=255)
   cataccttccttggATTCACGTATCACATTGTTGTTGTTTGGTTGCCAGCCAGAAGGGAGTTAGGTGGATATTTTGTCAATAATCCGGGTTCCCGATTCTCGCCATAACACCAAAGAATACTTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATTACGAAGGCATAACATGCTGTAGATCACAGCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCa                                                                                                                                                                                                                                       <  2:245433/239‑1 (MQ=255)
         tGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATcaaa                                                                                                                                                                                                                                  >  1:12987/1‑250 (MQ=255)
                       gTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGacac                                                                                                                                                                                                                   >  1:787748/1‑251 (MQ=255)
                       gTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGaca                                                                                                                                                                                                                    <  1:84703/250‑1 (MQ=255)
                       gTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGaca                                                                                                                                                                                                                    >  2:84703/1‑250 (MQ=255)
                                    gttTTTGGGTTGGCAGCCAGAATGGATTCAGGCTGAAATTTTGACAATAATCCGGGTTCGCGATGCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAATCAAAAATGTGACACTACTCACATTTaa                                                                                                                                                                                                      <  1:821851/251‑1 (MQ=255)
                                    gttTTTGGGCTGGCAGACAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTa                                                                                                                                                                                                       <  1:181519/250‑1 (MQ=255)
                                     tttttGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTaa                                                                                                                                                                                                      >  1:816864/1‑250 (MQ=255)
                                          gggCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGAGACACTACTCACATTTAAATGcc                                                                                                                                                                                                 >  1:347345/1‑250 (MQ=255)
                                           ggCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCTATTCTCGCCATAACACCAAAGAATAAGTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCAt                                                                                                                                                                                               <  1:403023/251‑1 (MQ=255)
                                               ggCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCAtttt                                                                                                                                                                                            <  1:640181/250‑1 (MQ=255)
                                                  aGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATAtt                                                                                                                                                                                                                                                 <  1:855067/194‑1 (MQ=255)
                                                     cAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCCTAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAAAGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACAATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGa                                                                                                                                                                                      >  2:171704/1‑250 (MQ=255)
                                                     aaGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGaa                                                                                                                                                                                     <  2:872464/250‑1 (MQ=255)
                                                      aGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAATAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGaaa                                                                                                                                                                                    <  1:393563/251‑1 (MQ=255)
                                                       gAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACAtt                                                                                                                                                                                                         >  1:848775/1‑229 (MQ=255)
                                                       gAAGGGAGTCAGGCTGATATTTTGAAAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACAtt                                                                                                                                                                                                         <  2:848775/229‑1 (MQ=255)
                                                        aaGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAAt                                                                                                                                                                                  <  1:258186/251‑1 (MQ=255)
                                                            gAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATcgcc                                                                                                                                                                              <  1:793036/251‑1 (MQ=255)
                                                                    cTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGCACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCttgtt                                                                                                                                                                      >  2:537163/1‑251 (MQ=255)
                                                                     tGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAAAATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTATGATTCAGATCACAAATATTCAACAAACCATCAAAACCAAAATGTGACACTACTCACATTTAAATGCCATTTATCACGAAACTCACCGCCttgttg                                                                                                                                                                     >  2:550502/1‑251 (MQ=255)
                                                                           tttGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCttttt                                                                                                                                                               >  1:790400/1‑251 (MQ=255)
                                                                             tGACAATAATCCGGGTTCACGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCtttttt                                                                                                                                                              <  1:341367/250‑1 (MQ=255)
                                                                              gACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTac                                                                                                                                                            >  1:861706/1‑251 (MQ=255)
                                                                              gACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTac                                                                                                                                                            >  1:591695/1‑251 (MQ=255)
                                                                                  agaatCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACaa                                                                                                                                                        <  1:119473/249‑1 (MQ=255)
                                                                                    aatCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGc                                                                                                                                                      >  1:159522/1‑251 (MQ=255)
                                                                                     atCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTAACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCg                                                                                                                                                     >  1:276502/1‑251 (MQ=255)
                                                                                              cgcgATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGtga                                                                                                                                             <  1:574122/250‑1 (MQ=255)
                                                                                              cgcgATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGATTTTTTACACAAGCGTTTtgtaat                                                                                                                                            >  2:186955/1‑248 (MQ=255)
                                                                                                  atgctcGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTAGGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAac                                                                                                                                        <  1:851337/248‑1 (MQ=255)
                                                                                                  aTTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAac                                                                                                                                        >  1:629074/1‑251 (MQ=255)
                                                                                                             aaCACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGtttt                                                                                                                                                 >  1:414675/1‑231 (MQ=255)
                                                                                                             aaCACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCAGCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGtttt                                                                                                                                                 <  2:414675/231‑1 (MQ=255)
                                                                                                               cacCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAAtt                                                                                                                            >  1:762299/1‑250 (MQ=255)
                                                                                                                    aaGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGACCGTCACGTCAATTACctct                                                                                                                      >  1:384762/1‑251 (MQ=255)
                                                                                                                     aGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCTCAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACGCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACctctc                                                                                                                     <  1:76626/251‑1 (MQ=255)
                                                                                                                       aataatTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGAAGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACctctct                                                                                                                    >  1:309275/1‑250 (MQ=255)
                                                                                                                             ttttAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTAccccc                                                                                                              <  1:24054/250‑1 (MQ=255)
                                                                                                                             ttttAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGAGTTTTGTGATGAACGTCACGTCAATTACCTCTCTAccccc                                                                                                              <  1:182982/250‑1 (MQ=255)
                                                                                                                              tttAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCt                                                                                                             >  1:224505/1‑250 (MQ=255)
                                                                                                                                 agagGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCtatag                                                                                                         <  1:604473/251‑2 (MQ=255)
                                                                                                                                                   aGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTAACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTAGCCCCTATATTTATGTGATTGATATcac                                                                                       >  2:579229/1‑251 (MQ=255)
                                                                                                                                                       aCATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGACGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATcacaca                                                                                    >  2:89193/1‑250 (MQ=255)
                                                                                                                                                           aaCGATTGTATGACGAAGACATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAgg                                                                               >  1:278313/1‑251 (MQ=255)
                                                                                                                                                               aTTGTATGACGAAGGCATAACATTCTGTAGATCACATCAGGTGAACGACGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAACTCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATCTTTATGTGATTGATATCACACCAAAGGCCGt                                                                           >  2:4396/1‑251 (MQ=255)
                                                                                                                                                                      gatgAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCAGTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTgg                                                                    <  1:116126/248‑1 (MQ=255)
                                                                                                                                                                              cATACCATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTAATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTa                                                             >  1:566601/1‑250 (MQ=255)
                                                                                                                                                                              aaTAACATGCTGTAGATCACATCAGGTGAACTCCGTCAGAAAATATCTGGAGATTCAGATCACAAATAATCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCAGTTTTAGCGAAAATCGCCGCCTTGTTGCTGTTTTACACAAGCATTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATGTATGTGATTGATATCCCACAAAATGCCGTCGACTGGACAGTTaa                                                            <  2:316116/250‑1 (MQ=255)
                                                                                                                                                                               aTAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAAc                                                           <  1:193869/251‑1 (MQ=255)
                                                                                                                                                                                aaaCATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAAcc                                                          <  1:550219/250‑1 (MQ=255)
                                                                                                                                                                                 aaCATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTTACCCTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGAGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAAcc                                                          <  2:187689/250‑1 (MQ=255)
                                                                                                                                                                                  aCATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCg                                                         >  1:560712/1‑250 (MQ=255)
                                                                                                                                                                                             aTCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTATCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGTCCGTCGACTGGACAGTTAACCGATTCAGTGcc                                               >  1:473926/1‑249 (MQ=255)
                                                                                                                                                                                             aTCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAg                                             <  1:494723/251‑1 (MQ=255)
                                                                                                                                                                                                acaTCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGAtt                                          >  1:75571/1‑251 (MQ=255)
                                                                                                                                                                                                            cGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAACAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATcc                               >  2:324208/1‑249 (MQ=255)
                                                                                                                                                                                                             gCCGGAAGAAAATATCTTGTTATTCAGATCCCAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCCCATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCTAGATTTCGCAGTATCTAc                             <  1:672258/251‑1 (MQ=255)
                                                                                                                                                                                                                gTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACa                            >  1:221473/1‑249 (MQ=255)
                                                                                                                                                                                                                      aaaTATCTTGCGATTCAGATCACAAATATTCAACAAACCCTCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCgg                    >  2:800024/1‑251 (MQ=255)
                                                                                                                                                                                                                              tgtgATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGTTAcc               <  1:89193/248‑1 (MQ=255)
                                                                                                                                                                                                                                tgATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTgc          >  2:864239/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTTACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTgc          <  1:550502/250‑1 (MQ=255)
                                                                                                                                                                                                                                    tCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCTATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGACCGTCGACTGGCCAGTTACCCGATTCAGTGCCAGATTTCGCCATATATACAAGGTCCGGCTACCTCTgccgcc      >  2:778383/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTgccgcc      >  1:57561/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aTCACAAAGATTCAACAAACCATCAAAAAAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCTTATATTTATGTGATTGATATCACACAAAAGGCCGTCGCCTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCggaaacatctgccaacacat  >  2:445756/1‑233 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
AAACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACAT  >  NC_000913/3771739‑3772221

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