Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 860,260 G→A R783C (CGT→TGT)  ybiW ← putative pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913860,2600GA98.4% 220.6 / NA 64R783C (CGT→TGT) ybiWputative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (36/27);  minor base C (1/0);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTACCCACTGAGCCAATGACCGCAGTAG  >  NC_000913/860021‑860502
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
aCGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCAGCGACACACACTACCAGa                                                                                                                                                                                                                                          >  1:23567/1‑250 (MQ=255)
      aTAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGAt                                                                                                                                                                                                                                   <  2:229389/251‑1 (MQ=255)
                 ccTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATCTTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCggg                                                                                                                                                                                                                         >  1:213757/1‑250 (MQ=255)
                 ccTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCggg                                                                                                                                                                                                                         >  1:232682/1‑250 (MQ=255)
                         gAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAAGTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGtttttt                                                                                                                                                                                                                <  1:110151/251‑1 (MQ=255)
                           aaaTGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCg                                                                                                                                                                                                              <  2:2415/251‑1 (MQ=255)
                           aaaTGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCg                                                                                                                                                                                                              <  1:430073/251‑1 (MQ=255)
                            aaTGAAACCGTACCAGTTTACTAACCATCTGACGCAGAGAACAACTTATCTTTCGAATGCAAGTAAGGATGACATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGCTATTCCAGCTTATGTTCAGTACAGGCGATGATATCTTCCTGAGCGTCTGTAGCGAGCGCTGGGACGAACCCGGACTAGCCCGCGCCCCCGACTACCAGATCGCGATACTGATCGGGCTGttttttagc                                                                                                                                                                                                             >  2:444275/1‑248 (MQ=255)
                               gAAACCGTAACAGTGCAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAAGAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCAt                                                                                                                                                                                                           <  1:426518/250‑1 (MQ=255)
                                aaaCCGTAACAGTGTAATATCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATcc                                                                                                                                                                                                         <  2:213757/251‑1 (MQ=255)
                                    cGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGCAAGGATGAAATTGATGATGTGAATGATTTAGCCAGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCcagc                                                                                                                                                                                                      >  2:147873/1‑250 (MQ=255)
                                            tgtAATAACAATGTGACGCAGCGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGCTGATATCGTCCTGAGCGTCTGGCGAGAGCGCGGTGCAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTcg                                                                                                                                                                                             >  2:208001/1‑251 (MQ=255)
                                             gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTcg                                                                                                                                                                                             <  1:5767/250‑1 (MQ=255)
                                                     aaTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGaaa                                                                                                                                                                                      >  1:284007/1‑249 (MQ=255)
                                                     aaTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACg                                                                                                                                                                                    >  2:168198/1‑251 (MQ=255)
                                                          gACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGCTATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGAtgttg                                                                                                                                                                              >  2:92460/1‑251 (MQ=255)
                                                            cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGAtgttgt                                                                                                                                                                             >  1:34176/1‑251 (MQ=255)
                                                              cAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTAc                                                                                                                                                                           >  2:313548/1‑251 (MQ=255)
                                                                gagCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGCTGATGTGAATGATCTAGCCCGGCGAAGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTg                                                                                                                                                                         >  1:184551/1‑251 (MQ=255)
                                                                  gCACAAATTATATTTCGAATGAAAGCAAGGATGAAATTGCTGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGTGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGCCACACACTACCAGATCGCGATCCTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGa                                                                                                                                                                        >  2:291829/1‑250 (MQ=255)
                                                                     aaaaTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAatat                                                                                                                                                                    <  1:104378/250‑1 (MQ=255)
                                                                      aaaTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATg                                                                                                                                                                   <  1:47409/251‑1 (MQ=255)
                                                                       aaTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGc                                                                                                                                                                  <  2:208396/251‑1 (MQ=255)
                                                                                          gTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCaa                                                                                                                                               <  2:281460/251‑1 (MQ=255)
                                                                                                 tGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGANGGt                                                                                                                                        >  2:405248/1‑251 (MQ=255)
                                                                                                  gAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGt                                                                                                                                        >  1:416623/1‑250 (MQ=255)
                                                                                                         atgatgGGAATGATTAAGCCCGGCGCCGACGCCGCCGGGCCGAGGAGATTACAGCGTATGTTCAGTACGGGCGATGATCTCGTCCTGAGCGTCTGGAGCGCGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTCCCAGATCGCGATACTGCTCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGTGACACGATGCTGTACTTAATATGCGAGCCTTTCTGCACTTCAACGAAGGTACgcagc                                                                                                                                 >  2:345020/1‑250 (MQ=255)
                                                                                                             tgtgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGat                                                                                                                            >  1:311144/1‑251 (MQ=255)
                                                                                                                            ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTg                                                                                                                                                                                                                               >  2:26264/1‑137 (MQ=255)
                                                                                                                            ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTg                                                                                                                                                                                                                               <  1:26264/137‑1 (MQ=255)
                                                                                                                               gCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTc                                                                                                          <  1:60357/251‑1 (MQ=255)
                                                                                                                                       gccgccGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATTCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGt                                                                                                  <  1:413830/251‑1 (MQ=255)
                                                                                                                                                ccGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCg                                                                                         <  2:311144/251‑1 (MQ=255)
                                                                                                                                                  gaggagATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGtt                                                                                       >  1:266381/1‑251 (MQ=255)
                                                                                                                                                  gaggagATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGtt                                                                                       >  1:131085/1‑251 (MQ=255)
                                                                                                                                                   aggagATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTg                                                                                      >  2:72291/1‑251 (MQ=255)
                                                                                                                                                       gaTTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCgg                                                                                  >  2:115444/1‑251 (MQ=255)
                                                                                                                                                            cAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCAGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTc                                                                              >  2:454291/1‑250 (MQ=255)
                                                                                                                                                                   tGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTg                                                                      >  2:368240/1‑251 (MQ=255)
                                                                                                                                                                    gtaaagcaCGGGCGATGATATGCTCCTGAGCGACTAGAGAGAGCACGGTGACTAACGGGGAATGGCTCGCGACACACACTACCAGATCGCGAACCTGATCGGGATGTTTTTTCGCATCAAGCAGCGTTCCGCGGGAAACGATGTTGTACTGAATAAGCAAGTCTTTATGCACTGCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCCTTGCAGGAATCacgttgtga                                                                     <  2:305846/244‑9 (MQ=255)
                                                                                                                                                                         ggagtggCGATGATATCGTCCTGATGTGCTGGAGAGAGCGCGGTGAAGAACGCGGAATATCCCGCGACACACACTACCAGATCGCGATACTGTTCCGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATCTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCATCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCGCCAGCGTTGCCGGATTCAGTTTCTGGtttag                                                                <  2:20214/246‑4 (MQ=255)
                                                                                                                                                                           aCGGGCGAGGATATCGTCCTGAGCGTCTGGAGAGAGCGAGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGAGATACTGATAGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGAACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGATTGTCAGATTCGTTATCCAGCGTTGCCGGATTAAGTTTCTGGTTGAga                                                               >  1:237877/1‑249 (MQ=255)
                                                                                                                                                                                 gatgatATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACAcgc                                                        >  2:442642/1‑251 (MQ=255)
                                                                                                                                                                                  atgatATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCCGCAGGATCATCAGTTTCTGCTTGTCAGCTTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACAcgcc                                                       >  2:129056/1‑251 (MQ=255)
                                                                                                                                                                                    gatATCGTCCTGAGCGTCTGGAGATAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGGACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCTGATTCGTTCTCCAGCGTTGCCGTATTCAGTTTCTGGTTGAGCAACAcgccgc                                                     <  1:305797/251‑1 (MQ=255)
                                                                                                                                                                                    cttatCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTAt                                                                                                                                                        <  2:349316/150‑1 (MQ=255)
                                                                                                                                                                                    cttatCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTAt                                                                                                                                                        >  1:349316/3‑152 (MQ=255)
                                                                                                                                                                                               gAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGAt                                                                                                                                                                                  >  2:241370/1‑115 (MQ=255)
                                                                                                                                                                                               gAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGAt                                                                                                                                                                                  <  1:241370/115‑1 (MQ=255)
                                                                                                                                                                                               gAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCt                                          <  2:232682/251‑1 (MQ=255)
                                                                                                                                                                                                         gagagaGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGc                                 >  1:235906/1‑251 (MQ=255)
                                                                                                                                                                                                           gagagCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGc                                 >  2:448619/1‑249 (MQ=255)
                                                                                                                                                                                                            agagCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGTGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCTGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGGCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACTCGCCGCCGAGAATCGCTGCCGTAGGCAGt                              <  1:424536/251‑1 (MQ=255)
                                                                                                                                                                                                            agagCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGt                              >  1:334800/1‑251 (MQ=255)
                                                                                                                                                                                                             gagCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGtt                             <  2:131085/251‑1 (MQ=255)
                                                                                                                                                                                                                cgcgGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTAc                          <  1:448619/251‑1 (MQ=255)
                                                                                                                                                                                                                    gTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTACCCAc                      <  1:390435/251‑1 (MQ=255)
                                                                                                                                                                                                                          aaagcgGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGACGTAGGCAGTTTACCCACTGAGcc                >  2:391702/4‑251 (MQ=255)
                                                                                                                                                                                                                            cgcgGAATAGCCCGCGACACACACTACCAGATCGCGTTACTGATCGGGAAGTTTTTTCCCATCCAGCAGCGTTTCGCGGGAAACGCTGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGCGCAACACGCTGCCGAGAATCGCTGCCGTAGGCAGTTTACCCACTGAGCCa               >  2:314334/1‑250 (MQ=255)
                                                                                                                                                                                                                               ggAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTACCCACTGAGCCAATg            >  1:303499/1‑250 (MQ=255)
                                                                                                                                                                                                                                  aTAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGCTCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTACCCACTGAGCCAATGACCg        >  1:339918/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTACCCACTGAGCCAATGACCGCAg     <  1:256306/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cccGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTACCCACTGAGCCAATGACCGCAGt    >  1:214964/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cgcgACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTACCCACTGAGCCAATGACCGCAGTAg  <  1:313548/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
ACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTACCCACTGAGCCAATGACCGCAGTAG  >  NC_000913/860021‑860502

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