Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,8500CT26.2% 83.2 / 12.9 42E40K (GAA→AAA) hdfRDNA‑binding transcriptional dual regulator HdfR
Reads supporting (aligned to +/- strand):  ref base C (9/22);  new base T (11/0);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 3.92e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGAAAGCTT  >  NC_000913/3947606‑3948099
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
aaGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGtt                                                                                                                                                                                                                                                         <  1:319792/250‑1 (MQ=255)
aaGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑AGAAGGGTTACACCCCGGTGATTTTccctttg                                                                                                                                                                                                                                                        >  1:21275/1‑246 (MQ=255)
                   gATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑AAAAGGGTCACACCCAGGTGATTTTTCAGTTTTCTGTTTCGAAAACTCAcc                                                                                                                                                                                                                                     >  1:208920/1‑250 (MQ=255)
                             cAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAAGTTCACACCCCGCTGAGTTTTCAGTTTTCTTAttactaaactacctgccgaatggg                                                                                                                                                                                                                           >  1:227114/1‑229 (MQ=255)
                              aaGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAAGATTACACCCCGCTGATTTTTCCGTTGTCTGATTCGCAAACTCACTGCTTACTgggg                                                                                                                                                                                                                          >  1:216396/1‑250 (MQ=255)
                                    cccAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGTGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGTTTCGAAAGCTCACTGCTGACTGGGTCAGAt                                                                                                                                                                                                                     <  1:378018/250‑1 (MQ=255)
                                    cacAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGAt                                                                                                                                                                                                                     <  1:231463/248‑1 (MQ=255)
                                          cACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑AAAAAGGTCACACCCCGGTGTTTTTCCCGTTGTCTGGTTCTAAAACCCCCTGCTGCCTgggggccgtagagggg                                                                                                                                                                                                              >  1:177560/1‑238 (MQ=255)
                                            cTGATTAAGCATACATTACCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑AAAAGGGTCACACCCAGTTGATTTTCCAGTTGCCTGATTAGAAAGCTCATTGCTGCCAGGGTCAGAAAGAGCGAtt                                                                                                                                                                                                            >  1:94903/1‑251 (MQ=255)
                                                  aagcagtcATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCATCGGTTGAACGGATATTGTTTCTGGGGCGTGTT‑AAAAGGTTCACACCCAGTTGATTTGCCAGTTGTCTGATTCGAAATCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGcc                                                                                                                                                                                                       <  1:145450/243‑1 (MQ=255)
                                                     cATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTc                                                                                                                                                                                                   <  1:100948/251‑1 (MQ=255)
                                                       tACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑GAAAGGGTCACACCCCGTTGATTTTTCAGTTGTCTGATTCGAAAACTCCCTGCTGACTGGGGCAGATAGAGCGATTCcgccgcccgc                                                                                                                                                                                                 >  1:263348/1‑247 (MQ=255)
                                                                acaACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCg                                                                                                                                                                                                  <  2:330313/241‑1 (MQ=255)
                                                                acaACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGGGAAAAGGT‑CACACCCCGTTGATTTTTCCGTTTTCTGATTCGTAAACTCCCTGCTGACTGGGGCACATAGAGCGATTCAGCCGCCCGACCaaactgg                                                                                                                                                                                        >  2:304010/1‑247 (MQ=255)
                                                                acaACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑AAAAAGTTCACACCCACGTGATTTTCCAGGTTTCTGATTCGAAAGCTCACTGCTGGCTGGGTCAGCTAGAGGGATTCAGCCGCTCg                                                                                                                                                                                                  >  1:330313/1‑241 (MQ=255)
                                                                      gggCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGCGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGtt                                                                                                                                                                                  <  2:206227/250‑1 (MQ=255)
                                                                        gcgcTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑AAAAAGGTCACACCCCGGTGATTTTTCAGTTGTCTGATTTGAAAACTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCg                                                                                                                                                                                >  1:201377/1‑251 (MQ=255)
                                                                                aCCGATAGAAAACTCGTTATGTCGTGAGGTATGCGTCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑Ctt                                                                                                                                                                       <  2:197331/251‑1 (MQ=255)
                                                                                    aTAGAAAACTCGTTATGTCGTGAGGTATGCGCCTACTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTGTCTGCATAAGGCAGTAGTTTTTCATCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGAGTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAg                                                                                                                                                                   <  2:245707/251‑1 (MQ=255)
                                                                                      aGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑a                                                                                                                                                                <  1:224740/251‑1 (MQ=255)
                                                                                              cGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATAATGTTTCTGTGGCGGGTG‑AAAAGGTTCACCCCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CGTCCAGG‑AAAGtttt                                                                                                                                                         <  2:368174/250‑1 (MQ=255)
                                                                                                      cGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAAtt                                                                                                                                                <  1:133834/250‑1 (MQ=255)
                                                                                                                    cACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTGGAAAAGGTT‑ACACCCCGTTGGTTTTTCCGTTGTTTTGTTCTAAAGCTCCCTGCTTTCTGGGGGAGATAAAGCGGTTTAGCCGCTCGCCCCAAGTGACGCGGTCGGCTAAACTCCCAGG‑AAGGTTTTT‑A‑CAATTCCGTATCCACCGTGc                                                                                                                                 >  2:164301/1‑251 (MQ=255)
                                                                                                                               gggCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCaaaa                                                                                                                      <  2:263348/251‑1 (MQ=255)
                                                                                                                                ggcggcCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAat                                                                                                                     <  1:380173/251‑1 (MQ=255)
                                                                                                                                           aGGTGCTCATGCGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTGGAAAAGGGTAAAACCCACTGCATTTTCC‑GTTGTCTGATTCGGAAAATCACAGCTGTCCGGGTGCGATAAAGAGATTCAGCCCCCCCACCACAGTGACGCGTTCGGCTAA‑ATTCCCGGCAAATTTTTT‑AACA‑TTCCGTCACCACCGATCCCTCCGCCAAATTAttttttgtt                                                                                                          >  2:393398/1‑245 (MQ=255)
                                                                                                                                                            tCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAAtgtttt                                                                                          <  1:334086/250‑1 (MQ=255)
                                                                                                                                                               gCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAAtgttttgt                                                                                        <  1:64745/249‑1 (MQ=255)
                                                                                                                                                                         gtagtTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑GAAAAGTTGCCACCCAGTTGATTTTTCAAGTGTGTTATTTCAAAGCTCACTGCTGACTGGGGCAGATAGAGGGATTTAGCCGCTCCACCACAATGAAGCGCTCGGCTAA‑CTTCCAGG‑GAAGTTTTTTAACA‑TTTCGGAACCACAGGGCTCTCCCCAAAATTATTTTTCGTTATGATTTATAAGTTTTGTTTTACAatcttt                                                                            >  1:416836/1‑249 (MQ=255)
                                                                                                                                                                           agtTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑GAAAGGGTCCCACCCCGATGAATTTTCAGTTTTCTGATTCGGAAGACAACCACTGACTGGGGCAGATAAAGCGATTCAGACGCGCGACCCAAGTGACGCGTTCGGCTAA‑CTTTCCGG‑AAAGTTTTT‑AACAATTTCGTGTCCACAGTGCTCTCCCCAACAATATTTGTTGTTATGATTTTAATgttttgttttttatcctgtc                                                                           >  2:213777/1‑240 (MQ=255)
                                                                                                                                                                            gtTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGGTCACACCCAGATTATTTTTCCGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCCCGACCCAAGTGACGCGGTCGGCTAA‑CTTCCCAG‑AAAGGTTTT‑AACAAATTCGTATCCACAGGGCTCTCCGCAAAATTATTTGTCGTTATGATTTTAATGTTTTTTTTTACACTCTGTCa                                                                          >  1:17710/1‑250 (MQ=255)
                                                                                                                                                                                 tCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGt                                                                     <  1:204417/250‑1 (MQ=255)
                                                                                                                                                                                   aCCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAAc                                                                  <  1:411080/251‑1 (MQ=255)
                                                                                                                                                                                     cGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑AAAAAGGTCACACCCAGTTGATTTTTCCGTTGGCTGATTCGAAAAATCACTGCTGACTGGGTCCGAAAGAGCGAGTCAGCCGCCCGACCAAAGAGACGCGGTCGGCTTA‑CTTCCAGG‑GAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTTTGATTTTAATGTTTTGTTTTTCAATCTGTCAAGCGTaaat                                                                 >  1:272038/1‑248 (MQ=255)
                                                                                                                                                                                                                      gtgTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑GAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGa                                >  2:7648/3‑250 (MQ=255)
                                                                                                                                                                                                                            aaaaGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAgg                          >  2:323011/1‑251 (MQ=255)
                                                                                                                                                                                                                             aaaGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGa                         <  1:164301/251‑1 (MQ=255)
                                                                                                                                                                                                                             aaaGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGa                         <  2:208920/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cacCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAgg                >  1:26478/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cccAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGAt              <  2:26478/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCg          >  2:235938/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGa        <  2:201377/250‑1 (MQ=255)
                                                                                                                                                                                                                                                aTTTTCCAGTTGTCTGATTCGAAATCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGac       >  1:94736/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     ccAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGAAAGCtt  <  1:283587/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
AAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAA‑CTTCCAGG‑AAAGTTTTT‑AACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGAAAGCTT  >  NC_000913/3947606‑3948099

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