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.90.0% 141.3 / 11.8 40coding (1251/1401 nt)xanQxanthine:H(+) symporter XanQ
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base . (14/22);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGG‑ATTAA‑CAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCC  >  NC_000913/3025353‑3025839
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
gCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGGTGATTGGCGGCTTCTTCACGACCATTCCCTCGGCCTGTCTGGGAGGCGCAATGGCGGTGGTGTTTTCCATGATTGCCATCGCAGGGTTTCGCATCCACATCACCAACGGTTTAAAGCGCCGTGGAACACTTATTGTCGCCCCTTCTTTAGGTTTTGTGCTTTGCGTCTCCTACGGATCC‑GACATTTTTa                                                                                                                                                                                                                                                   >  2:231437/1‑250 (MQ=255)
       gTTATGTCGGGCGCATCCTCGAGGTAATGCTAGTTATCCTCGGCGTATCTTCGGTGATTCGCGGCTTCTTCACGACCACTCCCTCGGCAGTTCTGGGATGCGCAATTACGTTGATGGTTTCCATGATTGCTGTCTCAAGGAATCGTATCATCATCGCCAACGGTTTAAAGCGCCGTGAAACACTTATCGTCGCCATTTCGTTAGGTTTAGGGATTGGCGACTTCTACG‑ACCCCGAAAATTT‑AAAATATTGc                                                                                                                                                                                                                                          <  2:474513/251‑1 (MQ=255)
            gTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATTATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGcc                                                                                                                                                                                                                                     <  1:25352/251‑1 (MQ=255)
                 gCGAACCGTCGCGTGAATGCGGGTTATCCGCGGCGTATTTCCTATGATTGGCGGTGTCTACACTGCCATTCCGTCGGCCGTTCTGGGAGGCGCAATGATGTTGAGGTTTTCCATTATTGCCATCGCAGGGATTCGCATCATCATCACCAACAGTTTAAAGCGCCGTGAAACACTTATTGTCGCCAGTTCTTTAGGTTTAGGGCTTGGCGGCTCCTACG‑ATCCCGAAATTTT‑AAAATATGGTTAGCCTCTAt                                                                                                                                                                                                                                <  2:178110/251‑1 (MQ=255)
                  cGAACCATCGCGTGAATGCTGTTTATCCTCGGCTTATTGCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCCATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCAACATCCCCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTCCG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTAtt                                                                                                                                                                                                                               <  1:484204/251‑1 (MQ=255)
                      ccATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATg                                                                                                                                                                                                                           >  1:434592/1‑251 (MQ=255)
                           gcgGTAATGCTGGTTATCCGCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTa                                                                                                                                                                                                                      <  1:442873/251‑1 (MQ=255)
                            cgGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAg                                                                                                                                                                                                                     <  1:127015/251‑1 (MQ=255)
                                            ccgcGGCTTATTTCCGATGATTGGCGGCTCCTTCACGTCCACTCCCTCTGCAGTTCTGGGAGGCGCAATGACGTGGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCACCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTAtctg                                                                                                                                                                                                     <  2:121874/248‑4 (MQ=255)
                                               cGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTtgt                                                                                                                                                                                                    <  1:94040/249‑1 (MQ=255)
                                                     aTTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCg                                                                                                                                                                                             <  1:358530/250‑1 (MQ=255)
                                                           gatgatTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAAc                                                                                                                                                                                     <  1:180409/251‑1 (MQ=255)
                                                              gatTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTTACTGc                                                                                                                                                                                  >  1:66807/1‑251 (MQ=255)
                                                                            tcccGACCATTCCCTCGGCAGTTCTGGGAGGCTCAATGACGTTGATGTTTTCCATGAGTGCCATCGAAGGGATTCGCATCATTATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCTATTTTATTAAat                                                                                                                                                                     <  2:336602/247‑1 (MQ=255)
                                                                                   cATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑GTATATGCGATTTTATTTAAtaatattct                                                                                                                                                             >  1:478100/1‑248 (MQ=255)
                                                                                             gCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGc                                                                                                                                                  <  1:101631/251‑1 (MQ=255)
                                                                                                               aTGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCt                                                                                                                                                                                                          >  2:322803/1‑179 (MQ=255)
                                                                                                               aTGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCt                                                                                                                                                                                                          <  1:322803/179‑1 (MQ=255)
                                                                                                                 gACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTTc                                                                                                                               <  1:299041/250‑1 (MQ=255)
                                                                                                                            tttCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTTCGCTTTTTTTAAATATTTTTCTCCC‑TGGTGGCTACCGTCAGGAAAACGTTCTTCCTGGGAttt                                                                                                                   >  1:304989/1‑250 (MQ=255)
                                                                                                                             ttCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTa                                                                                                                   <  1:157833/250‑1 (MQ=255)
                                                                                                                                 aTGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCa                                                                                                              <  2:110130/251‑1 (MQ=255)
                                                                                                                                      tGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAAATTAGATTTTTTTTAATAATATCCTCCC‑TGGTGGCTACCGACAGGAAAACGtttttactggtgtttcctccgcggg                                                                                                         >  2:133585/1‑228 (MQ=255)
                                                                                                                                              cAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATgg                                                                                                 <  1:484439/251‑1 (MQ=255)
                                                                                                                                                            atcatcaCCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTTTAAAATATTGCCAg                                                                                                                                                                                                                                       <  2:230756/106‑1 (MQ=255)
                                                                                                                                                            atcatcaCCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTTTAAAATATTGCCAg                                                                                                                                                                                                                                       >  1:230756/1‑106 (MQ=255)
                                                                                                                                                                        ggTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTTACTGCGATTTTATTAAATATTTTTCTCCCCTGG‑GGCTACCGGCAGGAAAAAGTTTTGACTGGTATTACCTCCACGGGAGAGGGGG‑ATTAA‑CCATAAAGGGGACAatgattg                                                                     >  2:359989/1‑249 (MQ=255)
                                                                                                                                                                        ggTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑GTAACTGCGGGTTTTTTTAATATTATTCTCCC‑TGGGGGCTACCCCCAGGAAAACGTTTTGCCTGGGATTACCTCAGCGGAAGAGAGGG‑ATTAA‑AACTAAAGGAGttcatgtttt                                                                     >  1:409127/1‑242 (MQ=255)
                                                                                                                                                                         gTTTACAGCGCCGTGAAACACTTATTTTCTCCACTGCTTTAGGTTTTGGGCTTGTCGTCTCCTACT‑ATCCCGAAATTTT‑AAAATATGTCCAGCCTCTATTTATGTATTAGTTGAAAACTCTATTTGTGCTGGCTGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTGCTGTCTGGTATTACCTCATCGGAAGAGATGG‑ATTAA‑CAGTATAGTAGTCAATGATGTc                                                                    <  2:261598/251‑1 (MQ=255)
                                                                                                                                                                           ttAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGG‑ATTAA‑CAGTAAAGGAGTCAATGATGTCa                                                                   <  1:191116/250‑1 (MQ=255)
                                                                                                                                                                             aaaGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑GTAAATGCGATTTTTTTAAATATTTTCCTCCC‑TGGGGGCTACCGACAGGGAAACGTTCTTCCCGGTATTACCTCCGCGGAAGAGAGGG‑AGTAA‑CAGTAAAGGAGGCAATTATGTCAggg                                                                >  1:60784/1‑250 (MQ=255)
                                                                                                                                                                                 cgcCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTTTGCTGGCGGGGTTAA‑TGCGATTTTATTAAATATTATCCTCCC‑TGGTGGGTTCCGACAGGAAAACGCTCTGCCTGGTATTTCCTCAGCGGAAGAGATGG‑ATTAA‑AAGTAAAGGGGGCAATTATGTCAGgaga                                                              >  1:165049/1‑250 (MQ=255)
                                                                                                                                                                                           cacTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGG‑ATTAA‑CAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTaaa                                                  <  2:66807/251‑1 (MQ=255)
                                                                                                                                                                                                        cACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTTTTTAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGGAGAGATGG‑GTTAA‑CAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGc                                     >  1:160813/1‑251 (MQ=255)
                                                                                                                                                                                                                        gggCTTGGCGTCTCCTACG‑ATCTCGAAATTTT‑AAAAGATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGGTCTGCCTGGTATTACCTCAGCGGAAGAGATGG‑ATTAA‑CAGTAAAGGAGTCAATGATGTCAGTAGAACACACGTTAAAAGCGGTACGAGGCAGTGTTATTGATGTct                     <  2:41383/251‑2 (MQ=255)
                                                                                                                                                                                                                         ggCTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAAACGCGGTTTTATTAAATATTATCCTCCC‑TGGTGGGTACCGACAGGAAAACGTTCTGCCTGGTATTACCCCCGCGGGAGAGATGG‑ATTAAACA‑TAAAGGGGTCAAAGTTGTCAGGGGAACACACGTTAAAAGCGGGACGAGGGAGTTTTATTGTTGTCAc                    >  1:471679/1‑251 (MQ=255)
                                                                                                                                                                                                                           cTTGGCGTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTCTTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGAATTTATTAAATATTATCCTCCC‑TGGGGGCTACCGACAGGGAAACTTTCTGCCTTGTATTACCTCAGCGGAAGAGATGGTATAAA‑C‑GTTAAAGAGACAATGATGTTAGGAGAAAACCCGTTAAAAACGGTAAGAGGGCGTTTTATTGATGTccacc                  >  1:478047/1‑247 (MQ=255)
                                                                                                                                                                                                                                 gtatcCTACG‑ATCCAGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAAATGCGAATTTATTAAATATTTTTCTCCC‑TGGGGGGTACCGACAGGAAAACGCTCTGCCTGGAATTAACTCAGCGGAAGAGATGG‑ATTAA‑CAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTAAGAGGACGTTTTATTGCTGTCAACCGTAcg             >  1:368115/4‑250 (MQ=255)
                                                                                                                                                                                                                                 gTCTCCTACG‑ATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAAATGCGAATTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACCGGAAAACCTTTTTCCTGGTATTACCTCAGCGGGAAAGATGG‑GTTAA‑AAGTAAAGGAGTCAATGATGTCAAGAGAACACACGTTAAAAGAGGTAAGGGGCAGTTTTATTTGTTTCACCCGTAcg             >  1:368117/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aTCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGGTGAAAACCCTATTTGTGCTTGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGG‑ATTAA‑CAGTAAAGTAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATTTCACCCGTACGATCGATAAccc  <  1:463958/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
GCTTCACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACG‑ATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGG‑TTAACTGCGATTTTATTAAATATTATCCTCCC‑TGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGG‑ATTAA‑CAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCC  >  NC_000913/3025353‑3025839

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