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

AACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGAAAGCTT  >  NC_000913/3947600‑3948099
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                            
aaCTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTc                                                                                                                                                                                                                                                              <  2:480344/251‑1 (MQ=255)
    gTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTTATTTTCCAg                                                                                                                                                                                                                                                           <  1:361277/250‑1 (MQ=255)
      aaGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGtt                                                                                                                                                                                                                                                         <  1:303534/250‑1 (MQ=255)
              ttGGGCATCATGATTTTGATCCACGCGTCCCAGCCACTGATTAAGCATTCATTCCCACAACGAGGCGCTGGCGCCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCTCCTCCTTACGGGCGGCCTGCCAGGTGCTCATGATCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGACGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGAtt                                                                                                                                                                                                                                                 <  2:398394/249‑1 (MQ=255)
                 ggCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑AAAAGGTTTACACCC‑AGTTGATTTTTCAGTTGTTTTATTCGa                                                                                                                                                                                                                                              >  1:586392/1‑250 (MQ=255)
                         gATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAAGGTCAAACCC‑AGTTGGTTTTTCAGTTGTCTGATTTGAAAACTCAcc                                                                                                                                                                                                                                     >  1:451232/1‑250 (MQ=255)
                         gATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGGGAAAAGGT‑CACACCC‑AGTTTGTTTTTCAGGTGTCTGATTCTAAAACTCAct                                                                                                                                                                                                                                     >  1:563163/1‑251 (MQ=255)
                          aTTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCActg                                                                                                                                                                                                                                    <  1:576260/251‑1 (MQ=255)
                           ttttGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑AAAAAGGTCACACCC‑AGTTGATTTTTCCGTTGTCTGATTTGGAAGCTCActtc                                                                                                                                                                                                                                   >  1:375675/1‑249 (MQ=255)
                                  aCAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGGTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTg                                                                                                                                                                                                                             <  1:197401/250‑1 (MQ=255)
                                     aGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTGGAAAAGGTT‑ACACCC‑CGTTGATTTTTCAGTTGTCTGATTCGAAAGATCACTGCTGGCTggggg                                                                                                                                                                                                                         >  1:571700/1‑249 (MQ=255)
                                          cccAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGAt                                                                                                                                                                                                                     <  1:181660/250‑1 (MQ=255)
                                                                            aGGCGCTGGCACTTATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCGCCTTACGGGCGTCCTGCCAGGTGCTCATTAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑tt                                                                                                                                                                                 <  2:361129/251‑1 (MQ=255)
                                                                            aGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAATTGACGCG‑tt                                                                                                                                                                                 <  1:398293/251‑1 (MQ=255)
                                                                            aGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGt                                                                                                                                                                                          <  2:215635/243‑1 (MQ=255)
                                                                            aGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGGTCACACCC‑AGTTGATTTTTCAGTTGTCTGATTCGAAAGCTCACCGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCaaaagagaacggg                                                                                                                                                                                  >  1:215635/1‑241 (MQ=255)
                                                                                       ccGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGGTCACACCC‑AGTTGGTTTTTCAGTTGTTTGGTTCGAAAACTCACCGATGGATGGGGCAGAAAGAGGGGTTCAGCCCCTCGAACAAAAGGACGCGGTTCGG‑TAACTTc                                                                                                                                                                      >  2:568435/1‑251 (MQ=255)
                                                                                           tAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAAGGTCACACCC‑AGATGTTTTTTCAGTTGTTTGTTTCGAAAAATCACTGCTGACTGGGTGCGATAGAGGGATTCAGCCGCCCGACCCAAGtaaagcgtgcgggtcactgccacg                                                                                                                                                                   >  2:451433/1‑228 (MQ=255)
                                                                                           gaGAAATCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑Agg                                                                                                                                                                 <  2:62280/250‑1 (MQ=255)
                                                                                            aGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑Agg                                                                                                                                                                 <  1:174283/250‑1 (MQ=255)
                                                                                                   tCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑GAAAAGTTCACACCC‑CGTTGATTTTTCAGTTGTCTGATTCGAAAACTCACTGCTGGCTGGGTCAGATAAAGCGATTCAGCCGCTCGCCCAAAATGACGCG‑TTCGGCTAAATTTC‑AGG‑AAAGattt                                                                                                                                                        >  1:405088/1‑251 (MQ=255)
                                                                                                    cGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGttttt                                                                                                                                                       <  1:535018/251‑1 (MQ=255)
                                                                                                    cGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGtttt                                                                                                                                                        <  1:447532/250‑1 (MQ=255)
                                                                                                        aTGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑GAAAAGTTTACACCC‑ACGTGATTTTTCCGTTGTCTGATTTGAAAGCTCACTGCTGACTGGGGCAGATAGAGCGATTCAGCCCCCCGAACAAAGTGACGCG‑GTTGGCTAACTTCC‑AGG‑AAAGttttttac                                                                                                                                                    >  1:110787/1‑247 (MQ=255)
                                                                                                                               cctTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGctct                                                                                                                              <  2:571700/249‑1 (MQ=255)
                                                                                                                                     gggCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑GAAAGGATTACACCC‑CGTTGGTTTTTTCGTTTTGTTATTTGAAAACCCACTGCTTACTTGGGCAGAGAGAGGGAGTTAGCCCCTCCACCAAAGTGGAGCG‑GTTGGCTAACTTCC‑AGGGACAGTTTTTAAAAAttacatttcacaccgctcctcccaaa                                                                                                                        >  2:510754/1‑226 (MQ=255)
                                                                                                                                               ccAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCTCCGGCAGCTGTTAAACGGATATTGTGTCTGTTGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATGTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCATATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TGCGGCTAACTTCC‑AGT‑AAAGTGTTTAACAATTCCGTATCCACAGTGCTCTCCTCAAAATTATTTGTCg                                                                                                            <  2:520461/251‑1 (MQ=255)
                                                                                                                                                cAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGt                                                                                                           <  1:310001/251‑1 (MQ=255)
                                                                                                                                                   gTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTAt                                                                                                        <  1:407180/251‑1 (MQ=255)
                                                                                                                                                     gCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATg                                                                                                       <  1:431648/250‑1 (MQ=255)
                                                                                                                                                              cGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGGGAAAAGGTTC‑CACCC‑CGTTGATTTTCCCGGTTTCTGGTTTGAAAACTCAATGCTGAATTGGTCAAAAAGAGCGATTCAGCCGCCCGACCAAAATGACGCG‑TTCGGCTAACTTCCCGGG‑AAGGTTTTTAA‑AAATCCCTATTCCCACTGCTCTCCGCCAAAATAttttttgtttttattttaaagt                                                                                             >  1:100962/1‑232 (MQ=255)
                                                                                                                                                                     gCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGGTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGCGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAAtgttttgtt                                                                                       <  2:95254/250‑1 (MQ=255)
                                                                                                                                                                            gCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGATGGTTTTTCAGGTTTTTGATTCGAAAACTCACTGCTGACTGGGGCAGATAGGGCGATTCAGCCGCTCGACCCAAATGGCGCG‑GTCGGCTTACTTCC‑AGG‑GAAGTTTTTAACAATTCCGTATCCACCATGCTCTCCGCAAAATTATTTGTCGTTTTGATTTTAATgttttgtttttcaatc                                                                               >  1:250713/1‑245 (MQ=255)
                                                                                                                                                                              agtagtTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCCCAGT‑GGTTTTTCAGTTTTCTGGTTTGAAAACTCCCTGCCGGCCGGGTCAGATAGAGGGATTTCGCCGCTCGCCCCAAGGGGAGCG‑GTCGGCTAACTTTC‑AGG‑GAAGTTTTTTAAAATTCCCTTTCCCCAGAGCTCTCCCCAAAATTATTTTTTGTTTTGATTTTAATTTTTTTTTTTTCActctt                                                                             >  2:125519/1‑250 (MQ=255)
                                                                                                                                                                                         cccGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATGCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAAc                                                                  <  2:234539/250‑1 (MQ=255)
                                                                                                                                                                                          ccGGCAGCGGTTAAACGGATTTTGTTTTTGTGGCGGGGG‑AAAAGGGTCACTCCC‑CGTTGATTTTCCAGTGGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACt                                                                 <  2:35977/251‑1 (MQ=255)
                                                                                                                                                                                              cAGCGGTTAAACGGATAGTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATa                                                             <  2:254680/251‑1 (MQ=255)
                                                                                                                                                                                                gCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACt                                                           <  2:451232/251‑1 (MQ=255)
                                                                                                                                                                                                   gTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCcg                                                        <  2:34787/251‑1 (MQ=255)
                                                                                                                                                                                                     tAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGtt                                                                                                          <  2:173372/199‑1 (MQ=255)
                                                                                                                                                                                                     tAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAAGTTTACACCC‑CGTTGATTTTTCCGTTGTCTGATTCGAAAGCTCACCGCTGACTGGGGCAGATAGAGCGATTCAGCCGCTCGACCAAAGAGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGGTTTTAACAAATCCGTATCCCCAATGCTCTCCGCAAAATTAttttttgtt                                                                                                          >  1:173372/1‑193 (MQ=255)
                                                                                                                                                                                                         cGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCAt                                                  <  1:304107/251‑1 (MQ=255)
                                                                                                                                                                                                              attttttCTGTGGCGGGTG‑AAAAGGTTCAAACCC‑AGTTGAATTGCCAGTTGTCGGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAATCGTAACTAATCCTCCGCGCCATAACt                                              <  2:370174/246‑1 (MQ=255)
                                                                                                                                                                                                              atTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACt                                              <  1:79798/250‑1 (MQ=255)
                                                                                                                                                                                                                    tCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCg                                        <  2:375675/250‑1 (MQ=255)
                                                                                                                                                                                                                            gggTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGa                                >  2:440680/1‑250 (MQ=255)
                                                                                                                                                                                                                            gggTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGa                                >  1:292035/1‑250 (MQ=255)
                                                                                                                                                                                                                            gggTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGa                                >  2:113486/1‑250 (MQ=255)
                                                                                                                                                                                                                             ggTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAAt                              <  1:146739/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTacac                                                                                 <  2:146900/189‑1 (MQ=255)
                                                                                                                                                                                                                                         tCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTacac                                                                                 >  1:146900/1‑189 (MQ=255)
                                                                                                                                                                                                                                                c‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATgg            >  2:276596/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      gATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGa        >  2:51547/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       aTTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGaaa      <  2:292035/251‑1 (MQ=255)
                                                                                                                                                                                                                                                            ccAGTTTTCTTATTCGAAAGCTCCATGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTGCC‑AGG‑AAAGTTTTTAACAATTCCTTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATTGCGGagatctt  <  2:528633/250‑7 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                            
AACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTG‑AAAAGGTTCACACCC‑AGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCG‑TTCGGCTAACTTCC‑AGG‑AAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGAAAGCTT  >  NC_000913/3947600‑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: ‑