Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,280,131 A→G 12.7% E51G (GAG→GGG)  yjcE → putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,280,1310AG12.7% 172.1 / 2.7 71E51G (GAG→GGG) yjcEputative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (39/20);  new base G (8/2);  total (49/22)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCT  >  NC_000913/4279899‑4280363
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                         
ttttGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGG‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTg                                                                                                                                                                                                                                      <  1:645467/238‑1 (MQ=255)
 tttGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGGGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGTTGCGCTACTGGCGTGGCCGCGGTTTGGTTTGCATGTG‑GAGTTTGa                                                                                                                                                                                                                                     >  1:188621/1‑238 (MQ=255)
  tttggATATTTTGAAAAAAATTAGAATTCAGGCATAACGTTTTTCTCAGATCAAAAAAAAAATCGCAAACAGGGAACGCATGGAAATCTTCTTCCCCATACTGATAATGCCCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATc                                                                                                                                                                                                                                   <  2:454589/235‑1 (MQ=255)
  ttGCGATATTTGGAAAAAAATGAGAATTCCGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATc                                                                                                                                                                                                                                   <  1:264406/239‑1 (MQ=255)
  ttGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GGGTTTGAt                                                                                                                                                                                                                                    >  1:227872/1‑238 (MQ=255)
   tGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGTTGCGCTACTGGCGTGGCCGACGTTTGGTTTGAATGGG‑GAGTTTGAt                                                                                                                                                                                                                                    >  1:520146/1‑237 (MQ=255)
   tGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GGGTTTGAt                                                                                                                                                                                                                                    >  1:248045/1‑237 (MQ=255)
   tGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGAt                                                                                                                                                                                                                                    >  1:398636/1‑237 (MQ=255)
        tatTCTGAAAAAAATTGGAATTCAGGCATAACGGCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCGTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑Ac                                                                                                                                                                                                                            <  1:570512/239‑1 (MQ=255)
        tatTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGa‑a                                                                                                                                                                                                                             <  1:196928/238‑1 (MQ=255)
             tGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCttt                                                                                                                                                                                                                       <  1:673261/239‑1 (MQ=255)
               aaaaaaaTGAGACTTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGCTAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCACGCTTCCGCTTATGCAAATCGCCATCGGTGCGTTACTGGCGTGGCCGACGTATGGTTT‑CATGTGTGAGGTTGATCCTGAGACTCtttt                                                                                                                                                                                                                      >  2:65607/1‑239 (MQ=255)
                    aaTGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GGGTTTGATCCTGA‑ACTCCTTTTAGTCt                                                                                                                                                                                                                >  1:661320/1‑239 (MQ=255)
                                    tAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGtt                                                                                                                                                                                                >  1:448976/1‑239 (MQ=255)
                                    tAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGCCCGACGTTTGGTTTGCCTGTG‑GAGTTTGGTCCTGA‑ACCCTTTTTAGTCTTGTTTATCccgcctt                                                                                                                                                                                                 >  1:388991/1‑236 (MQ=255)
                                       cGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCt                                                                                                                                                                                             >  1:267374/1‑239 (MQ=255)
                                             tctcAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATc                                                                                                                                                                                                                                   <  1:370945/196‑1 (MQ=255)
                                                   aTCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGCCGTTTGGTTTGCATGTG‑GGGTTTGATCCTGA‑ACTCTTTTTAGGCTTGTTTATCCCGCCGTTGCTGTTCGCttgtgg                                                                                                                                                                                 >  1:826828/1‑234 (MQ=255)
                                                               ctcccaacgcgggACCGCAAGGAAACTTTCTTAACCATTTTGATAATGGCCCTGGGGGGTCCGGTGT‑CTGGGGGGGTCACG‑CGTGGCATCCCCGGTCAGACCCCGCTGCCGTTTATGCAACTCGCAATCGGTGCGCTACGGGCGTGGCGGACGGTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGcc                                                                                                                                                                     <  2:730064/229‑1 (MQ=255)
                                                                 cGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑Ac                                                                                                                                                                                                                            >  2:109637/1‑182 (MQ=255)
                                                                 cGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑Ac                                                                                                                                                                                                                            <  1:109637/182‑1 (MQ=255)
                                                                 cGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGATGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCCTGTG‑GAGTTTTATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCCCTGATGGCTGGGAAACCCCg                                                                                                                                                                    >  2:328102/1‑238 (MQ=255)
                                                                           gAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAAtt                                                                                                                                                         >  1:203332/1‑239 (MQ=255)
                                                                             ccGCATGGAACTCTGCTTCACCATACTGATAGTGCCCCTCGTGGTCTCGCTGC‑CCGGGGTGGTCACC‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAAt                                                                                                                                                          <  1:43531/235‑1 (MQ=255)
                                                                                cATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAAAGCCGACCCGTGAATTTCTCg                                                                                                                                                    >  1:713438/1‑239 (MQ=255)
                                                                                  tGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GGGTTTGATCCTGG‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGGTGGCTGGAAAACGCCGAACCGGGAAttttcggaa                                                                                                                                                  >  1:231852/1‑233 (MQ=255)
                                                                                       aTCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTTTCCCGCCGTTGCTGTTCGCTGTTGGGTGGAAAACGCCGCCCCCGGAATTtctccacaa                                                                                                                                                >  2:43531/1‑231 (MQ=255)
                                                                                               cACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTATCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCgaga                                                                                                                                       >  1:745873/1‑237 (MQ=255)
                                                                                                  cATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGCGCGAttt                                                                                                                                  >  1:333723/1‑239 (MQ=255)
                                                                                                         aTAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCTGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATTTG‑GCGTTTGTTCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAAGCCCGACCCGTGAATTTTCCGAACATGGGCGGGGGTttttgtg                                                                                                                              >  2:837886/1‑233 (MQ=255)
                                                                                                         aTAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGcc                                                                                                                            >  1:649776/1‑238 (MQ=255)
                                                                                                         aTAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGcc                                                                                                                            >  1:652648/1‑238 (MQ=255)
                                                                                                         aTAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGCCTTGTTTTTCCCGCCGTTGCTGTTCGCCGATGGCGTGAAAACGCCGACCCGGGAAATTTTCAAACATGGTCGAGAGATTTTCGGccc                                                                                                                           >  2:107869/1‑238 (MQ=255)
                                                                                                           aaTGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGGTTGATCCTGA‑ACTCTTTTTAGTCTTTTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGAACCGTGAATTTTTCGAACATGGTAGGTGGATTTTCGGCCTc                                                                                                                          >  1:523374/1‑238 (MQ=255)
                                                                                                              gACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGGGGGATTTTTGCCCTCtcgc                                                                                                                      >  1:568097/1‑239 (MQ=255)
                                                                                                              gACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGGTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTAACTTTTTGCAACATGGTCGCGGGATTTTTGGCCTccc                                                                                                                        >  1:568111/1‑235 (MQ=255)
                                                                                                                   tCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTCCTGGCGTGGCCGCCGTTTGGTTTGCATGTG‑GGGTTTGGTCCTGA‑ACGCTGTTTAGTACTTTTTTTCCCGCCGTTGCTGTTCGCTGTTAGCCGGAAAACACCGGCGCGTAAAGTTTCCGGAAATGGGGGGGAGAATTTCGGCCTcgcggggggg                                                                                                                 >  2:751653/1‑233 (MQ=255)
                                                                                                                   tCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GCGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTTTCGAAAATGGGCGAGAGCTTTTCGGCCACtcgctgtc                                                                                                                  >  1:344315/1‑236 (MQ=255)
                                                                                                                   tCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGGAAACGCCGACCCGTGAATTTCTCGAACATGGTCGGAGGATTTTCGGCCTccgcctggc                                                                                                                  >  1:892510/1‑230 (MQ=255)
                                                                                                                   tCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgctggc                                                                                                                  >  1:250421/1‑238 (MQ=255)
                                                                                                                    cGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGACCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GGGTTTGATCCTGA‑ACTCTTTTTAGTCTTTTTTATCCCGCCCTTGCTTTTCGCGGCTGGCCGTAAAACGCCGACCCGTTAATTTCTCGAACATGGTCGCGGGGTTTTTGGCCCCcgcctggcgc                                                                                                                >  2:2982/1‑239 (MQ=255)
                                                                                                                     cgtttctcGCTTTCCCGGGGTGGTCCCG‑CGTTTCCGGCCTTTTCAGATCCCGCTTCCGTTTATGCAAATCCCCATCGGTGCTCTCGTGGCGTGGCCGACTTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgctggcgc                                                                                                                <  2:713438/235‑1 (MQ=255)
                                                                                                                           tcgcctgccgggggggGTCTCT‑GGTTTCAATCCTTTCCAGATCCCTCTTCCGCTTAGGCAAATCGCCATGGGTGCGCTATGGGCGTGCCCGACGTTGGTTTGGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTCAGTCTGGTTCATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtg                                                                                                           <  2:826828/228‑1 (MQ=255)
                                                                                                                           tcGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGGGATTTTCGGCCTCGCCCTGGGGTtgggggg                                                                                                         >  1:582988/1‑238 (MQ=255)
                                                                                                                               gtt‑CCGGGGTTGTTTCTCCGTGTCATGCCTTTTGAGTCCCCGTTTCCGTTTAGGCAAATCGCCATCGGTGCGCTACTGGCGGGGCCGACTTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTTTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTggtggtggt                                                                                                      <  2:248045/237‑1 (MQ=255)
                                                                                                                                    cGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGGGCTGTGGGGGGGCACCg                                                                                                 >  1:454354/1‑239 (MQ=255)
                                                                                                                                       ggtggtCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTTTCGAACATGGGCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGGGGTCACCGTgg                                                                                              >  1:581353/1‑239 (MQ=255)
                                                                                                                                        gtggtCACC‑CGTGTCATGCCCTTTCAAATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTTCTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTgg                                                                                              <  1:140766/238‑1 (MQ=255)
                                                                                                                                            tCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTTTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGGGGGGGTTACCGTGGTCg                                                                                           >  1:863371/1‑237 (MQ=255)
                                                                                                                                             cACT‑CGTGTTCTGCCCTTTCAGTTCCCGTTTCCTGTTTTGCAAATCGCCAGCGGTGCGCTACGGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCa                                                                                        <  1:525770/239‑1 (MQ=255)
                                                                                                                                             cACC‑CGTGTCATGCCCTTACAGACCCCCGTTCCGCTTAGGCAAATCGCCATCGGTGCGCTACTGGCGTGTCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCa                                                                                        <  1:220103/239‑1 (MQ=255)
                                                                                                                                              aCT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGGGGTAACCGTGGTCGGCAt                                                                                       >  1:797691/1‑239 (MQ=255)
                                                                                                                                                         tccccTTTCAGGCCCTGTTTCCGGTTATGCAAATCGCCATGGGTGCGCTACGGGGGGGGCCGATTTTTGGTTTGCTTGTG‑GGTTTTGATCCTGA‑ACTTTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCtt                                                                             <  2:203332/237‑1 (MQ=255)
                                                                                                                                                           cctttgctgttcccgttcccgtTTTTGCAATTGGCCTGCGGGGCGCTGTTGGCGTGCCGTACTTTTGTTTGTCTGGGG‑GAGTTTGTTCCTGA‑ACCTTTTTTAGTCTTGTTTCTCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTAt                                                                           <  2:388991/229‑1 (MQ=255)
                                                                                                                                                              tttCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGGGGGGGTCACCGTGGTCGGCATTGGCTTCCTTAtttt                                                                        >  1:243830/1‑238 (MQ=255)
                                                                                                                                                              ttccgaaccccGCTTCCGTTTATGCAAATCGCAATCGGTGCGCTACTGGCGTGGCGGCCGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTa                                                                        <  1:76797/231‑1 (MQ=255)
                                                                                                                                                                    aTCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGCCCCGTGAATTTCTCGAACATGGTCGCGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGCTGGCATTGGCTTCCTTTTTTTCTgggg                                                                  >  1:533/1‑238 (MQ=255)
                                                                                                                                                                        cgcttccgcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGTAAACGCCGACCCGTTAATTTCTCGAACATGGTCGAGAGATTTTCCGCCTCGCGCTGGCGCTGGTGGTGGTCCCCGTGGTCGGCATTGGCTTCCTTCTTTACTTGggggtg                                                              >  2:872127/1‑239 (MQ=255)
                                                                                                                                                                          cttccgcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGGGATTTTCGGCCTCGGCCTGGCGCTGGTGGTGGTCACCTGGGTCGGCATTGGCTTCCTTATTTAGTggggggggag                                                            >  2:419831/1‑235 (MQ=255)
                                                                                                                                                                          cttccgcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGGGGTGcc                                                            >  1:113685/1‑239 (MQ=255)
                                                                                                                                                                          cttccgcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAAAGCCGACCCGTGAATTTCTCGAACATGGTCGAGGGATTTTCGTCCTCGCGCTGGGCGTGGTGGTGGTGACCTCGGTCGGCAGTGGCTTCCTTATTTACTGggtgtgtcc                                                            >  1:123609/1‑234 (MQ=255)
                                                                                                                                                                               cttatcgCAAATCCCCATCGTTACCTTCCTGGCGGCTCCGGCGTTTGTTTTTCTTGTG‑TGTTTTGATCCTGC‑ACTCTTTGTAGTCTTTTTTCTCCCGCCGTTGCTGTTCGGTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCa                                                       <  2:188621/233‑1 (MQ=255)
                                                                                                                                                                                                              gCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGGGCCGGGCATTCCCCTGATCCCTGCCTTTTCGCTggcggc                        >  1:165992/1‑239 (MQ=255)
                                                                                                                                                                                                              gCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGGGGTCACCGTGGTCGGCATTGGCTTCCTTCTTTACTGGGTGGTGCCGGGCCTTCCGCTGCTCCCTGCCTTTTCCCTgggggc                        >  1:565541/1‑239 (MQ=255)
                                                                                                                                                                                                                   gccggctgttggttggGATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTg                    <  1:675928/230‑1 (MQ=255)
                                                                                                                                                                                                                   gCCGGCGTTTGGGTTGCAAGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACTCCGACCCGTGAATTTCTCGCACAGGGTCGCGCGATTTTCGGCCTCCCGTCTGCGCTAGTGGTGGTCACCGGTGTCGGCATTGGCTTTCTGATTTACTGGGTGGGGCCGGGCATTCCCCTGGTGCCCTCATTTGCGCCGGgggaggggg                   >  2:478991/1‑236 (MQ=255)
                                                                                                                                                                                                                   gCCGACGTTTGGTTTGCATGTG‑GGGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGt                     >  1:582466/1‑237 (MQ=255)
                                                                                                                                                                                                                   gCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGGGCCGGGCATTCCGCTGATCCCTGCCTTTCCGCTTGCGGGGGTgg                   >  1:158180/1‑238 (MQ=255)
                                                                                                                                                                                                                     cGACTTTTGGTTTGAATGTG‑GAGTTGGATCCTTA‑ACTCTTTTGAGTTTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCt                  <  1:751653/238‑1 (MQ=255)
                                                                                                                                                                                                                                    aTGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGGGATTTTTGGCCTCGCGCTGGCGCTGGTGGGGGTCACCGTGGTCGGCATTGGCTTCCTTGTTTCCTGGGTGGTGCCGGGCTTTCCGCGGATCCCTGCCTTTGCGCTGGCGGCGGGGCCTTCTcctcgcgctatt  >  2:31851/1‑229 (MQ=255)
                                                                                                                                                                                                                                          gAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTggtggtggt                                                                                                      >  1:62102/1‑134 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                         
TTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACCCTCGTGGTCTCGCTGT‑CCGGGGTGGTCACT‑CGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTG‑GAGTTTGATCCTGA‑ACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCT  >  NC_000913/4279899‑4280363

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