Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 533,846 T→C 8.6% intergenic (+20/‑70) allR → / → gcl glyoxylate‑inducible transcriptional repressor of all and gcl operons/glyoxylate carboligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913533,8460TC 8.6% 252.3 / 11.9 82intergenic (+20/‑70)allR/gclglyoxylate‑inducible transcriptional repressor of all and gcl operons/glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base T (42/32);  new base C (3/4);  total (46/36)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-01

CCACCTTGCTGAAGGACCTGGAACAAGCGCGTGAACTGGGCTATACCGTAGATAAAGAAGAGCATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCATGTGGAAGGTGCTTCGCACATGGCGGAAGGTTATACCCGCGCAACGGCA  >  NC_000913/533579‑534113
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                           
ccACCTTGCTGAAGGACCTGGAACAAGCGCGTGAACTGGGCTATACCGTAGAGAAAGGAGAGCATGTTGTAGGTCTGAATGGCATAGCTTCAGCAATTTACGATGATGTCGTTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGc                                                                                                                                                                                                                                                                      <  2:505539/275‑1 (MQ=255)
   ccTTGCTGAAGGACCTGGAACAAGCGCGTGAACTGGGCTATACCGTAGATAAAGAAGAGCATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCAAATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCggg                                                                                                                                                                                                                                                                   >  2:424281/1‑275 (MQ=255)
    cTTGCTGAAGGACCTGGAACAAGCGCGTGAACTGGGCTATACCGTAGATAAAGAAGAGCATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCCGCCGCATCCCGCCCTGCGGAGCGGGt                                                                                                                                                                                                                                                                  >  4:158808/1‑275 (MQ=255)
        cTGAAGGACCTGGAACAAGCGCGTGAACTGGGCTATACCGTAGATAAAGAAGAGCATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                                                                                                                                              >  1:134023/1‑275 (MQ=255)
                                      ggCTATACCGTAGATAAAGAAGAGCATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCGACGGGCCTTCATCAAGACTGACAGAAGCTCGTTTTGTAAGTCAGGGTGAGCTGGTCAGAGACAACGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCGGACCTGCGGAGCGGATTTTTTTGACAAAATGTGAAAGTTGGAAAAAtttt                                                                                                                                                                                                                                >  2:75455/1‑275 (MQ=255)
                                                tAGCTAAAGAAGAGCATGTTGTAGGTCTGAATTGAATAGCTTCAGCAATTTACGATGATGTCCGTAGCGCTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGATGGTCAGAGACACAGCAAGCGATATCAGCAAGGCGTTGGGACTGAAAGAACATCATTAATGTCTGTCGCATCCCGCACTGCGGAGCGGGTTTTTTTGACAAAATTGGAAGGTTGGACAACTTTTCCaataaata                                                                                                                                                                                                                      >  4:197333/1‑275 (MQ=255)
                                                tAGATAAAGAAGAGCATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGACAGTTGGAAAAATTTTCCaataaata                                                                                                                                                                                                                      >  2:367208/1‑275 (MQ=255)
                                                      aagaagAGCATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGt                                                                                                                                                                                                                 <  2:158809/274‑1 (MQ=255)
                                                       agaagaGCATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAATTTTGAAAGTTGGAAAAATTTTCCAATGAATAGGGGTAg                                                                                                                                                                                                               >  2:88824/1‑275 (MQ=255)
                                                            agCATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATa                                                                                                                                                                                                          >  4:318101/1‑275 (MQ=255)
                                                              cATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGATTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGAANCCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGGACATCCATAATGTCTGTCGCATCCCGCTCGGCGGACGNGGTTTTTTTGACAAAATTTGAAAGTTGGAANNATTTTCAAATAAATAAAGGGAGGAATaaa                                                                                                                                                                                                        >  3:73839/1‑275 (MQ=255)
                                                              cATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaa                                                                                                                                                                                                        >  1:101852/1‑275 (MQ=255)
                                                                                  cATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGt                                                                                                                                                                                    >  2:13065/1‑275 (MQ=255)
                                                                                       cTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGCAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCCAACATGAGAGCCGTTgcag                                                                                                                                                                               >  2:89708/1‑272 (MQ=255)
                                                                                                 ttaaggttgggTCGGGATTGTGGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGGCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgt                                                                                                                                                                      <  4:367207/264‑1 (MQ=255)
                                                                                                              ggTAGTGTTCTTGCCGCTATCTCCATCTCCGGGACTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTAAGGGTGAGTTGGTCAGAGACACCGCCCGCGATAACAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAAGTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGagaaa                                                                                                                                                        <  1:73839/275‑1 (MQ=255)
                                                                                                                      tgttgCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAAAAAGAAGGTAt                                                                                                                                                >  2:387809/1‑275 (MQ=255)
                                                                                                                                atctccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCg                                                                                                                                                                                                                                                                     <  5:879962/148‑1 (MQ=255)
                                                                                                                                     catctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttt                                                                                                                                                                                                                                                                >  6:60964/1‑148 (MQ=255)
                                                                                                                                      ccttccGGGACTTCATCAAGACTGACAGANTTTTGTTTTGTCAGGCCGGGTGAGCNGGTCCGAAACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATACATAATGTCTGTCGCATCCCGCTCTGCNGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGgt                                                                                                                                <  3:101852/272‑1 (MQ=255)
                                                                                                                                       tctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                                                                                                                                              <  5:60964/148‑1 (MQ=255)
                                                                                                                                                 ttCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggga                                                                                                                                                                                                                                                                  >  6:1095401/1‑133 (MQ=255)
                                                                                                                                                 ttCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggga                                                                                                                                                                                                                                                                  <  5:1095401/134‑2 (MQ=255)
                                                                                                                                                 ttCATCAAGACTACAAGAAGATCGTTTTGTCAGTCAGGGTGGGCTGGTCAGGGACACCTCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGc                                                                                                                     <  4:13065/275‑1 (MQ=255)
                                                                                                                                                  tcatcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTg                                                                                                                                                                                                                                                           >  7:86971/1‑140 (MQ=255)
                                                                                                                                                  tcatcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTg                                                                                                                                                                                                                                                           <  8:86971/140‑1 (MQ=255)
                                                                                                                                                     tcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGc                                                                                                                                                                                                                                                                      >  7:73694/1‑126 (MQ=255)
                                                                                                                                                     tcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGc                                                                                                                                                                                                                                                                      <  8:73694/126‑1 (MQ=255)
                                                                                                                                                      caAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                                                                                                                                             <  6:349978/134‑1 (MQ=255)
                                                                                                                                                      caAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                                                                                                                                             >  5:349978/1‑134 (MQ=255)
                                                                                                                                                      caAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggg                                                                                                                                                                                                                                                                   >  6:409704/1‑128 (MQ=255)
                                                                                                                                                      caAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggg                                                                                                                                                                                                                                                                   >  5:1042741/1‑128 (MQ=255)
                                                                                                                                                      caAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggg                                                                                                                                                                                                                                                                   <  6:1042741/128‑1 (MQ=255)
                                                                                                                                                      caAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggg                                                                                                                                                                                                                                                                   <  5:409704/128‑1 (MQ=255)
                                                                                                                                                          aCTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                                                                                                                                             >  8:283328/1‑130 (MQ=255)
                                                                                                                                                          aCTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                                                                                                                                             <  7:283328/130‑1 (MQ=255)
                                                                                                                                                            tGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGGGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGcaatcaat                                                                                                           <  4:247359/274‑1 (MQ=255)
                                                                                                                                                            cgaaagaagaTCTGTTTGGCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGTCGTTGGGACTGAAAGCACATCAATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGcaatca                                                                                                             <  4:88823/268‑1 (MQ=255)
                                                                                                                                                              aCAGAAGATCGTTTTGTCAGTCAGGGTGAACTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAACAATTTTCCAATAAATAGCGGTAGGAATAAAATGGCAAAATTGAAAGCCGTTGTCGAGGCATTGTAGGTGCTGGAGAAAGAAGGTATCCATACCCCCTTCGGTTTTCCTGCGGCTGCAATCAAAcc                                                                                                         >  4:453563/1‑274 (MQ=255)
                                                                                                                                                               cAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGaa                                                                                                                                                                                                                                               >  5:345077/1‑139 (MQ=255)
                                                                                                                                                               cAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGaa                                                                                                                                                                                                                                               <  6:345077/139‑1 (MQ=255)
                                                                                                                                                                     aTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGt                                                                                                                                                                                                                                            >  8:452080/1‑136 (MQ=255)
                                                                                                                                                                     aTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGt                                                                                                                                                                                                                                            <  7:452080/136‑1 (MQ=255)
                                                                                                                                                                      tCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                                                                                                                                             <  5:813250/118‑1 (MQ=255)
                                                                                                                                                                      tCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                                                                                                                                             >  6:813250/1‑118 (MQ=255)
                                                                                                                                                                      tCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGa                                                                                                                                                                                                                                                          <  7:643629/121‑1 (MQ=255)
                                                                                                                                                                      tCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGa                                                                                                                                                                                                                                                          >  8:643629/1‑121 (MQ=255)
                                                                                                                                                                      tCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                                                                                                                                           <  7:200990/136‑1 (MQ=255)
                                                                                                                                                                      tCGTGTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATACATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAgat                                                                                                                                                                                                                                           >  8:200990/1‑134 (MQ=255)
                                                                                                                                                                          tttGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAg                                                                                                                                                                                                                                             >  6:122150/1‑130 (MQ=255)
                                                                                                                                                                          tttGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAg                                                                                                                                                                                                                                             <  5:122150/130‑1 (MQ=255)
                                                                                                                                                                                   cagGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                                                                                                                                              <  7:733783/104‑1 (MQ=255)
                                                                                                                                                                                   cagGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                                                                                                                                              >  8:733783/1‑104 (MQ=255)
                                                                                                                                                                                       gTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCg                                                                                                                                                                                                                                                                     >  7:573948/1‑93 (MQ=255)
                                                                                                                                                                                       gTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCg                                                                                                                                                                                                                                                                     <  8:573948/93‑1 (MQ=255)
                                                                                                                                                                                              ggTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAACATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGTTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAACGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTCACCAGCCATGCGTTAGCAcaacgttat                                                                        >  4:9067/1‑267 (MQ=255)
                                                                                                                                                                                                tCAGAGACACCGCCCGCGATATCAGCACGGCGTTTGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTa                                                                                                                                                                                                                <  7:666276/137‑1 (MQ=255)
                                                                                                                                                                                                tCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTa                                                                                                                                                                                                                >  8:666276/1‑137 (MQ=255)
                                                                                                                                                                                                          cgcccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCt                                                            >  4:92600/1‑275 (MQ=255)
                                                                                                                                                                                                              cgcgATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCg                                                                                                                                                                                                                                                                     <  6:1093601/70‑1 (MQ=255)
                                                                                                                                                                                                              cgcgATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCg                                                                                                                                                                                                                                                                     >  5:1093601/1‑70 (MQ=255)
                                                                                                                                                                                                                                   ggACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                                                                                                                                             <  6:1083720/57‑1 (MQ=255)
                                                                                                                                                                                                                                   ggACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                                                                                                                                             >  5:1083720/1‑57 (MQ=255)
                                                                                                                                                                                                                                   ggACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCCTGTGGAAGGGTCTTCg                                   >  2:71847/1‑275 (MQ=255)
                                                                                                                                                                                                                                     aCTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCgg                                                                                                                                                                                                                                                                         <  8:601208/43‑1 (MQ=255)
                                                                                                                                                                                                                                     aCTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCgg                                                                                                                                                                                                                                                                         >  7:601208/1‑43 (MQ=255)
                                                                                                                                                                                                                                             cacaTCCATAATGTCGGTCGCAACCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCATGTGGAAGGTGCTTCGCACATGGCgg                         <  4:424280/275‑1 (MQ=255)
                                                                                                                                                                                                                                              acatacatAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCATGTGGAAGGTGCTTCGCACATGGCgg                         <  1:10820/274‑1 (MQ=255)
                                                                                                                                                                                                                                                 tCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaataaata                                                                                                                                                                                                                      <  7:733698/82‑1 (MQ=255)
                                                                                                                                                                                                                                                 tCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaataaata                                                                                                                                                                                                                      >  8:733698/1‑82 (MQ=255)
                                                                                                                                                                                                                                                    aTAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGTGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTa                                                                                                                                                 >  5:917207/1‑148 (MQ=255)
                                                                                                                                                                                                                                                        tgtctgtcGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                                                                                                                                              >  7:391194/1‑35 (MQ=255)
                                                                                                                                                                                                                                                        tgtctgtcGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                                                                                                                                              <  8:391194/35‑1 (MQ=255)
                                                                                                                                                                                                                                                        tgtctgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCTTGTGGAAGGTGCTTCGCACATGGCGGAAGGTTATAcc              >  2:392786/1‑275 (MQ=255)
                                                                                                                                                                                                                                                         gtctgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCaa                                                                                                                                                                         >  7:814932/1‑119 (MQ=255)
                                                                                                                                                                                                                                                         gtctgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCaa                                                                                                                                                                         <  8:814932/119‑1 (MQ=255)
                                                                                                                                                                                                                                                          tctgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCATGTGGAAGGTGCTTCGCACATGGCGGAAGGTTATACCcg            >  2:147348/1‑275 (MQ=255)
                                                                                                                                                                                                                                                            tgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGTGCTGCAATCCATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCCTCTGGAAGGTGCTTCGCCCATGGCGGAAGGTTATACCcgcg          >  2:181634/1‑275 (MQ=255)
                                                                                                                                                                                                                                                            tgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCATGTGGAAGGTGCTTCGCACATGGCGGAAGGTTATACCcgcg          >  4:113416/1‑275 (MQ=255)
                                                                                                                                                                                                                                                            tgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAACTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCATTTGGAAGGTGCTTCTCACATGGCGGAAGGTAATACCcgca          >  1:152789/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                 cATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATTAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCtt                                                                                                                                    <  6:917207/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cccTCTCTGCGGAGCGGGTTTTTTTGACCATATTGGAAATGTGGAAAAATTTTCCCATTAATAGGGGTTGCAAGCAAATGGCACAAATGAGCGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGGGCTGCAGGCAATCCGTTATGCTCAGCGATGCGTAAGCACGGCGGTCTTCGTCACATTCTGGCGCGTCATGTGGAAGGTGCTTCGCACATGGCGGAAGGTTATACCCGCGCAACGGCa  <  4:71846/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                           
CCACCTTGCTGAAGGACCTGGAACAAGCGCGTGAACTGGGCTATACCGTAGATAAAGAAGAGCATGTTGTAGGTCTGAATTGCATAGCTTCAGCAATTTACGATGATGTCGGTAGTGTTGTTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAGCTGCAATCAATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCATGTGGAAGGTGCTTCGCACATGGCGGAAGGTTATACCCGCGCAACGGCA  >  NC_000913/533579‑534113

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: