Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,1480GA27.3% 76.6 / 15.2 44E28K (GAA→AAA) fabRDNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base G (27/5);  new base A (0/12);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 2.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGT  >  NC_000913/4160909‑4161381
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           
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)
                                                                                                                                                                                                                                              aaaGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGAGCGGGAAGCGGGCAtt                                                                                                                                   >  2:111605/1‑108 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           
GTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGG‑AGCAATCCAGTATGTTCATTCTCTG‑GTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGT  >  NC_000913/4160909‑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: ‑