Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,1420CA27.1% 100.9 / 14.1 48Q26K (CAA→AAA) fabRDNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base C (33/2);  new base A (0/13);  total (33/15)
Fisher's exact test for biased strand distribution p-value = 5.44e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCC  >  NC_000913/4160909‑4161388
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
gTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACccgc                                                                                                                                                                                                                                       >  1:61458/1‑251 (MQ=255)
                     tttttATTGGATTTTTATCGACGATCCTCTCCGTTTAAGCGGCAGTTTTCCGTTGTACATAAAAGACCCGGCAAAAAAATGGCAGTAAATATTTTTTATTGCTTTCCCTTTCATCCACAATAC‑TGGAGCAACCCAGTATGTTCATTCTCTGTTATAGTGCCAGCAGTATTTTTGA‑AAAGGATTCA‑GACATCGTG‑ATGGGCGAAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGc                                                                                                                                                                                                                   <  2:156130/249‑1 (MQ=255)
                     gTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGc                                                                                                                                                                                                                   >  1:446622/1‑250 (MQ=255)
                        ttATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAttt                                                                                                                                                                                                               >  1:33470/1‑251 (MQ=255)
                           ttGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc                                                                                                                                                                                                            >  1:157308/1‑251 (MQ=255)
                              gATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa                                                                                                                                                                                                          >  2:453528/1‑250 (MQ=255)
                               aTTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTACAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTAACGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAt                                                                                                                                                                                                        >  2:77133/1‑251 (MQ=255)
                               aTTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAt                                                                                                                                                                                                        >  1:465692/1‑251 (MQ=255)
                                      ccgacgaTACTTTCCGTTTAAGCGGCGGGTTCTCGCTGTACGTAAAAGAACCGGCCAAAGAAGTGCAGTAAATATGTTTTATGGTGTTCCCTTCCATTCACAATAC‑TGGACCAATCCAGTATGTTCATTCTCTGGTATAGGGCAAGCAGTACTGTTGG‑CAAGGATTCA‑GACATCGTG‑AGGGGCTTAA‑GAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTg                                                                                                                                                                                                  <  2:201561/249‑1 (MQ=255)
                                       cgacgaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCCTTTAGCCAATTAAGTgat                                                                                                                                                                                                >  2:325795/1‑249 (MQ=255)
                                        gacgaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                                                                                                                                                                                               >  1:146862/1‑251 (MQ=255)
                                           gaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc                                                                                                                                                                                            >  1:92691/1‑251 (MQ=255)
                                              aCTCTCCGTTGAAGCGGCAGGTTTCCGCTGTAAGTAAAAGAACCGGCAAAAGAATGGCAGTAAATTTGTTTTATTGCGTTACGGTCAATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCGGGTATAGTGCAAGCAGGATTTTGGG‑AAAGGATCCA‑GAAATCGTG‑AGGGGGGTAA‑GAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc                                                                                                                                                                                          <  1:85193/250‑1 (MQ=255)
                                                     gTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGc                                                                                                                                                                                  >  2:382645/1‑251 (MQ=255)
                                                                        gCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGa                                                                                                                                                                >  1:228095/1‑250 (MQ=255)
                                                                             aagaaaaaGACCCGGAAAAAGAATTGCATAAAATATTTTTTTTTGCGTTACGGTTCATT‑ACAAAACCTGGAGCAATCCAGTATGTTCATTCTCGGTTATAGGGCCAGCAGTGTTTTGGG‑AAAGGATCCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgc                                                                                                                                                          <  2:54056/247‑1 (MQ=255)
                                                                                          ggcaaaaaaatagaaGTAAATATTTTTTTTTGCGTTACCTTTCATTCACAATAC‑GGGAGCAATCCATTTTGTTCATTCGTTGGTATAGTCCCAGAAGT‑CTTTTGGGCAAGGATCAA‑AACATCGTG‑AGGGGCTAAA‑GAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGc                                                                                                                                             <  2:464561/237‑1 (MQ=255)
                                                                                             cAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAtt                                                                                                                                          >  2:468691/1‑251 (MQ=255)
                                                                                                ataatggCAGAAAATATGTTTTATTGTGTTACCGTTCATTCAAAATCC‑TGGAGAAACCCAGTATGTTCATTCTCTGGTATAGGGCCAGCAGTTCTTTGGG‑CA‑GGATTCAAGACATCGTG‑AGGGGCGTAA‑GGGCACAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCt                                                                                                                                       <  2:228095/245‑1 (MQ=255)
                                                                                                 gAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTc                                                                                                                                      >  2:306548/1‑251 (MQ=255)
                                                                                                  aaTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTcc                                                                                                                                     >  2:373367/1‑251 (MQ=255)
                                                                                                         gTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACctct                                                                                                                              >  2:504300/1‑251 (MQ=255)
                                                                                                           aaaTATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTg                                                                                                                                                                 >  1:208505/1‑214 (MQ=255)
                                                                                                             atatGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCttttt                                                                                                                          >  2:232307/1‑251 (MQ=255)
                                                                                                                         gctttACCTTTCATTCACAATAC‑GGGAGCAACCCTGTATGTTAATTCTCTGGTATAGTCCCACCAGTTCTTTGGG‑AAAGGATTCA‑GACATCGTG‑AGGGGCGAAA‑GGGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTATTATCGGCATTTc                                                                                                               <  2:446622/247‑1 (MQ=255)
                                                                                                                              ccctttcattcaCAATAC‑TGGAGCATTCCAGTTTGTTCATTCTCTGGTATGGTGCCAGCGTTACTTTGGG‑AAAGGATTCA‑GACATCGTG‑AGGGGCTTAA‑GAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGAc                                                                                                         <  2:278414/247‑1 (MQ=255)
                                                                                                                              acgtttcatacaAAATAC‑TGGAGAAACCAAGTATGTTCATTCCCTGGTATATTGCGAGCATTTCTTTTGG‑AAAGGATTCA‑AATATCGGG‑A‑GGGCTTAAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGAc                                                                                                         <  1:325795/247‑1 (MQ=255)
                                                                                                                                  ttcattcaAAATGC‑TGGAGCAACCAATTATGTTCATTCTTTGTTATAGTGCCAGCAGTACTTTTGG‑AAAGGATTCA‑GACATCGGG‑AGGGGCGTAA‑GAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTa                                                                                                      <  1:291997/250‑1 (MQ=255)
                                                                                                                                     catccacaatctgggataaagcctgTATGTTCATTCTCGGGTATAGTGCAAGAATTATTTTTGA‑CAAGGATTCA‑GACATGGTG‑ATGGGCGAAA‑GAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGan                                                                                                   <  2:252081/227‑2 (MQ=255)
                                                                                                                                     attcaCAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACg                                                                                                  >  1:406808/1‑251 (MQ=255)
                                                                                                                                                        aTCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATgg                                                                                >  1:286491/1‑251 (MQ=255)
                                                                                                                                                            aatattttCATTCTCTGTTATTGTGCCAGCATTTTTTTTGG‑AAAGGTTTCA‑GACATCGGG‑GTGGGCGGAA‑GAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtga                                                                            <  2:200833/245‑1 (MQ=255)
                                                                                                                                                               aTGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTc                                                                                                               >  1:161709/1‑213 (MQ=255)
                                                                                                                                                               aTGTTAATTCCCGGGTATGGTGCCGGCATTATTTTTGG‑AAAGGATTCA‑GCCATCGTG‑AGGGGGGTAA‑GAGCCAAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgatg                                                                          <  1:306548/250‑1 (MQ=255)
                                                                                                                                                                  ttCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgag                                                                        >  2:33207/1‑249 (MQ=255)
                                                                                                                                                                  ttCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgag                                                                        >  1:67397/1‑249 (MQ=255)
                                                                                                                                                                       tctcTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGtt                                                                 >  1:300037/1‑251 (MQ=255)
                                                                                                                                                                       tctcTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGCGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTTGACGAACTGGGTCTGACCCTGGTtgatgtgatggttt                                                                 >  1:461179/1‑242 (MQ=255)
                                                                                                                                                                                    gCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc                                                     >  2:117206/1‑250 (MQ=255)
                                                                                                                                                                                               ttttGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCAGGCGcc                                         >  2:334815/1‑251 (MQ=255)
                                                                                                                                                                                                tttGG‑AAAGAATTCA‑GACAT‑GTGGATGGGCGTAA‑GAGCAAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGcc                                         <  1:108353/250‑1 (MQ=255)
                                                                                                                                                                                                                 gACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTa                          >  2:169082/1‑250 (MQ=255)
                                                                                                                                                                                                                    aTCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGc                      >  1:186951/1‑251 (MQ=255)
                                                                                                                                                                                                                        gggagggtctaaa‑GAGCGAAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTTAACGCAGCTTCGCCAGCCTGAGTGTGCGTGAAGTGGCGCGTGAAGCGCGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCa                    <  1:158991/239‑1 (MQ=255)
                                                                                                                                                                                                                              ggCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCgg             >  1:166021/1‑251 (MQ=255)
                                                                                                                                                                                                                               gCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCggg            >  1:62897/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gagacaaacaaAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGtgtg      <  1:468691/245‑1 (MQ=255)
                                                                                                                                                                                                                                        gcgcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGGGAt    >  2:421055/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATcc  >  2:303344/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
GTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATAC‑TGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGG‑CAAGGATTCA‑GACATCGTG‑ATGGGCGTAA‑GAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCC  >  NC_000913/4160909‑4161388

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