Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,025,601 (T)5→4 coding (1251/1401 nt) xanQ → xanthine:H(+) symporter XanQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,025,6010T.94.6% 149.6 / 4.2 37coding (1251/1401 nt)xanQxanthine:H(+) symporter XanQ
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (10/25);  total (11/26)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

CACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAG  >  NC_000913/3025357‑3025844
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
cACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAatat                                                                                                                                                                                                                                               <  2:807390/251‑1 (MQ=255)
    gtATGGCGGGCGAACCATCGCGGTAATGCTGGTTAGCCGCGGCTTATTTCCGATGATTGGCGGCTTCTGCACGACCATTCCCTCGGCAGGTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCAGCGCAGGGATTCGCATCAACATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGcc                                                                                                                                                                                                                                           <  2:593421/250‑1 (MQ=255)
        ggcGGGCGAACCGTCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGc                                                                                                                                                                                                                                        <  2:713892/248‑1 (MQ=255)
            ggCGAACCATTGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCGTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTa                                                                                                                                                                                                                                   <  1:159151/251‑1 (MQ=255)
                      cgcgGGAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTAt                                                                                                                                                                                                                          <  2:791766/250‑1 (MQ=255)
                       gcgGTTATGCTGTTTATCCTCGGCTGATTTCCGATGATTGGCGGCTTCTTCACGTCCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTa                                                                                                                                                                                                                        <  2:771996/251‑1 (MQ=255)
                       gcgGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGTCCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTa                                                                                                                                                                                                                        <  2:139121/251‑1 (MQ=255)
                             aggCTGTTTATCCTCGGCGTATTTCCGATGATTGGCGGCTTCTTCTCGACCAGTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTTATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCAGCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATGTATGTATTAGTTGaa                                                                                                                                                                                                                  <  2:61358/249‑1 (MQ=255)
                                  ggTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCCTTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGACAGCCTCTATTTATGTATTAGTTGAAAAccc                                                                                                                                                                                                             >  2:629173/1‑251 (MQ=255)
                                            ggCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCt                                                                                                                                                                                                   >  1:466264/1‑251 (MQ=255)
                                              cgtATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTgg                                                                                                                                                                                                 <  1:797021/249‑1 (MQ=255)
                                                tagttCCGATGATTGGCGGCTGCTCCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGc                                                                                                                                                                                                <  2:28831/247‑1 (MQ=255)
                                                        atgatTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACt                                                                                                                                                                                      <  1:560365/251‑1 (MQ=255)
                                                                cggcTTCTTCACGACCATTCCCTCGGCCGTTCTGGGAGGCGCACTGACGTTGATGTTTTCCCTGATTGCCATCGCAGGGATTCGCATCATCATCCCCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGCTCCCGAAATTTT‑AAATTATTGCCCGCCTCTATTTATGTATTAGTTGACAACCCTATTTGTGCTGGCGGG‑tttatctcgatttt                                                                                                                                                                              >  2:53733/1‑239 (MQ=255)
                                                                  gcTTCTTCACGACCATGCCCTCGGCAGGTCTGGGAGGCGCAATGACGGTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACTCTTATTTTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACTATCACGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTttat                                                                                                                                                                            <  2:285975/251‑1 (MQ=255)
                                                                      cttcACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAACTATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGGtttttttta                                                                                                                                                                         >  1:269994/1‑245 (MQ=255)
                                                                          aCGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTatattgggtttttttaaatatt                                                                                                                                                                    >  2:667614/1‑230 (MQ=255)
                                                                                   cccccGGCAGGTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCt                                                                                                                                                           <  2:771282/247‑1 (MQ=255)
                                                                                       cGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGGTTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑GTTACTTCGATTTTTTTaaaatatttctttcctggtg                                                                                                                                                       >  2:551469/1‑234 (MQ=255)
                                                                                                          aaTGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCATCGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTTAGGTTTTATTATATTTTTTCCTcccctgggggtcaccacaggaaaaagg                                                                                                                                    >  2:643608/1‑227 (MQ=255)
                                                                                                          aaTGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCAATTTGTGCTGGCGGA‑TTTATTGCGATTTTTTTAAAAATTATTCTCCCCGGTGGGTACCGCCCGGAAAACGt                                                                                                                                    >  2:28913/1‑251 (MQ=255)
                                                                                                                gTTGATGTTTTCCATGATTGCCATCGCAGGGCTTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGGTTAA‑TGCTATTTTTTTGAATATTTTTCTTCCTGGTTGCtttcatcagggatacattttgc                                                                                                                               >  2:346889/1‑230 (MQ=255)
                                                                                                                      ggttacactGATTGCCGTCGCAGGGATTCGGATCATCAGCCCCAACTGTTTAAAGCGCCGTTAAACTCTTATTGTCGCCACTCCTTTAGTTTTAGGGCTTGGCGTATCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCGCTATTTATGTATTAGTTGAAAACCCTATTGGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGGTCTGCCTGGtc                                                                                                                         <  2:494808/242‑2 (MQ=255)
                                                                                                                               gATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCt                                                                                                                                                                                                            <  1:681415/159‑1 (MQ=255)
                                                                                                                               gATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCt                                                                                                                                                                                                            >  2:681415/1‑159 (MQ=255)
                                                                                                                                        cGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGAt                                                                                                      <  1:2768/251‑1 (MQ=255)
                                                                                                                                                   gCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGt                                                                                           <  1:629173/251‑1 (MQ=255)
                                                                                                                                                       catcatcaCCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTaaa                                                                                        <  1:519736/250‑1 (MQ=255)
                                                                                                                                                          catcaCCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGa                                                                                     <  1:346889/251‑1 (MQ=255)
                                                                                                                                                                cAACGGTTTAAAGCGCCGTTAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCATCGGAAGAGATGGATTAACAGTAAAGGAGTCAatg                                                                              <  2:258515/251‑1 (MQ=255)
                                                                                                                                                                                cGTGAAACTCTTATTGTCGCCACTTATTTAGGTTTAGGGTTTGGCTTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTGGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGAATTTGTTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGAAGGATGAACAGTAAAGGAGGCAATGATGGCAGGAGAacaca                                                              <  2:338886/251‑1 (MQ=255)
                                                                                                                                                                                cGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAacaca                                                              <  1:667614/251‑1 (MQ=255)
                                                                                                                                                                                    aaacccctatggTCGCCACTTATTTAGGTTTAGGGCTTTGCGTTGCCGTCGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTAGTTATGTATTAGTTGAAATCCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTTGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGTGTCAATGAGGTCAGGAGAACCCACGtt                                                          <  2:328415/240‑1 (MQ=255)
                                                                                                                                                                                       cacTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTaaa                                                       <  1:46917/251‑1 (MQ=255)
                                                                                                                                                                                              tGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCGACGATCCCGAAATTTT‑AAAAGATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGt                                                 <  1:53733/250‑1 (MQ=255)
                                                                                                                                                                                                 cGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAGATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTTTGCTGGCGGG‑TTAACTTCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACCGGAAAAAGTTCTGCCTGGTATTTCCTCAGCCGGAGGGAATGATTAACAGGAAAGGAGTCAATTCTGTCTGGAGGACACACGGTaaaacaggaacg                                              >  1:539087/1‑242 (MQ=255)
                                                                                                                                                                                                                      gCTTTGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCAcc                        <  2:171876/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gaaattttAAAATATTGCCAGCCTCTATTTATGTATTAGTGGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTCCCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAag  <  2:616988/247‑1 (MQ=255)
                                                                                                                                                                                                                                                attttAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAAATTTGATTTTTTTTAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTTTGCCTGGTATTTCCCCCGCGGAAGAGATGGGTTAACAGTAAAGGAGTCACTGATGTCAGGGGAACACACCTTAAAAGCGGGACGAGGGCGTTTTATTTATTTTCCCCCTACGATCGCTAACACgggag  >  2:53500/2‑245 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
CACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAG  >  NC_000913/3025357‑3025844

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