Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,1420CA31.8% 71.0 / 13.0 44Q26K (CAA→AAA) fabRDNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base C (27/3);  new base A (0/14);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 5.92e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGT  >  NC_000913/4160896‑4161381
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
ctttgctttgTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaa                                                                                                                                                                                                                                                >  1:457427/1‑250 (MQ=255)
             gTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACcc‑‑gc                                                                                                                                                                                                                                >  1:439593/1‑251 (MQ=255)
               ttcttcCTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACcc‑‑gcc                                                                                                                                                                                                                               <  2:254602/250‑1 (MQ=255)
                          ccccTGACGTTTTATTGGTTTTTTATCGACGATTCTCTCCTTTTAAGCGGCAGGTTTCCGCTGTACTAAAAAGACCCGCCCAAAGAATTGCAGAAAATATTTTTTATTGCGTTACCTTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTTCTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCC‑‑GCCGTTCGCtggtgg                                                                                                                                                                                                                   <  1:295137/250‑1 (MQ=255)
                             cTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGa                                                                                                                                                                                                                  >  1:220388/1‑249 (MQ=255)
                                                       cgaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCAAATTAAGTGCTGa                                                                                                                                                                                       >  2:315257/1‑250 (MQ=255)
                                                             tctcCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc                                                                                                                                                                                >  2:446342/1‑251 (MQ=255)
                                                                         cGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTg                                                                                                                                                                    >  1:95154/1‑251 (MQ=255)
                                                                                 tcccGCTTTACGAAAAAGAGCCGGCAAAAAAATTGAAGTAAATTTTTTTTATTGCTTTCCTGTTCATTCAAAAT‑CTGGGAGAAATCCAGTATGTTCATTCCTTG‑GTATAGTGCCAGCAGTTCTTTGGGCAAGGATTCAGACATCTTGAGGGGGGTAAGAGCGCAAAAAAAAAAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGc                                                                                                                                                             <  1:77545/248‑1 (MQ=255)
                                                                                    cGCTGTACGTAAAAGAACCGGCCAAAGAATTGAAGTAAATATTTTTTATTGCTTTCCCGTTCATTCACAATACGGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGGGCAAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGa                                                                                                                                                         <  1:40530/251‑1 (MQ=255)
                                                                                           cGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgcg                                                                                                                                                  >  2:209026/1‑251 (MQ=255)
                                                                                              aaaaGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGt                                                                                                                                               >  1:413108/1‑251 (MQ=255)
                                                                                                aaGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTg                                                                                                                                              >  2:458629/1‑250 (MQ=255)
                                                                                                  gAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGa                                                                                                                                             >  1:118388/1‑249 (MQ=255)
                                                                                                          aaaaGAATTGCATTAAATATGTTTTATTGCGTTCCCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGGGCCGGCAGTATTTTTGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGAAACAAAAAAAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAtt                                                                                                                                   <  1:383866/250‑1 (MQ=255)
                                                                                                             aGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCt                                                                                                                                >  1:384628/1‑251 (MQ=255)
                                                                                                              gAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTc                                                                                                                               >  1:117478/1‑251 (MQ=255)
                                                                                                                aTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTccc                                                                                                                             >  1:211281/1‑251 (MQ=255)
                                                                                                                     attaaattttttttttTGCGTTACTTTTCATTAACAATACTGG‑AGAAATCCGGTATGTACATTCT‑TGTGATTATTACCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGTTGGGCGTAAGACCACAAAAAAAAAAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACctc                                                                                                                        <  1:267283/241‑1 (MQ=255)
                                                                                                                      gTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTACAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGCGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCCTTGCTCCCACctct                                                                                                                       >  2:504504/1‑251 (MQ=255)
                                                                                                                         gatatGTTTTATAGCGTTACCGTTCATTCCCAATACTGG‑AGCACTCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGCCAAGGATTCAGGCATCGTGATAGGCGTAAGAGAGCAACAAAATGAAAAAACCC‑‑GCCGTTCGAAGTTGGAAGCCGCATTTAGCCAATAAAGTGCTGAACGCAGCTTCGCCAGCATGAGTTTACGTGCAGTGGCGTGTGAAGCGGGAATTggtacgacatctttt                                                                                                                    >  2:359594/2‑237 (MQ=255)
                                                                                                                         aaTATGTTTTCTTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCtttt                                                                                                                    >  2:178978/1‑251 (MQ=255)
                                                                                                                             tGTTTTTTTGCGTTCCCGTTCATTCACAATACGGG‑AGCAACCCAGTATTTTCATTCTCTG‑GTATCGCGCCAGCAGTCCTTTGGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCACAAAAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCTAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTAt                                                                                                                 <  1:488488/250‑1 (MQ=255)
                                                                                                                                   aTTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAtt                                                                                                          >  1:116572/1‑251 (MQ=255)
                                                                                                                                           cccGTTCATTCAAAATACTGG‑AGAAATCCAGTATGTTCATTCTCGG‑GTATAGTGCCGGGAATCCTTTTGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGAc                                                                                                  <  1:209026/250‑1 (MQ=255)
                                                                                                                                                        aaTACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGACGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGAATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTGGACGAACTgg                                                                                     >  2:10964/1‑251 (MQ=255)
                                                                                                                                                                    aaTCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATg                                                                          >  1:319585/1‑251 (MQ=255)
                                                                                                                                                                      tCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGt                                                                        >  1:293815/1‑251 (MQ=255)
                                                                                                                                                                        cggaaTGTTCATTCTCTG‑ATATAGTGCAAGCAGTATTTTGGGAAAGGTTTCAGACAGCGGGAGGGGCGAAAGAGCGCAAAAAAAAGAAAAACCCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtg                                                                      <  1:315257/247‑1 (MQ=255)
                                                                                                                                                                           tATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgatg                                                                   >  1:63216/1‑251 (MQ=255)
                                                                                                                                                                                         g‑gTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATTGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATg                                                     >  2:371224/1‑251 (MQ=255)
                                                                                                                                                                                             ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCt                                                   >  1:229836/1‑250 (MQ=255)
                                                                                                                                                                                               aGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACg                                                >  1:144165/1‑251 (MQ=255)
                                                                                                                                                                                                    cagcagTACTTTGGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAACCCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAAc                                           <  1:10964/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaaggttacAGACGTCGTGATGGGGGAAAGAAC‑CAAAAAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGgcgc                            <  2:118388/243‑1 (MQ=255)
                                                                                                                                                                                                                   aaaGGATTCAGACATCGTGATGGGCGAAAGGGCGAAACAAAAAAAAAAACCCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGgcgc                            <  1:483485/250‑1 (MQ=255)
                                                                                                                                                                                                                       gATTCAGACATCGTGATGGGCGTAAGAGCGAAAAAAAAAAAAAAACCCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCa                        <  1:335110/251‑1 (MQ=255)
                                                                                                                                                                                                                               cagggggaggggCGAAACAGCGAAAAAAAAAGAAAAACCCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTagcg                <  2:293815/242‑4 (MQ=255)
                                                                                                                                                                                                                                 gcgtgaggggaggaaGAGCGCAAAAAAAAAAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGcc              <  2:457427/238‑1 (MQ=255)
                                                                                                                                                                                                                                      gtggggaaaggagCGCAAAAAAAAAAAAAACCCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAgg         <  1:178978/241‑1 (MQ=255)
                                                                                                                                                                                                                                      gtatgggatagtgagcAAAAAAAAAAAAAAACCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGATACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCaaa           <  2:319538/237‑1 (MQ=255)
                                                                                                                                                                                                                                        gggCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCg       >  1:141216/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGTCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCCAAGGCggg     >  1:355500/1‑249 (MQ=255)
                                                                                                                                                                                                                                             tAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAgt  >  1:109139/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
CTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGT  >  NC_000913/4160896‑4161381

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