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

ACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTACCCACTGAGCCAATGACCGCAGTAG  >  NC_000913/860021‑860502
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
aCGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCAGCGACACACACTACCAGa                                                                                                                                                                                                                                          >  1:23567/1‑250 (MQ=255)
      aTAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGAt                                                                                                                                                                                                                                   <  2:229389/251‑1 (MQ=255)
                 ccTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATCTTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCggg                                                                                                                                                                                                                         >  1:213757/1‑250 (MQ=255)
                 ccTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCggg                                                                                                                                                                                                                         >  1:232682/1‑250 (MQ=255)
                         gAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAAGTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGtttttt                                                                                                                                                                                                                <  1:110151/251‑1 (MQ=255)
                           aaaTGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCg                                                                                                                                                                                                              <  2:2415/251‑1 (MQ=255)
                           aaaTGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCg                                                                                                                                                                                                              <  1:430073/251‑1 (MQ=255)
                            aaTGAAACCGTACCAGTTTACTAACCATCTGACGCAGAGAACAACTTATCTTTCGAATGCAAGTAAGGATGACATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGCTATTCCAGCTTATGTTCAGTACAGGCGATGATATCTTCCTGAGCGT‑CTGTAGCGAGCGCTGGGACGAACCCGGACTAGCCCGCGCCCCCGACTACCAGATCGCGATACTGATCGGGCTGttttttagc                                                                                                                                                                                                             >  2:444275/1‑248 (MQ=255)
                               gAAACCGTAACAGTGCAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAAGAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCAt                                                                                                                                                                                                           <  1:426518/250‑1 (MQ=255)
                                aaaCCGTAACAGTGTAATATCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATcc                                                                                                                                                                                                         <  2:213757/251‑1 (MQ=255)
                                    cGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGCAAGGATGAAATTGATGATGTGAATGATTTAGCCAGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCcagc                                                                                                                                                                                                      >  2:147873/1‑250 (MQ=255)
                                            tgtAATAACAATGTGACGCAGCGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGCTGATATCGTCCTGAGCGT‑CTGGCGAGAGCGCGGTGCAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTcg                                                                                                                                                                                             >  2:208001/1‑251 (MQ=255)
                                             gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTcg                                                                                                                                                                                             <  1:5767/250‑1 (MQ=255)
                                                     aaTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGaaa                                                                                                                                                                                      >  1:284007/1‑249 (MQ=255)
                                                     aaTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACg                                                                                                                                                                                    >  2:168198/1‑251 (MQ=255)
                                                          gACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGCTATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGAT‑GGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGAtgttg                                                                                                                                                                              >  2:92460/1‑251 (MQ=255)
                                                            cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGAtgttgt                                                                                                                                                                             >  1:34176/1‑251 (MQ=255)
                                                              cAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTAc                                                                                                                                                                           >  2:313548/1‑251 (MQ=255)
                                                                gagCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGCTGATGTGAATGATCTAGCCCGGCGAAGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTg                                                                                                                                                                         >  1:184551/1‑251 (MQ=255)
                                                                  gCACAAATTATATTTCGAATGAAAGCAAGGATGAAATTGCTGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGTGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGCCACACACTACCAGATCGCGATCCTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGa                                                                                                                                                                        >  2:291829/1‑250 (MQ=255)
                                                                     aaaaTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAatat                                                                                                                                                                    <  1:104378/250‑1 (MQ=255)
                                                                      aaaTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATg                                                                                                                                                                   <  1:47409/251‑1 (MQ=255)
                                                                       aaTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGc                                                                                                                                                                  <  2:208396/251‑1 (MQ=255)
                                                                                          gTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCaa                                                                                                                                               <  2:281460/251‑1 (MQ=255)
                                                                                                 tGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGANGGt                                                                                                                                        >  2:405248/1‑251 (MQ=255)
                                                                                                  gAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGt                                                                                                                                        >  1:416623/1‑250 (MQ=255)
                                                                                                         atgatgGGAATGATTAAGCCCGGCGCCGACGCCGCCGGGCCGAGGAGATTACAGCGTATGTTCAGTACGGGCGATGATCTCGTCCTGAGCGT‑CTGGAGCGCGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTCCCAGATCGCGATACTGCTCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGTGACACGATGCTGTACTTAATATGCGAGCCTTTCTGCACTTCAACGAAGGTACgcagc                                                                                                                                 >  2:345020/1‑250 (MQ=255)
                                                                                                             tgtgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGat                                                                                                                            >  1:311144/1‑251 (MQ=255)
                                                                                                                            ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTg                                                                                                                                                                                                                               >  2:26264/1‑137 (MQ=255)
                                                                                                                            ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTg                                                                                                                                                                                                                               <  1:26264/137‑1 (MQ=255)
                                                                                                                               gCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTc                                                                                                          <  1:60357/251‑1 (MQ=255)
                                                                                                                                       gccgccGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATTCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGt                                                                                                  <  1:413830/251‑1 (MQ=255)
                                                                                                                                                ccGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCg                                                                                         <  2:311144/251‑1 (MQ=255)
                                                                                                                                                  gaggagATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGtt                                                                                       >  1:266381/1‑251 (MQ=255)
                                                                                                                                                  gaggagATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGtt                                                                                       >  1:131085/1‑251 (MQ=255)
                                                                                                                                                   aggagATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTg                                                                                      >  2:72291/1‑251 (MQ=255)
                                                                                                                                                       gaTTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCgg                                                                                  >  2:115444/1‑251 (MQ=255)
                                                                                                                                                            cAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCAGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTc                                                                              >  2:454291/1‑250 (MQ=255)
                                                                                                                                                                   tGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTg                                                                      >  2:368240/1‑251 (MQ=255)
                                                                                                                                                                    gtaaagcaCGGGCGATGATATGCTCCTGAGCGA‑CTAGAGAGAGCACGGTGACTAACGGGGAATGGCTCGCGACACACACTACCAGATCGCGAACCTGATCGGGATGTTTTTTCGCATCAAGCAGCGTTCCGCGGGAAACGATGTTGTACTGAATAAGCAAGTCTTTATGCACTGCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCCTTGCAGGAATCacgttgtga                                                                     <  2:305846/244‑9 (MQ=255)
                                                                                                                                                                         ggagtggCGATGATATCGTCCTGA‑TGTGCTGGAGAGAGCGCGGTGAAGAACGCGGAATATCCCGCGACACACACTACCAGATCGCGATACTGTTCCGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATCTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCATCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCGCCAGCGTTGCCGGATTCAGTTTCTGGtttag                                                                <  2:20214/246‑4 (MQ=255)
                                                                                                                                                                           aCGGGCGAGGATATCGTCCTGAGCGT‑CTGGAGAGAGCGAGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGAGATACTGATAGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGAACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGATTGTCAGATTCGTTATCCAGCGTTGCCGGATTAAGTTTCTGGTTGAga                                                               >  1:237877/1‑249 (MQ=255)
                                                                                                                                                                                 gatgatATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACAcgc                                                        >  2:442642/1‑251 (MQ=255)
                                                                                                                                                                                  atgatATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCCGCAGGATCATCAGTTTCTGCTTGTCAGCTTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACAcgcc                                                       >  2:129056/1‑251 (MQ=255)
                                                                                                                                                                                    gatATCGTCCTGAGCGT‑CTGGAGATAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGGACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCTGATTCGTTCTCCAGCGTTGCCGTATTCAGTTTCTGGTTGAGCAACAcgccgc                                                     <  1:305797/251‑1 (MQ=255)
                                                                                                                                                                                    cttatCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTAt                                                                                                                                                        <  2:349316/150‑1 (MQ=255)
                                                                                                                                                                                    cttatCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTAt                                                                                                                                                        >  1:349316/3‑152 (MQ=255)
                                                                                                                                                                                               gAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGAt                                                                                                                                                                                  >  2:241370/1‑115 (MQ=255)
                                                                                                                                                                                               gAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGAt                                                                                                                                                                                  <  1:241370/115‑1 (MQ=255)
                                                                                                                                                                                               gAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACACACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCt                                          <  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)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
ACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGT‑CTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATTCAGTTTCTGGTTGAGCAACACGCCGCCGAGAATCGCTGCCGTAGGCAGTTTACCCACTGAGCCAATGACCGCAGTAG  >  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: ‑