Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,223,6440CT21.4% 170.7 / 17.1 70intergenic (‑16/‑184)iclR/metHDNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base C (13/42);  new base T (14/1);  total (27/43)
Fisher's exact test for biased strand distribution p-value = 1.22e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGT  >  NC_000913/4223406‑4223893
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                            
ggcGCACGAAACCTGGCTGTTTCAGCGTGGTTAGCAGGCTGTGGGTCGTGGAATTGGGTAACCCCGCTTGTGGCGCCAGTACCGTGAGTGCCACTCTGCCATTGGATTCGGCAATCCGCTCCAGTAATTTCAGGCCACGCTGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACTGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCGCTTCTTTCTGTATCGTGGaa                                                                                                                                                                                                                                                 <  2:205936/249‑1 (MQ=255)
    cACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtca                                                                                                                                                                                                                                             <  2:468784/251‑1 (MQ=255)
      cGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcatt                                                                                                                                                                                                                                          <  1:551587/251‑1 (MQ=255)
      cGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcatt                                                                                                                                                                                                                                          <  2:7659/251‑1 (MQ=255)
            ccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcattttcatt                                                                                                                                                                                                                                    <  1:206509/251‑1 (MQ=255)
            ccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcattttcatt                                                                                                                                                                                                                                    <  2:325416/251‑1 (MQ=255)
               gctgTGGCATCGGGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCTGCTTGTTGCGCCAGTTCAGTGAGGTCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACTTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCTTTTCGCGGGAATAGGTGCTAGCATGACAGTCTCTTGTGTTTGTATCGTGGAAATCATTTTCAttttt                                                                                                                                                                                                                                 <  2:154219/251‑1 (MQ=255)
                  gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGTCGCGTTTCGCGGGGAAGGGGGGGACCATGACAGGCTCTATTTTTTGTAtggggggaagcatttttaattttttt                                                                                                                                                                                                                                >  1:112225/1‑226 (MQ=255)
                            ggTTAGCAGGCGGTGGGCCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAt                                                                                                                                                                                                                   <  1:485124/251‑1 (MQ=255)
                                   aGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAg                                                                                                                                                                                                            <  1:428231/251‑1 (MQ=255)
                                   aGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTTCCGCGCTTCGCCGGGGAGGGGGCGGCCCCGACAAGCTCTTTTTTTTTTATCGTGGGAATCTTTTTTATTTTTTTTTTTAGCTAATGCCATa                                                                                                                                                                                                             >  2:318618/1‑250 (MQ=255)
                                        gTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCTGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATTACAGTCTCGTTTTTCTGTATCGTGGAAATCATTTTCATTTTGATTGTTAGCTAATGCAATAGTTact                                                                                                                                                                                                       <  2:291752/251‑1 (MQ=255)
                                         ggggTCGTGGAATTGGTTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTTCCATTGGATTCGGCAATCCACGCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTACTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTtctg                                                                                                                                                                                                      <  1:258351/250‑1 (MQ=255)
                                           ggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGGTTTGCGGGGATGGGGGGGACCATGACAGGCTCTTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTTTTTGCTAATGAAATAgattctga                                                                                                                                                                                                     >  1:82515/1‑243 (MQ=255)
                                                      tGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGAt                                                                                                                                                                                          <  1:38935/250‑1 (MQ=255)
                                                       gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATg                                                                                                                                                                                         <  2:110451/250‑1 (MQ=255)
                                                         gTAACCCGGCTTGTTGCGCCAGTGCCGTGATTGCCACTCTTCCATTGGATTCGGCAATCCACGCGAGTAATTTCAGGCCACGCGTTAAAGTCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTTCGACCATGTCAGTCTCTTTTTTCTGTATCGTGCAAATCATTTTCATTTTTATTTTTAGCTAAGGTAATAGTTACTGATCTGAGCCGATGAGt                                                                                                                                                                                      <  2:448816/251‑1 (MQ=255)
                                                            ccccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGTCCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTaa                                                                                                                                                                                   <  2:500581/250‑1 (MQ=255)
                                                             cccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGCGGGGATGGGGGCGGACCGGGCAGGCTCTTTTTTTTGTATTGTTGGAAACCTTTTTATTTTTTTTTTTTGCTAATGCTAAAATTACTTAACTGAttctgtgggttta                                                                                                                                                                                   >  1:427609/1‑238 (MQ=255)
                                                                 gCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGGATGGGGGCGGCCATGACAGTCTCTTTTTTTTGTATCGGGGGAATAATTTTCATTTTTTTTTTTATCTAATGGCATAGTTACTGAAATTATCCGATgtggtaatattg                                                                                                                                                                              >  1:283451/1‑240 (MQ=255)
                                                                    tgttgCGCCAGTTCCGCGAGTGACACACTGCCATTGCATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCAGGATTTTTGGCGGGTTTCGACGGGATGGGTGCTGCCCCAGCAGTGCCTTTTTTTTGTTTCGTGGAAATCAATTTTCCTTTTTTTGTTATCTAATGCCATAATTTCTGAAATTATCCgcttatttacaatttaac                                                                                                                                                                           >  2:224854/1‑234 (MQ=255)
                                                                        gcgcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTTAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAAt                                                                                                                                                                       <  1:541705/251‑1 (MQ=255)
                                                                             atttCCGTGAGTACCACACTGCCATTGGATTCGGCAATCCACTCCAGGAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGTCAGTCTCTTTGTTCTGTATCGTGGAAATCTTTTTCCTTTTTATTGTTAGCTAAAGCAATAGTTACTGAACTGATCCGATGAGTGGATGTTGAACACAtctcct                                                                                                                                                                  <  2:502686/249‑3 (MQ=255)
                                                                             aGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAt                                                                                                                                                                  <  1:229811/251‑1 (MQ=255)
                                                                               ttCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGGATGGGTGCGACCATGACAGTCTCTTTTTTTTGTATCGTGGGAATCATTTTCATTTTTTTTTTTTGCTAATGCAATAATTACTGAAATGATCCGATGAGTTAATTTTTAAAAAATATCATGt                                                                                                                                                                >  1:195109/1‑251 (MQ=255)
                                                                                            acacTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCTGCGTTTCGCGGGGATGGGGGCGGCCATGCCAGTCTCTTTTTTTTGTATTGTGGGAATCCATTTCTTTTTTTTTTTTTTTTAATGCAATAATTATTTTACTGCTCCCCTGAGTTTATGTTGTAAAAAACTCATGTTTCgtggggggcg                                                                                                                                                   >  2:355006/1‑245 (MQ=255)
                                                                                              acTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCCTGACAGTCTCTTTTTTTTGTTTCGTGGGAATCATTTTTATTTTTTTTTTTAGCTAATGCAATAGTTACTGAACTGATCCCAGGAGTTAATGTTGAACAAATCTCATGTTGCGGGGTGGTCGCt                                                                                                                                                 >  1:380082/1‑251 (MQ=255)
                                                                                                 gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGCGGGGATGGGGGGGGACCTGGCAGTTTCTTTTTTTTGTTTCGTGGAAATAATTTTTATTTTTTTTGTTAGCTTATGCAATAGTTACTGAACTGATCCGATGAGATTATGTTGTACAAATATCATGTTTCgtggggggtgctttt                                                                                                                                              >  2:457351/1‑240 (MQ=255)
                                                                                                    aTTGGATTCGGCATTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGGCTGAACCTGTCCAGTCGCTGTTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGATCAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc                                                                                                                                           <  1:574996/251‑1 (MQ=255)
                                                                                                            cGGCAATCTGCTCCAGTAATTGCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCAGTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGc                                                                                                                                   <  2:383682/251‑1 (MQ=255)
                                                                                                               cAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGtt                                                                                                                                <  1:435214/251‑1 (MQ=255)
                                                                                                               cAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGtt                                                                                                                                <  1:93682/251‑1 (MQ=255)
                                                                                                                aaTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGttt                                                                                                                               <  2:195447/251‑1 (MQ=255)
                                                                                                                   ccACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTGTTCTGTATCGTGGAAATCATTTTCAGTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTAt                                                                                                                             <  1:436295/250‑1 (MQ=255)
                                                                                                                          aGTAATGTCAGGCCACGCGTTCAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTTCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTAt                                                                                                                     <  2:195367/251‑1 (MQ=255)
                                                                                                                             aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCCCGCTTCGCGGGGATGGGGGCGACCCTGACAGTCTCTTTTTTTTTTATCGTGGAAATCATTTTTATTTTTTTTTTTTGCTAATGCAATAGTTACTGAACTGATCCGCTGAGTTTATGTTGTACAAAAATCATGTTGCGTGGTGGGCGCTTTTTCCACAGATGCGTTTTTGCCAGTATGGt                                                                                                                  >  1:559963/1‑251 (MQ=255)
                                                                                                                              aTTTCAGGCCAAGCGTTAAAGACTGAACCGGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTGTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGATTATGCCAGTATGGtt                                                                                                                 <  2:37108/251‑1 (MQ=255)
                                                                                                                                  cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgtt                                                                                                             <  1:457351/251‑1 (MQ=255)
                                                                                                                                              tAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttatt                                                                                                  <  1:241706/250‑1 (MQ=255)
                                                                                                                                                 aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTCTTCCCGGGATTGGGGGGGCCATGAAAGTCTCTTTTTTTTTTTTCGTGGAAATCATTTTCTTTTTTTTTGTTAGCTTAATCAATAGTTAATGAACACATTCCATGGGGTTATGTTTAACAAAACTTATGTTTCGCGGGGGTCGCTTTTTACACAGATGTGTTTATGACCGTGTGGGTTGTTTTATTTTTTTTaaa                                                                                               >  2:135221/1‑250 (MQ=255)
                                                                                                                                                      gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTTCCGCGTTTCGCGGGGATGGGGGGGGCCATGACAGTCTCTTTTTTCTGTATCGTGGGAATCATTTTTATTTTTTTTTTTAGCTTATGCAATAGTTACTGAACTGATCCGATGAGTTTATGTTGGACAAATCTCATTTTGCGTGGGGGTCGCTTTTTCCCCAGATGCGTTTTTGCCCGGTTGGTTTGTTTGATTTTTTTTTAATATggg                                                                                         >  1:454460/1‑251 (MQ=255)
                                                                                                                                                        aCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGtt                                                                                       <  2:554320/251‑1 (MQ=255)
                                                                                                                                                          cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCATGACAGTCTCTTTTTTTTGTATTGTGGGAAACATTTTTATTTTTATTGTTAGCTAATGCAATAGTTACTGAAAATAACCGATGGGTTAATGGTGAAAAAAACTCATTTTGCGTGGGGGTCGCTTTTATCACAGAAGCGTTTATGCCCGTATGGGTTGTTGAATTTTTATTAAATTTTGGTTgg                                                                                     >  1:99498/1‑250 (MQ=255)
                                                                                                                                                           tGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTTGCGGGGATGGGGGCGGCCCAGACAGTCTCTTTTTTCTGTTTCGTGGGAATCATTTTTCTTTTTTTTTTTTGCTAATGCAATAGTTACTGAACTGATCCGATGTGTTAATTTTTTAAAAATCTCATGt                                                                                                                                                                >  2:235537/1‑175 (MQ=255)
                                                                                                                                                           tGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGt                                                                                                                                                                <  1:235537/175‑1 (MQ=255)
                                                                                                                                                                aGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGGTTCGCGGGAATGGGGGCGGCCCTGACAGTCTCTTTTTTTTGTATCGTGGAAATCATTTTTATTTTTTTTTTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTTATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTTCCACAGATGCGTTTATGCCCGTATGGTTTGTTGAATTTTTTTTTAATCTGGGTTGAGCgggt                                                                               >  1:32594/1‑248 (MQ=255)
                                                                                                                                                                 gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAt                                                                                                                                                                                                                                               <  2:177324/92‑1 (MQ=255)
                                                                                                                                                                       ggTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGa                                                                          <  1:304369/249‑1 (MQ=255)
                                                                                                                                                                          gCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCaa                                                                      <  1:331460/250‑1 (MQ=255)
                                                                                                                                                                          gCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAg                                                                     <  2:356975/251‑1 (MQ=255)
                                                                                                                                                                                       cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCaaa                                                        <  2:454460/251‑1 (MQ=255)
                                                                                                                                                                                          gTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTg                                                     >  2:1864/1‑251 (MQ=255)
                                                                                                                                                                                                              cgGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTgcgc                                     >  1:343663/1‑249 (MQ=255)
                                                                                                                                                                                                                 ggAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGtt                                 <  1:355006/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGtt                                 <  1:407480/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaTGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCATTTaa                               <  2:492027/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATTCCAGTATGGTGTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaaa                              <  1:462656/251‑1 (MQ=255)
                                                                                                                                                                                                                      gggTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATg                            >  1:77278/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGCACAACTGCGTGCGCAGTTAAATg                            >  2:559700/1‑250 (MQ=255)
                                                                                                                                                                                                                        gTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGa                           >  2:13367/1‑250 (MQ=255)
                                                                                                                                                                                                                             gACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGt                       >  1:558268/1‑250 (MQ=255)
                                                                                                                                                                                                                             gACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGt                       >  1:294680/1‑250 (MQ=255)
                                                                                                                                                                                                                               aCCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAACTCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGACAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCCAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAt                     >  1:484258/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg                >  1:140916/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aCAGTCTCTTTTTTCTGTATCGTGGAAAGCATGTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt               <  2:85313/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACGGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATTCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAGCTGCGTGCGCAGTTAAATGAACGTATTCTGGTg              <  2:159267/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ctgttttCTGTATCGTGGAAAGCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCTCAGATGCGTGTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTgt          <  2:564653/246‑2 (MQ=255)
                                                                                                                                                                                                                                            ctTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTg           <  1:290005/248‑1 (MQ=255)
                                                                                                                                                                                                                                                  tCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGt  >  2:538223/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGt  >  1:354610/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                            
GACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGT  >  NC_000913/4223406‑4223893

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