Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,433,459 C→T G189D (GGC→GAC)  qorB ← NAD(P)H:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,433,4590CT98.6% 215.3 / NA 71G189D (GGC→GAC) qorBNAD(P)H:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (54/16);  total (54/17)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGCCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATG  >  NC_000913/4433215‑4433709
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
tGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGt                                                                                                                                                                                                                                                       >  1:128111/1‑251 (MQ=255)
ggggTGCCCAATCAATTTGCTAAGTTATTTTCGTTCACCAAACAGCCCCCCTTTCGAGGCGCCACCGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGt                                                                                                                                                                                                                                                       <  1:458805/250‑1 (MQ=255)
 gggaggCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTg                                                                                                                                                                                                                                                      >  1:304746/5‑251 (MQ=255)
      cccAATCATTTTGTTAGGTGTTTGCCTGTCACAAACCGACCCGCCTTTCGATGCCCCAACGCAAGAACCCGCCACAATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTGCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCgg                                                                                                                                                                                                                                                 <  1:267754/250‑1 (MQ=255)
         aatcaatTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCtt                                                                                                                                                                                                                                              >  1:325481/1‑251 (MQ=255)
         aatcaatTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGCCAGCTCTGCCGCTAACTGTGTCAACGTCCAGTCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGcct                                                                                                                                                                                                                                              >  2:282493/1‑249 (MQ=255)
         aatcaatTTGAAAAGTGTTTTGCTGTCATCAACCAGACCGCCTTTCGATGCGCCAACGTCAGAATCAGCCAGCATATCCGCCAGTCCGTCGGGCCGTCCGACGCTTTTCAGCGCCGCCGAGCAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCCCTGTTTGGTCACCTCTGCGGCAAACTGGTCCAACGTCCCGGCACTATCGCCCGCCAGTTCGTAAACATTGTCTTCGAGAcagacgt                                                                                                                                                                                                                                              >  2:470312/1‑245 (MQ=255)
                gggtgatgtgtgtggcttccaccaAACAGCCCGCTTTTGGAGGCGCAAACGTCAGAATCCCCCAGAATACCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGTTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCTCTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAAt                                                                                                                                                                                                                                       <  1:425431/229‑1 (MQ=255)
                   cTAAGTGTTATGCTGTCATAAAACATACCGCCTTTCGATGCGCCAACGTAAGAATCCGCCAGCATATCAGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGAGGCGAAATCGGCTTCGCTCTGATTTTGATAGGAAACCTGTTTGCAGCTCTGTTTGGTCAGCTCTGCAGCTAACTGTGTCAACGTCCAGGACCTCACGCCCGCCCTTTCGTAAACCTTGTCTTCGTGACCGGCtaccgtaatcac                                                                                                                                                                                                                                    >  2:370956/1‑240 (MQ=255)
                    tAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCCGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCAcc                                                                                                                                                                                                                                   >  2:484421/1‑250 (MQ=255)
                        tgtTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGACAGCATATCCGCCAGTCCGTCGGGAAGTCCGACGCTTTTCAGTGCAGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTg                                                                                                                                                                                                                               >  1:359324/1‑251 (MQ=255)
                             tGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGc                                                                                                                                                                                                                          >  2:276453/1‑250 (MQ=255)
                                tGTCATCAAACAGACCGCCTTTCGAGGCGCCAACGTCAGAATCCGCAAGCATATCCGCAAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCCGATTTTGATAGGTAACCTGTTTGCCGCTCTGATTGGTCCGCTCTGCCGCTAACTGTGTCAAATTCCAGGCCCTATCGCAAGCCAATTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCAc                                                                                                                                                                                                                                    >  2:531216/1‑238 (MQ=255)
                                gtccaacaaactgcccGCCTTTCGAGGCCCCAACGTCAGAACCCCCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAAGCGGCTTCGCTCAGATTTTGATAGGTAACCGGTTTGCCGCTCTTTTTTGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCAc                                                                                                                                                                                                                                    <  1:531216/224‑1 (MQ=255)
                                  tcatcaCACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCTCCATCATATCCTCCACTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCCGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTCATCACGCGTGCCGCAGCTGc                                                                                                                                                                                                                    >  2:184707/1‑251 (MQ=255)
                                       aaaCAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCAt                                                                                                                                                                                                               >  2:285779/1‑251 (MQ=255)
                                          cAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCTGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCCCGCGTCCCGCAGCTGCCGCATAAt                                                                                                                                                                                                            >  2:154967/1‑251 (MQ=255)
                                               cGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGGTAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAAACTTGTCTTCGTGACCGGCTTCGCTACTACCGCCTGCCGCAGCTGCCGCATAATCTGcc                                                                                                                                                                                                       >  2:87350/1‑251 (MQ=255)
                                                             gcCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGACCGCTTTGCTgggg                                                                                                                                                                                          >  2:333075/1‑248 (MQ=255)
                                                                    tCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGAtttt                                                                                                                                                                                   >  1:219043/1‑251 (MQ=255)
                                                                              ccAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATcgccc                                                                                                                                                                         >  1:321866/1‑251 (MQ=255)
                                                                               cAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATcgcccc                                                                                                                                                                        >  1:440623/1‑250 (MQ=255)
                                                                                 acATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGACATCGGCTTCGCTCCGCTTTTGATAGGTACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCACCGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGAGACCGGCTTCGCTAATCCCGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGCTTCTGAGGCGATTTTGCCATcccccgccc                                                                                                                                                                     >  2:26654/2‑250 (MQ=255)
                                                                                   atatCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCCGTTCGTAAACCTTGTCTTAGTGACCGGCTTCGCTACTCACCCGTGCCGCAGCTGCCGCCTAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCccgccgg                                                                                                                                                                    >  2:220417/1‑250 (MQ=255)
                                                                                     atCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCAAGCGGGCCGCAGCTGCCGCCTAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCggccgccc                                                                                                                                                                  >  2:428374/1‑249 (MQ=255)
                                                                                        cGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCCTACTCTGCCCGCGTTGCTGAGGCGATTTTGCCATCCCCCGCCGCAcca                                                                                                                                                                >  2:526396/1‑249 (MQ=255)
                                                                                         cccggtcctgcgggccgtccgccctttttCAGTGCCGCGGCAAATTCGCCTTCCCTCAGTTTTTGATAGGTACCTTGTTTGCCCCTCTGTTGGGTCAGCTCTGCCGCTAACTTTCTCAACGCCCAGGCACTATCCCCCGCCAGTTCGTAACCCTTGTCTTCGTGACCGGCTTCGCTAATCACGCTTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGAt                                                                                                                                                              <  1:106126/226‑1 (MQ=255)
                                                                                           cAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCAAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTGGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAACCCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATaa                                                                                                                                                            <  1:523794/251‑1 (MQ=255)
                                                                                               ccGTCGGGCAGTCCGCCGCTTTTCAGTGCCGCGACGAAATCGGCTTCTCTCAGCTTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCACCGTCCAGGCCCTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCAGCAGCTGCCGCATAATCTGCCCGCGTTGCGGAGGCGATTTTGCCATCGCCCGCCCCCCCGATAAAtcc                                                                                                                                                        >  2:296876/1‑249 (MQ=255)
                                                                                                    gggCAGTCCGACGTTTTTCAGTGCCGCGGCAAAACCGGCTTCGCTCAGATTTTGATAGGTAACCGTTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCgt                                                                                                                                                   <  1:484421/251‑1 (MQ=255)
                                                                                                    gggCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCgt                                                                                                                                                   >  2:447405/1‑251 (MQ=255)
                                                                                                       cAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGtt                                                                                                                                                >  1:201364/1‑251 (MQ=255)
                                                                                                        ggTCCGAGGCTTTTCGGGGCCGCGGCGAAATCGGCTCGGCTCAGATTTGGTTAGGTAACCTGTTGGCCGCTCTTTTTGGCCACCCCGGCCGCAAACTGTGGCAACGTCCAGGCACTACCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTc                                                                                                                                               <  1:448565/250‑1 (MQ=255)
                                                                                                              aCGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCCTAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATCCGCCGCGTTCCAgggc                                                                                                                                         >  2:287141/1‑248 (MQ=255)
                                                                                                               cGCTTTCCAGTGCCGCGGCGAAATGGGCTTCGCTAAGTTTTTGATAGGAACCCTGTTTCCCGCTCTGTTTGGTCAGCTCTGCCGCTACCTGTTTCAACGCCAAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGtgct                                                                                                                                        <  1:361212/251‑1 (MQ=255)
                                                                                                                  ttttCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCAGCACTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGcc                                                                                                                                     >  1:352049/1‑251 (MQ=255)
                                                                                                                      cAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGAGTTCCAGTGCTGCCgggg                                                                                                                                 >  2:476706/1‑251 (MQ=255)
                                                                                                                             gcggcgAAATCGGCTTCGCTCATATTTTGATAGTAAACCTGTTTGCCGCTCTGTTGGGTCAGCCCTGCCGCTAACTGTGTCAACGCCCAGGCCCTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGc                                                                                                                          <  1:125520/251‑1 (MQ=255)
                                                                                                                                 cgAAATCGGCTTCGCTCAGATTTTGATAGGGAAACTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCAGCAGCTGCCGCATCATCTGCCCGCGTTGCTGAGGCGCTTTTGCCATCGCCCGCCGCACCGATAAATACGCCGAGTTACAGTGCTGACGGGGCGTTGGCgcgg                                                                                                                      >  2:179458/1‑248 (MQ=255)
                                                                                                                                                tCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAc                                                                                                       >  1:308233/1‑251 (MQ=255)
                                                                                                                                                    cttttGATCGGTAACCTGATTGACGCTCTGTTTGGCCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGACATCGCCCGCCGCCCCGATAAATACGCCGTGATCCAGTGCTGCCGGGGCGCTCGCGCGGACGGATTCGCTGTACCatc                                                                                                   >  2:533761/2‑249 (MQ=255)
                                                                                                                                                         gATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTg                                                                                              >  1:438972/1‑251 (MQ=255)
                                                                                                                                                          aTAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgc                                                                                             >  1:115671/1‑251 (MQ=255)
                                                                                                                                                           tAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCAGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTg                                                                                              >  1:5263/1‑249 (MQ=255)
                                                                                                                                                             ggTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATa                                                                                                                                                                                                              >  1:252280/1‑135 (MQ=255)
                                                                                                                                                                aaCCTTTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCACCGTCCAGGCCCTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCcgcg                                                                                                                                                                                                   <  2:346908/143‑1 (MQ=255)
                                                                                                                                                                aaCCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCcgcg                                                                                                                                                                                                   >  1:346908/1‑143 (MQ=255)
                                                                                                                                                                 aCCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCgcagc                                                                                        >  1:176538/1‑249 (MQ=255)
                                                                                                                                                                  ccTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGTGCAGCAgg                                                                                     >  2:46052/1‑251 (MQ=255)
                                                                                                                                                                  ccTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGATGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAgg                                                                                     >  2:266935/1‑251 (MQ=255)
                                                                                                                                                                                 gTTTGGTCAGCTCTCCCGCTAACTGTGCCAACGCCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGGGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGc                                                                        <  1:98340/249‑1 (MQ=255)
                                                                                                                                                                                  ttgggtcaccccggCCGCTAACTTTGTCAACGCCCAGCCCTTACCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGTTGCGCTAATCCCGCGTGCCGCAGCTGCCGCATAATCGGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCa                                                                      <  1:333075/237‑1 (MQ=255)
                                                                                                                                                                                     ggtcaccgctGCCGCTAATTGGGTAAACTCCAGGGCACTATCGCCCGCCGGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCTGGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAAt                                                                  <  1:24437/243‑1 (MQ=255)
                                                                                                                                                                                     ggTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGAGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCCGCAGGGTGTAAACGATGCCAGAAt                                                                  >  1:371971/1‑251 (MQ=255)
                                                                                                                                                                                              tGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATccgccagc                                                          >  2:95809/1‑250 (MQ=255)
                                                                                                                                                                                              ggCCGCAAACGGTGTCACCGTCCAGGCACTACCGCCCGCCAGTTCGTAACCCTTGTCTTCGTGACCGGCTTCGCTAATCCCGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCa                                                         <  1:95809/250‑1 (MQ=255)
                                                                                                                                                                                                  gCTAACTGTGTCAACGTCCAGTCACTATCGCCCGCCAGTTCGTAAACCTTGTATTCGTGACCGGCTTCGCTAATCCCGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATCCCCCGTGTTCCCATGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTCCCCGCCGTTGGGCAGCAGGGTGTGAACGCTGCCAGAATCAGCACGCAtttt                                                     >  2:330271/1‑251 (MQ=255)
                                                                                                                                                                                                  gCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGAGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCACAAGGGTGTAAACGATGCCAGAATAAGCCCGAAtttt                                                     >  2:308705/1‑251 (MQ=255)
                                                                                                                                                                                                    tAACTGTGTCACCGTCCAGGCACTACCGCCCGCCAGTTCGTAACCCTTGTCTTCGGGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAcc                                                                                                      <  1:27864/200‑1 (MQ=255)
                                                                                                                                                                                                    tAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTtct                                                   >  2:221497/1‑251 (MQ=255)
                                                                                                                                                                                                    tAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGGTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAcc                                                                                                      >  2:27864/1‑200 (MQ=255)
                                                                                                                                                                                                        tgtgTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATg                                                                         >  2:76763/1‑225 (MQ=255)
                                                                                                                                                                                                          tgtCAACGTCCACGCACTATCGCCCGCCCGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTTCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCAGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCTCTGCACCACCCGTTGCGTCGCAGGGTGTAAACGATGCCCGAACCAGCCAGACTTTTCTCCctct                                             >  2:191498/1‑251 (MQ=255)
                                                                                                                                                                                                             cAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATaa                                                                                                                                                            >  1:131162/1‑137 (MQ=255)
                                                                                                                                                                                                                aCGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGAtg                                        >  1:261040/1‑251 (MQ=255)
                                                                                                                                                                                                                              aTCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGCTAAATCCGCCGTGTTCCAGTGCTGCCGGCGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTTCGCAGCAGGGCGTAAAAGATGCCAGAATCAGCAAGCCTTTTTTCCGTCTGGATGTGCTCATCGGCg                            >  2:323488/1‑250 (MQ=255)
                                                                                                                                                                                                                              aTCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTCCCAGCCGTTGAGAAGCAGGGAGTAAACGGTGCAGGAATCAACCAGCATTTTCTCCGTCTCGAtgtgaccctcggcgg                           >  2:572015/1‑240 (MQ=255)
                                                                                                                                                                                                                                 gcccgccAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACTCGTGCCGCAGCTACCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGc                        >  1:69362/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ccgccAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCgg                     >  1:373868/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ttCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCTCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCAGGGGCGCTGGCGCGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGAGAAAACGATGCCAGAATCAGCCAGCACTTTCTCCGTCTCGTTGTGCTACTCGGCTAGGCCGAGCgggg               >  2:434342/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    tggCCTTCGTGCCCGGCTACGCTAATCACGCGTGCCGCAGCTGCCGCATAACCTGCCCGCGTTGCTGAGGCGATTTTGCCTGCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTAt         <  1:21859/245‑1 (MQ=255)
                                                                                                                                                                                                                                                       ccTTCGTGTCCTCCTTCGCTAATCCCGCTGCCCGACGCTCCCGCAGAATCTGCCCTCGTTTCTGTGGCTTTTTTCCCATCGCCCGCCGCCCCGATAAATACGCCGTTTTCCAGTTCTGCCGGGGCGCTGGCGCGGTAGTTTGCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCATAATCAGCCAGCATTTTCTCCGTCTCGATGTTCTCATCGGCGAGGCCGATCGGGGAGGTATCTGCAtt  <  1:394249/251‑2 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
TGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGCCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATG  >  NC_000913/4433215‑4433709

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