Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,344,501 G→A G378D (GGC→GAC)  melB → melibiose:H(+)/Na(+)/Li(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,344,5010GA98.1% 150.2 / NA 52G378D (GGC→GAC) melBmelibiose:H(+)/Na(+)/Li(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/27);  minor base T (0/1);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGCCGGAGCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGGCGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCC  >  NC_000913/4344265‑4344741
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
ggCCGGAGCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAACCTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTg                                                                                                                                                                                                                                     >  1:102399/1‑251 (MQ=255)
       gCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCttt                                                                                                                                                                                                                              <  1:19144/251‑1 (MQ=255)
        cATCTATTGTTCCGGTGTTAAGCTGTGGTGTTCTACTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGAACGCTGTGAACGTATAGCTTACTCCGTGCAGACTATGGTGGTGAAGGAAGGTTCAGCCTTTGCGGCtttt                                                                                                                                                                                                                             >  2:300285/1‑251 (MQ=255)
            tatgatcccGGTGGTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATACCAACGTCGGCCTAATGGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTTGGTATTACAGGTCATCATGGTGGCAGATATCGTTGAGTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTAGGGTGGTGAAGGACGTTTCATCCTTTGCGGCTTTTTTc                                                                                                                                                                                                                          <  2:242371/243‑1 (MQ=255)
            tATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTACTGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGAACGCTGTGACAGTATCGCTTACTCCGTGCAGACAATGGTGGTGAAGGACGGTTAGGCCTTTGCGACTTTTTTc                                                                                                                                                                                                                          >  2:201315/1‑250 (MQ=255)
                         ttAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGtt                                                                                                                                                                                                            >  2:176836/1‑251 (MQ=255)
                          tAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGtt                                                                                                                                                                                                            <  2:501332/250‑1 (MQ=255)
                               tgtgGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGgat                                                                                                                                                                                                      >  1:375967/1‑251 (MQ=255)
                                       tctcCTGTTAATGGCATTACTGAGCTATCACACCGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGTATGTGGGAACCTCGCTTTTCTGGGTATTACCTGTCATCATGGTGGCAGATATCGTTGATTCCGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTTCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTTCTGCTTTTTTCATTGCGGTTGGGTTAgtgtatgttgtct                                                                                                                                                                                              >  2:197147/1‑239 (MQ=255)
                                                   ggCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGTCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCt                                                                                                                                                                                                                                <  2:426937/205‑1 (MQ=255)
                                                   ggCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCt                                                                                                                                                                                                                                >  1:426937/1‑205 (MQ=255)
                                                             aGCTATAACAACGTCGTCCTCAGTGTGATTGCGGGTATTGTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTAAGCCGTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTc                                                                                                                                                                        <  1:197147/251‑1 (MQ=255)
                                                                 aTCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACg                                                                                                                                                                    >  1:283846/1‑251 (MQ=255)
                                                                 aTCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACg                                                                                                                                                                    <  2:482669/251‑1 (MQ=255)
                                                                  tCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGTATGTTGAACAGTCTTCGc                                                                                                                                                                   >  1:236793/1‑251 (MQ=255)
                                                                              gttttgtgtgATTGCTGGTATTGTGCGGAATGTGGGAACGGTGGTTGTCTTGGTATTATAGGTCAGCATGGTGGCATATATCGTTGATTACTGTTAATATAAACTGCACTTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTTGTGAAGGTCGGTTCAGCCTTTGCGTCTTTTTTCATCGCGGTTGTGTTAGTGATGATTGGCTATGTACTGAATGTTGAACAGTCTATGCAAGcccgagta                                                                                                                                                        <  2:327272/244‑6 (MQ=255)
                                                                              ccTCATTGTGATTGCGGGTATTTTGCTGAAGGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGGTGAACAGTCTACGCAAGCCCTAT‑TAg                                                                                                                                                      <  2:102399/251‑1 (MQ=255)
                                                                                    tgtgATTGCGGGTATTTTGCTGAATGTGGGACCGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGAATATGGTGGTGAAGGACGGTTCAGCCTTTGCGTCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGTCCTTGTACCGATGGTTGAACAGTCTACCCAAGCCCTGT‑TAggatttc                                                                                                                                                >  2:283574/1‑246 (MQ=255)
                                                                                       gATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGt                                                                                                                                             <  1:298217/251‑1 (MQ=255)
                                                                                          tGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGttt                                                                                                                                           <  1:176836/250‑1 (MQ=255)
                                                                                            cGGGTATTTTGCTGAACGTGGGAACGGCTCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGTCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TTGGTATGCAGTttat                                                                                                                                         <  1:494009/250‑1 (MQ=255)
                                                                                             gggTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTttatt                                                                                                                                        >  1:99791/1‑250 (MQ=255)
                                                                                               gTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATg                                                                                                                                     <  1:597917/251‑1 (MQ=255)
                                                                                                 aTTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGAt                                                                                                                                   >  1:151193/1‑251 (MQ=255)
                                                                                                         gAATTTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATTGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCGAT‑TAGGTATGCAGTTTATTATGATTGCTCTAc                                                                                                                           <  1:546503/251‑1 (MQ=255)
                                                                                                           aTGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGCACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTAcc                                                                                                                          >  2:154042/1‑250 (MQ=255)
                                                                                                           aTGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGAAGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTAcc                                                                                                                          >  1:245936/1‑250 (MQ=255)
                                                                                                                              cTGGGTATTACAGGTCATCATGGTGGCAGATCTCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCCGTTTATTATGCTTGCTCTACCAACTCTGTTTTTCATGGTa                                                                                                       >  1:437532/1‑250 (MQ=255)
                                                                                                                              cTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTa                                                                                                       >  1:354355/1‑250 (MQ=255)
                                                                                                                                gggaaTTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGGGAATATAAAGTGCACGTACGCGGTGAAAGAATAGCTTACTCAGAACAGAATATGGTGGTGAAGGACGAATCAGCATTTGAGGCTTTTTTCATTGCGGTTGTGTTAGGGAGGAGTGGCTATGTACCGAATGTTGAAAAGTCGACGCAAGACC‑ATGTTGGGATGCATGTTCTTATGATTGCTCTAACAACACTGTTTTTccagggaa                                                                                                      >  2:607426/5‑242 (MQ=255)
                                                                                                                                gggTAGTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAAc                                                                                                     <  1:300285/250‑1 (MQ=255)
                                                                                                                                             cagcatGGTGGCAGATAGCGTTGATTACGGCGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGGGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCGTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTAGGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTAGGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGAGTCTCTAc                                                                                        <  2:412989/247‑1 (MQ=255)
                                                                                                                                                            tatCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGcc                                                                        >  1:409640/1‑251 (MQ=255)
                                                                                                                                                             atCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCt                                                                       <  1:201315/251‑1 (MQ=255)
                                                                                                                                                                  tGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAAt                                                                   <  2:409640/250‑1 (MQ=255)
                                                                                                                                                                   gATTTCGGTGAATAGAAACTGCACGTACGCTGTGAAAGTATCGCTGTCTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCGATCGTCTCAATg                                                                  <  2:36381/250‑1 (MQ=255)
                                                                                                                                                                                   caCTGCACGTACGCTGTGACAGTATCACTTACGCCGGGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGGATGTTGCACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGCTTGCTCTAACAACTCTGTTTTTCATTGTAACGCTGATTATCTACTTACGTTTCTATCGCCTCAATGGTGATACGCttggaaa                                                  >  2:358474/2‑244 (MQ=255)
                                                                                                                                                                                    aCTGCACGTACGCTGGGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGAGACTTTTTTCCTTGCGGTTGTGTTAGGGCGGATTGGCTATGAACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TCGGTATGCAGTTTATTATGATTGCTCAACCAACTCTGTTTTTCATGGGAACGCTGATTCTCTACTTCCATTTCTATCGCCTCCATGGTGCCACGCTGCGAAgg                                                 >  2:674658/1‑250 (MQ=255)
                                                                                                                                                                                        cacttACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTAGGGTGTTGAAGGACTGTTCATCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTTCCGAATGTTGAACATTCTACGCAAGTCCTAG‑TAGGTATGCTGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTGCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGgatcca                                            <  2:356067/247‑1 (MQ=255)
                                                                                                                                                                                                    gAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGACTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATctgc                                 <  1:562828/250‑1 (MQ=255)
                                                                                                                                                                                                    gAAAGTACCGCTTACTCCGTGCAGCATTTGGTGGTGAAGGACGGTTCGGCCTTTGCGGGTTTGTGCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TATGTATGCTGTAGATTATGATTGCTTTACCAACGCTGGTTTTCATGGTAACGCTGATTCTCTTCTGCGGTTTCTATCGCCTCAATGTTGACACGCTGCGCAGGATCCAGATCCATctgc                                 <  1:286458/250‑1 (MQ=255)
                                                                                                                                                                                                      aaGTATCGCTTACTCCGTGCAGACTTTGGTGGTGAAGGACGGTTCTGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTGCCGAAAGTTGAACAGTCTACGCAAGCCCTAG‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGTTAACGCTGATTCTCTACTTCCGTGTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTTCTgg                              <  2:245936/251‑1 (MQ=255)
                                                                                                                                                                                                       agtctagCTTACACCGTGAAGACTCTGAAGGTGAAGCACGGTTCAGCCTTTGCGGATTTTTTCATTGCGGTTGTGTTAGGGAGGAGTGGCTATATACCGAATGTTGAAAAGACCACGCCAACTCTAT‑TAGGTATGGAATTTATTAAGATTGCTCTCCCAAATCTGTTTATCCTGGTAACGCTGGATCTATACTTCCGTTTCTATCGCCTCAATGGTGACACTCGGCAAAGGATCCAGCTCAATCTGCTgg                              >  2:450553/7‑250 (MQ=255)
                                                                                                                                                                                                            cGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAat                        >  2:33937/1‑251 (MQ=255)
                                                                                                                                                                                                                tACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAatat                      <  1:33937/249‑1 (MQ=255)
                                                                                                                                                                                                                    ccGTGCAGACTATGGTGGTTAAGGACGGTTCAGCCTTTTCGGATTTTTTTAGTGCGGTTGTGTTAGGGATGATTTGCTATGTACCGTATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATTCAGTTTATTATGATTGCTCTACCAACTCTGTTTGTCAGGGTAACGCTGATTCTCTTCGTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATAGCGCaaa                <  2:666055/251‑1 (MQ=255)
                                                                                                                                                                                                                          aGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCg          <  1:612012/251‑1 (MQ=255)
                                                                                                                                                                                                                               aTGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCCg      <  1:9741/250‑1 (MQ=255)
                                                                                                                                                                                                                                 ggtggtGAAGGACGGTTCAGCCTATGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCCgtg    >  2:294990/1‑248 (MQ=255)
                                                                                                                                                                                                                                  gtggtgAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGcc  <  1:294990/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aaGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCa                  <  1:179354/229‑1 (MQ=255)
                                                                                                                                                                                                                                        aaGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCa                  >  2:179354/1‑229 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
GGCCGGAGCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGGCGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTAT‑TAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCC  >  NC_000913/4344265‑4344741

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