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,5010GA94.4% 164.2 / NA 56G378D (GGC→GAC) melBmelibiose:H(+)/Na(+)/Li(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (24/27);  minor base C (3/0);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCATTTTATGGGCCGGAGCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGGCGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCC  >  NC_000913/4344255‑4344736
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
aatagatgctccgcGGAGCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGt                                                                                                                                                                                                                                         >  1:147229/14‑251 (MQ=255)
        tGGGCCGGAGCATCTATTCTTCCGGTGATAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGGATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGATCAGCCtt                                                                                                                                                                                                                                 >  1:94759/1‑251 (MQ=255)
           gCCGGAGCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCATCATTGTGATTGCGGGTATTTTGCTGAATGTGGAAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGc                                                                                                                                                                                                                              >  2:178239/1‑251 (MQ=255)
                      tATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATAGTTGATTACGGTGAATATAAACTGAACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGCCGGTTCAGCCTTTGCGGCTTTTTTct                                                                                                                                                                                                                   >  2:373356/1‑250 (MQ=255)
                            tcCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGc                                                                                                                                                                                                               <  1:79372/249‑1 (MQ=255)
                             ccGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGCCGGTTCAGACTTTGCGGCTTTTTTCATTGCGGt                                                                                                                                                                                                            >  2:144455/1‑251 (MQ=255)
                                   ttAAGCTGTGGTGTTCTCCTGTTAATGGCCTTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCGGACTATGGTGGTGAAGGACGGTTCAGCTTTTGCGGCTTTTTTCATTGCGGTTGGGtt                                                                                                                                                                                                      >  2:90671/1‑251 (MQ=255)
                                           gggtgTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGgatga                                                                                                                                                                                              <  1:90671/250‑1 (MQ=255)
                                                   tcCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGCCGGTTCAGCCTTTGCGGCTTTTTTCa                                                                                                                                                                                                                   >  2:108547/1‑222 (MQ=255)
                                                   tcCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCa                                                                                                                                                                                                                   <  1:108547/222‑1 (MQ=255)
                                                       gTTAATGGCATTAATGAGCTATCTCAACGTCGTCTTCATTGTGATTGCGGGTATTTTGCTGAAGGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTa                                                                                                                                                                                   <  2:287822/250‑1 (MQ=255)
                                                       gTTAAGGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTAc                                                                                                                                                                                  <  2:147229/251‑1 (MQ=255)
                                                                   aaTGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTTGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAAc                                                                                                                                                                      <  1:356495/251‑1 (MQ=255)
                                                                           aTCACAACGGCGGCCTCATTGTGATTGCGGGTATTTTGCTTAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTAc                                                                                                                                                               <  2:125673/250‑1 (MQ=255)
                                                                                  cgtcgtCCTCATTGTGATTGCGGGTCTTTTGATGAATGTGGGAACGGCGCTTTTCTGGGTATTGCAGGTCATCATGGCGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTCCTCCGTGCAGACTATGGGGGGGCAGGACGGTTCAGCCTTTGCGGCTTTTTTAATTGCGGTTGTGTTAGGGGGGATTGGCTATGTACCGGATGtttaaaggtttaagcgagcac                                                                                                                                                       >  2:242831/1‑232 (MQ=255)
                                                                                     cgtcCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGcctatt                                                                                                                                                    >  2:360299/1‑247 (MQ=255)
                                                                                                 gATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGcc                                                                                                                                                        >  1:102992/1‑235 (MQ=255)
                                                                                                 gATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGcc                                                                                                                                                        <  2:102992/235‑1 (MQ=255)
                                                                                                      cGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTg                                                                                                                                                                                                                               >  1:4310/1‑159 (MQ=255)
                                                                                                      cGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTg                                                                                                                                                                                                                               <  2:4310/159‑1 (MQ=255)
                                                                                                       gggTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTttatt                                                                                                                                   >  1:477975/1‑250 (MQ=255)
                                                                                                         gTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGCAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATg                                                                                                                                >  1:189718/1‑251 (MQ=255)
                                                                                                           aTTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGa                                                                                                                               <  1:405654/250‑1 (MQ=255)
                                                                                                                gCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATCGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGctc                                                                                                                         <  1:182455/251‑1 (MQ=255)
                                                                                                                   gAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGGGAATATAAACGGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTAGGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATGATGTATGCAGTTTATTATGATTGctct                                                                                                                        <  2:75366/249‑1 (MQ=255)
                                                                                                                       gggggAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTTATTACTGTGAATATTAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTCTGCGGCTTTTTTCATTTCGGTTGTGTTAGGGTTGATTGGCTATGTACCGAATGTTGAGCAGTCTACGCAAGCCCTAGTAGGTATGCAGTTTATTATGATTGCGCTGCCaa                                                                                                                   <  2:270425/248‑1 (MQ=255)
                                                                                                                        tgGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACAAAct                                                                                                                 >  1:491858/1‑251 (MQ=255)
                                                                                                                          ggAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAActct                                                                                                               >  1:381261/1‑251 (MQ=255)
                                                                                                                               ggCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCttt                                                                                                                                                                                                                                <  1:121178/133‑1 (MQ=255)
                                                                                                                               ggCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCttt                                                                                                                                                                                                                                >  2:121178/1‑133 (MQ=255)
                                                                                                                                gcgtttttttGGGTATTACAGGTCATCATGGTGGCAGTAATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGATCAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTTCCAACTCTGtttt                                                                                                          <  2:392440/246‑1 (MQ=255)
                                                                                                                                           ggTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGCGAATATAAAATGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGACTTTTTTAATTGCGATTGAGTTAGGGATGATTGGCTATGTACCGTATGTTGAACAGACGACGCAAGCCCTATTAAGTATGCTGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTCACGc                                                                                              >  1:432086/1‑251 (MQ=255)
                                                                                                                                             tATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTg                                                                                            >  1:268381/1‑251 (MQ=255)
                                                                                                                                                                       atCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTGCCGTTTCTATCGCCt                                                                  <  1:373356/251‑1 (MQ=255)
                                                                                                                                                                       atCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTAATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCt                                                                  >  1:489641/1‑251 (MQ=255)
                                                                                                                                                                         cGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTc                                                                 <  1:83162/250‑1 (MQ=255)
                                                                                                                                                                             gATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATg                                                             <  1:152254/250‑1 (MQ=255)
                                                                                                                                                                                tACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCGTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTTAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTg                                                          <  1:434685/250‑1 (MQ=255)
                                                                                                                                                                                  cGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGaca                                                       <  1:292980/251‑1 (MQ=255)
                                                                                                                                                                                       aaTATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTAt                                                                       >  1:372756/1‑230 (MQ=255)
                                                                                                                                                                                       aaTATAAACTGCACGTACGCTGTGAAAGTATCGCTTACGCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTAt                                                                       <  2:372756/230‑1 (MQ=255)
                                                                                                                                                                                            aaaCTGCACGTACGCTGTGAAAGGTTCTCTTACTCCGTGCAGACTATGGTGGTGAAGGTCGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGAGGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAg                                             <  2:381261/251‑1 (MQ=255)
                                                                                                                                                                                                 gCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGgatcca                                       <  1:360299/251‑1 (MQ=255)
                                                                                                                                                                                                      tACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGAAGgatcc                                        >  2:444231/1‑246 (MQ=255)
                                                                                                                                                                                                         gCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGGTGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCCGCGCCGGATCCAGATCCAt                                >  2:356953/1‑251 (MQ=255)
                                                                                                                                                                                                         gCTGTGAAAGTATCGCTTACTCCGGGCAGACTATTCTGTTGAAGGGCGGTTCAGCCTGTGCGGCTTTTTTCATTGCGGTGGTGTTAGGGATGATTGCCTATGTACCTAATCTTGAACAGTCTTCGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGGTTTTCATTGTAACGCTGGTTCTCTTCTTCCGTTTCTATCGCCTCAATGGTGCCACGCTGCGCAGGATCCAGATCCAt                                <  1:242831/251‑1 (MQ=255)
                                                                                                                                                                                                           tgtgAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCAttt                              >  2:183835/1‑249 (MQ=255)
                                                                                                                                                                                                                 aGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTgg                         >  2:458619/1‑250 (MQ=255)
                                                                                                                                                                                                                    gtcgttgtctccgtgtgggttATGGTGGTGAAGGACGGTCCAGGCTTTGGGGCTTTTTACATTGCGGTTGAGTTAGGGATGATTGGCTATGTACCGAATGTTTAACAGTCTACGCAAGCTCTAATACGGATGCAGTTTATTATGATTGCTCTTCCAAAGCTGGTTTTCATGGTAACGCTGATTCTTTACTGCCTTGTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCAACTGCGGGATa                      <  2:51216/230‑1 (MQ=255)
                                                                                                                                                                                                                        cctACTCCGTGCAGACTAGGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTCTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAatat                 <  2:491858/249‑1 (MQ=255)
                                                                                                                                                                                                                         ttACTCCGTGCAGACTATGGTGGTGACGGAAGGTTCCGCCTTTGCGCCTTTTTTCCTTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGCACCGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTTTGATTGCTCTACCAACTCTGTTTTTCCTGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCACTGGTGACCAGCTGCGCAGGATCCAGATCCTTCTGCTggcaaattatc                >  2:302221/1‑242 (MQ=255)
                                                                                                                                                                                                                              ccGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCaa            <  2:171203/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTGTTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTa        <  2:394670/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTa        <  1:183835/250‑1 (MQ=255)
                                                                                                                                                                                                                                   cAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGAGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGGATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGAAcc      >  1:404356/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cTATGGTGGTGACGGACGGTTCACCCTTTGCGGCTTTTTTCATTGCGGTTGGGTTAGGGATGATTGGCTATGTACCGAATATTGAACAGTCTACGCAAGCCATATTAGGTATGAAGTTTATTATGATTACTCTACCAACTATGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACACTGCGCAGGATCAAGTTCCATCTGCTGGATAACTATCGCAAGGTaacgacc  >  2:349885/1‑245 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
GCATTTTATGGGCCGGAGCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCGCTTACTCCGTGCAGACTATGGTGGTGAAGGGCGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCC  >  NC_000913/4344255‑4344736

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