Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,223,5970CT28.0% 83.9 / 15.9 50R11K (AGA→AAA) iclRDNA‑binding transcriptional repressor IclR
Reads supporting (aligned to +/- strand):  ref base C (5/31);  new base T (14/0);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 1.24e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGA  >  NC_000913/4223356‑4223844
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
gACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGtt                                                                                                                                                                                                                                                 <  2:112066/251‑1 (MQ=255)
  ccATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTc                                                                                                                                                                                                                                               <  1:229825/251‑1 (MQ=255)
      aaatgCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTTAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGATCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCgg                                                                                                                                                                                                                                           <  1:1985/247‑1 (MQ=255)
      aaaGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTGC‑GCGTTTCTCgg                                                                                                                                                                                                                                           >  1:8068/1‑251 (MQ=255)
               tgCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGCTTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTTACGCGTTTTGCggggatagggg                                                                                                                                                                                                                                  >  2:392789/1‑243 (MQ=255)
               tgCGCAGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGATGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCACTCCACTCCAGTAATTTCAGGCCACGCGTTACAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTGCcgccgtttgcgggggtggggg                                                                                                                                                                                                                                  >  2:166182/1‑233 (MQ=255)
                       ttGCCCAATATCCCAGTTCGCCACCCTGAAGCACGCAACCCTGCTGTTGCATCGTGGTTAGAAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGACATTGGATTCGGCAATCCACTCCAGTAATTTCGGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGAGCGGTGGCCACGGCGGG‑GTTTTTGTCGCGtttttacggtataggggggtacat                                                                                                                                                                                                                           >  2:439346/1‑229 (MQ=255)
                            cAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt                                                                                                                                                                                                                     <  1:71579/251‑1 (MQ=255)
                            cAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGGTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGGGGCAACGGCGGG‑TGTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt                                                                                                                                                                                                                     <  2:446806/251‑1 (MQ=255)
                              atggcccAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct                                                                                                                                                                                                                   <  1:119208/247‑1 (MQ=255)
                              aTGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGTTGGTGCGGTGGCAACGGCGGG‑GTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct                                                                                                                                                                                                                   <  1:203983/251‑1 (MQ=255)
                              aTGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct                                                                                                                                                                                                                   <  1:411087/251‑1 (MQ=255)
                                gTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct                                                                                                                                                                                                                 <  2:215236/251‑1 (MQ=255)
                                   ccAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTTCCCCGGTTTGCGGGGATGGGTggggccctggcaggctttttt                                                                                                                                                                                                              >  2:177510/1‑231 (MQ=255)
                                    cAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttttt                                                                                                                                                                                                             <  2:307470/251‑1 (MQ=255)
                                         cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGt                                                                                                                                                                                                        <  2:233463/251‑1 (MQ=255)
                                           ccAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTAt                                                                                                                                                                                                      <  2:139296/251‑1 (MQ=255)
                                             aaCCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCg                                                                                                                                                                                                    <  1:444811/251‑1 (MQ=255)
                                              aCCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGt                                                                                                                                                                                                   <  2:434856/251‑1 (MQ=255)
                                               cctgacgtccgaatgccggctgTTGCATCGTGGTGAGCAGGCGGGTGGG‑CGTGGGATTGGGTAACCCGGCTTGTTGCGCCAGTGCCGTGAGTGCCACACGGGCATTGGATTCGGCAATCCACGCCAGGAAGGTCAGGCCACGCGTTAAAGCCTGAACCTGTCCAGGCGCGGGTGCGGTGGCCACGGCGGG‑TTTTCTGCCGCGTTTCTCGGGAATGGGTGCGATCATGACAGTCTCTTTTTTCTGTATCGTg                                                                                                                                                                                                  <  2:272692/233‑1 (MQ=255)
                                                  gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGaa                                                                                                                                                                                               <  2:348167/251‑1 (MQ=255)
                                                        cGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTGCCGCGGTTTGCGGGGAGGGGGGGGGACATGACAGGCTCTTTTTTTTGTATCGGGGaaaaacat                                                                                                                                                                                         >  1:345224/1‑246 (MQ=255)
                                                         gAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TATTTTTCCCCGTTTTTTGGGAATGGGGGTGACCATTTCAGTCTCTTTTTTTTTTATCGTGGAAAtaattt                                                                                                                                                                                        >  1:342900/1‑251 (MQ=255)
                                                         gAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGACGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTGCCGCGTTTTGCTGGGATGGGTGCGACCATGGCAGTCTCTTTTTTTTGTATCGTGGAGAtcattt                                                                                                                                                                                        >  1:221429/1‑251 (MQ=255)
                                                            aCCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTTTGCCCCGTTTCGCGGGAATGGGTGCGGCCATGACAGTCTCTTTTTTTTTTTTCGTGGGAAtcattttct                                                                                                                                                                                     >  1:117662/1‑250 (MQ=255)
                                                              ccTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGGTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGTCGGG‑TTTTCTTCCGCGGTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc                                                                                                                                                                                                           <  1:274242/227‑1 (MQ=255)
                                                                 gctgTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTTCCGCGTTTCGCGGGGATGGGTGGGGCCATGACAGTCTCTTTTTTCTGTATCGGTGAAATCCTTTTTAttttt                                                                                                                                                                                >  1:357559/1‑251 (MQ=255)
                                                                            gTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTa                                                                                                                                                                     <  2:117662/251‑1 (MQ=255)
                                                                                       gctg‑TGGGCCGTGGAATTGGGTAGCCCGGCGTGTTTCGCCAGTCCCTTGAGAGCCACACTCCCATTGGAGTCGGCAATCCACTCCAGTAATTTCAGTCCACGCGTTAAATACGGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGTCGGG‑TTTTCTGCCGTGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGtt                                                                                                                                                          <  2:8068/248‑1 (MQ=255)
                                                                                                     gAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTTTGCCGCGTTTCGCGGGGATGGGGGGGACCCTGACAGGCTCTTTTTTCTGTTTTGGGGAAATCCTTTTTCTTTTTTTTGTTTTCTCATGCAATAATTAATTAACTGATcc                                                                                                                                             >  1:323258/1‑251 (MQ=255)
                                                                                                                               cAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTGACGCGTTTCGCGGGAAAGGGAGCTACCCTGTCAGTCTCTTTTTTCTTTATCGTGGAAATCATTTTCATTTTTATTGTTTTCTAATGCAATAGGTACTGAACTGATCCGCTGAGGT‑AATGTTGAACAAAtctcc                                                                                                                  >  1:54131/1‑250 (MQ=255)
                                                                                                                                               acacTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCg                                                                                                  <  1:177510/251‑1 (MQ=255)
                                                                                                                                                            aTTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGa                                                                                     <  1:134934/251‑1 (MQ=255)
                                                                                                                                                                   aaTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGGTTCGCGGGAATGGGGGCGGCCATGACAGTCTCTTTTTTTTGTATCGTGGAAATCATTTTTATTTTTATTTTTAGCTAATGCAATAGTTACTGAACTGGTCCGATGAGTT‑AATGTTGAACAAAACTCATGTTGCGTGGGGGGGGCTTTTACCCCAGATGCggtt                                                                              >  1:179715/1‑248 (MQ=255)
                                                                                                                                                                       cACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAt                                                                                                                                                                   <  2:329745/163‑1 (MQ=255)
                                                                                                                                                                                aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATgg                                                                  <  1:11348/250‑1 (MQ=255)
                                                                                                                                                                                          cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGCCGCGTTTCGCGGGAATGGGGGCGACCATGACAGTCTCTTTTTTTTGTATTGGGGGAATCATTTTCATTTTTTTTTTTTGCTTATGCAATAGGTACTTAACTGATCCGATGGGTT‑AATGTTTAAAAAATCTCATTTTTCGTGGGGGTCGCTTTTTCCCCAAATGCGTTTTTGCCCGTTTGGTTTGTTTAAtt                                                       >  1:2336/1‑251 (MQ=255)
                                                                                                                                                                                               gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTTTGCCGCGTTTCGCGGGGAAGGGGGCGACCATGACAGTCTCTTTTTTTTGTATCGTGGGAATAATTTTTATTTTTTTTGTTAGCTAAAGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAAAAAAACTCATGTTGCGTGGGGGTCGCTTTTTCCACAAATGCGGTTATGCCCGGATGGGTTGTTTGATTTttat                                                  >  1:290737/1‑251 (MQ=255)
                                                                                                                                                                                                 tAAAGACTGTACCTGTCCAGTCGCTGGTGCGGTGGCATCGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTTCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCTCTTTTACCACAGATGCGTTTATGCCATTATGGTTTGTTGAATTTttatt                                                 <  1:72640/250‑1 (MQ=255)
                                                                                                                                                                                                   aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTTTTCCGCGGCTTGCGGGGGTGGGGGCGGCCCTGAAAGGCTCTTTTTTTTTTATAGTGGGAATAATTTTTATTTTTTTTTTTAGCTTATGCCATAATTACTTAAATTATCCCATGAGTTTATTTTTTAAAAAATC‑CATGTTTGGGGGTGGTGGGTTTTTCCACAAATGCGTTTTTGCCCGGATTGTTTTTTGAAttttttttta                                               >  2:273634/1‑245 (MQ=255)
                                                                                                                                                                                                               gtccatgcGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAg                                   <  2:462225/243‑1 (MQ=255)
                                                                                                                                                                                                               gTCCATTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCGGTATCGTTGAAATGATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACGGATCCTAGGAGTT‑AAGGTTGAACAAATCGCATGTTGCGTGGTGGTTGCTTTTACCACAGATGCGTTTACGCCAGTATTGTTTGTTGAATTTTTATTAAATCTGGGTTGAGc                                  <  1:439346/251‑1 (MQ=255)
                                                                                                                                                                                                               gTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTTCTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGc                                  <  2:221429/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTGCCGCGTTTTCCGGGGAAGGGCGGGACCATGGCCGTCTCTTTTTTTTGTATAGCGGAAATCATTTTCTTTTTTATTTTTTTATAAATCAATAGTTACTGTAATGATCCGATGAGTT‑AATGTTGTAAAAATCTCCTGTTGCGTGGTGGGCGCTTTTTCCACAGATGCGTTTATGTCAGTATGGGTTGTTGAATTTTTTTTAACTTTGGGTTGAGCgg                                >  1:210058/1‑250 (MQ=255)
                                                                                                                                                                                                                       gCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCggg                          <  2:210058/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggTGCGGTGGCAACGGCGGG‑TTTTCTTCCGCGCTTCCCGGGGAAGGGGGCGGCCAAGACAGGCTCTTTTTTCTGTATCGTGGGGATCATTTTCATTTTTTTTGTTAGCTTTTGCAATATTTACTGAACTGATCCGCTTAGTT‑TATTTTGCACAAATCTCATGTTGCGTGGTGGTCGCTTTTTCCACCGCTGCGTTTCTGCCAGTATTGTTTTTTGAATTTTTTTTTAATCTTGGtttgggggggcgggggg                       >  2:101389/1‑236 (MQ=255)
                                                                                                                                                                                                                                            gg‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt     <  2:400536/251‑1 (MQ=255)
                                                                                                                                                                                                                                                tttCTGCCGCGTCTAGCGGGAATGGGTGCGACCATGACAGCCTCTTTTTTCTGTATAGTGGAAATCATTATCATTTTTATTGTTAGCTAATGCAATCGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTgg   >  2:149428/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGa  <  1:101389/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAGAGTTACTGAACTGAGCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACATATGCGTTTATGCCAGTATGGTTTGGTGAATTTTTATTAAATCTGGGGTGAGCGTGTCGGGAGCAAGTGTGATCATCAAAGTGGa  <  2:276207/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
GACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGG‑TGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGA  >  NC_000913/4223356‑4223844

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