Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,1480GA21.7% 85.1 / 10.7 46E28K (GAA→AAA) fabRDNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base G (35/1);  new base A (0/10);  total (35/11)
Fisher's exact test for biased strand distribution p-value = 2.70e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTC  >  NC_000913/4160923‑4161395
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        
cccTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCCGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCtggtgg                                                                                                                                                                                                                                 >  2:280695/1‑250 (MQ=255)
    gACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCCGCAGTACTTTTGG‑CAAGG‑ATTAAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACC‑GCCGTTCGCTGGTGGAAGcc                                                                                                                                                                                                                            >  1:263147/1‑251 (MQ=255)
     aCGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCg                                                                                                                                                                                                                           >  1:298561/1‑251 (MQ=255)
         tttATTGGCTTTTTATCGACGCTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGCATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGACAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCAtt                                                                                                                                                                                                                       >  2:184211/1‑251 (MQ=255)
                   ttttATCGACGATCCTCTCCGTTTAAGCGGCAGGTTTCCGCGGTACGTAAAAGAACCGGACAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCTTTCACAATAATGGAGAAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGAAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAGCAAACCA‑GCCGTACGCTGGTGGCAGCCGCATTTAGCCAAtt                                                                                                                                                                                                              >  2:267212/1‑250 (MQ=255)
                   ttttATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACATAAAAGAAACGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTAATTCTATGGT‑ATAGTGCCAGCAGTACATTTGG‑CCAGG‑ATTCAGTCATCGTGATGGGCGAAAGAACGCAAAAAAAAGAAAAAACCC‑GCCGTTCGCTGGCGGAAGCCGCATTTAGCCAAtt                                                                                                                                                                                                              >  2:263357/1‑250 (MQ=255)
                        acgacgaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGGAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGA‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTTCCGGGCGTAAGAGCGCAACAAAAAGCCAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGc                                                                                                                                                                                                        >  2:352982/2‑251 (MQ=255)
                                           aaGCGGCAGGTTTCCGCTGTACGTCAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCATACATCGTGCTGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGCAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGc                                                                                                                                                                                     >  2:427343/1‑251 (MQ=255)
                                                        ccgttgtagggaaaaGACCCGGCAAAAGAAGTGCAGAAAATATTTGTTATTGCGTTCCCGGTCATTCACAAGACGGGCGCAATACAGGTTGTTCAT‑GTCTGGTGTTAGTGCCA‑CTGTTCCTTTGG‑GAAGGCATCCGTAAAACGGGATGGGAGTAGGAGCGAAACAAAAAAAAAAAAGCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTTCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTg                                                                                                                                                                        <  1:163014/240‑1 (MQ=255)
                                                                     aaGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGAGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAACCAATTAAGTGCTAAACGCAGCTACGCCAGCCTGAGTTTGCGTGAGGTGGCGCGTg                                                                                                                                                            >  2:257392/1‑250 (MQ=255)
                                                                       gAACAGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTAACGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGAGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGcgcggga                                                                                                                                                           >  2:354429/1‑246 (MQ=255)
                                                                                 aaGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGc                                                                                                                                               >  1:386067/1‑251 (MQ=255)
                                                                                  aGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCt                                                                                                                                              >  1:308697/1‑251 (MQ=255)
                                                                                    aaTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg                                                                                                                                                                                                                                                          >  2:221102/1‑142 (MQ=255)
                                                                                    aaTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg                                                                                                                                                                                                                                                          <  1:221102/142‑1 (MQ=255)
                                                                                    aaTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTcc                                                                                                                                            >  1:88992/1‑251 (MQ=255)
                                                                                           gTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CACGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGACGCATTTAGCCAATTAAGTGATATACGCAGCTTCGCCAGCCTGAGTTTGCGGGAAGGGGCGCGTGACACGGGCCTTGCTCCCACctct                                                                                                                                     >  2:348893/1‑251 (MQ=255)
                                                                                           gTAAATATGTTCTATTGCGTTACCGTTAATTCACAATCCTGGCGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTCATGAGCGTAAGAGCGACACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGATTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGCGAAGAGGGCATTGCTCCCACctct                                                                                                                                     >  2:446091/1‑251 (MQ=255)
                                                                                             aaaTATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCg                                                                                                                                                                                                                                        >  1:341227/1‑150 (MQ=255)
                                                                                                   gTTTTATTGCGTTCCCTTTCATTCAAAATCTTGGAACAATCCGTTATGTTCATGCTCTGGT‑ATAGGGCCAGCAGTATTTTGGG‑AAAGA‑ATTCAGACTCCGTGATGGGCATAAGAGCGCAAAAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATc                                                                                                                              <  2:358158/250‑1 (MQ=255)
                                                                                                   gTTTTATTGCGTTACCTTTTTTTCAAAATATTGGACCCCTCCAGTATGTTCTTTCTCGGTT‑ATCGGCCCGGCAGT‑CTTTTGGGCAAGG‑ATTCAGACATCGTGCCGGGCGAAAGAGCGAAAAAAAAAAAAAAAACCC‑GCCGTTCGCTGGTGGAAGCTGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATc                                                                                                                              <  1:352982/250‑1 (MQ=255)
                                                                                                        aTTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCNGGCATTGCTCCCACCTCTTTTTATCGGCAtt                                                                                                                        >  2:239291/1‑251 (MQ=255)
                                                                                                           ggttaacgtttatttcaCAATACTGGACCAAACCAGTCTTTTCATTCTGTGGT‑ATAGTGCAAGCTGTGTTTTTGG‑AAAGG‑CTACAGCCGTCGTGGGGGGCGAAAGAGCGCACCAAAAAAAAAAACCCC‑GCCGTGCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAttacc                                                                                                                     <  2:128345/238‑4 (MQ=255)
                                                                                                            cGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTccc                                                                                                                    >  1:98263/1‑250 (MQ=255)
                                                                                                             ttttcgtttcattcaCAATACGGGACCATACAAGTATTGCCATTCTGTGTT‑ATATTGCCAGCATTATTTTGGG‑CAAGG‑ATCCAGAAATCGGGGGGGGGGAAAGAGC‑CAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc                                                                                                                     <  2:98263/242‑1 (MQ=255)
                                                                                                              ttACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTAAGCAAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCTTTTCTCCCACCTCTTTTTCTCGGCATTTcctcg                                                                                                                  >  2:125788/1‑248 (MQ=255)
                                                                                                               tACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcgcg                                                                                                                  >  1:244247/1‑250 (MQ=255)
                                                                                                                                ccgggAGCAACCCAGTATGTTCTTTCTCTGGT‑ATGGGGCCAGCATTACTTTGGG‑CAAGG‑ATTCAGACATCGGGAGGGGGGTAAGAGCGCAAAAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGt                                                                                                 <  1:101086/247‑1 (MQ=255)
                                                                                                                                  gggAGCAATCCATTAGGTTCATTCTCGGGT‑ATGGTGCCAGCAGTATTTTGGG‑CAAGG‑ATTCAAACATGGTGAGGGGCGAAAGAGCCAAAAAAAAAGAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTg                                                                                              <  2:344396/250‑1 (MQ=255)
                                                                                                                                     agaaaTCCAGTTTGTTAATTCTCTGTT‑ATAGTGCCAGCATTACTTTTGG‑AAGGG‑ATTCAGACATCGTGAGGGGGGAAAGAGCGCAAAAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGAcc                                                                                           <  1:239291/248‑1 (MQ=255)
                                                                                                                                                aTGTTCATACTCTGGT‑ATAGTGCCAGCAGTACTTTTAG‑CCAGG‑ATTCAGACATCGTGATCGGCGTAAGAGCGCAACACAAAGACAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTACTGAACGCAGCTAAGCCAGCCTGAGTTCGCGTGCAGTGGCGCGTGCAACGAGCATAGTTCCCACCTCTTTTTTTCGGCATTACCGCGAGGACGACGAACTgggcattcacttcgttgctg                                                                                 >  2:56461/1‑233 (MQ=255)
                                                                                                                                                     cATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGACGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTCGCGGCATTTCCGCGACGTAGACGAACTTGGGCGGACCATGGTTGAGGGGAGCg                                                                           >  2:363516/1‑251 (MQ=255)
                                                                                                                                                              gT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGCGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGAGACGTAGACGAACTGGGTCTGACCATGGCAGTTGAGAGCGGTTTAATGc                                                                  >  2:445747/1‑251 (MQ=255)
                                                                                                                                                                         agcagTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGAGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGACGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATAGGCATTTCCGCGACGTAGCCGAAATGGGTCTGCCCATGGTTGATGAGAGCGGTTTACAGATACGCCAAc                                                         >  2:257008/1‑250 (MQ=255)
                                                                                                                                                                             gTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtg                                                                                    >  1:361093/1‑219 (MQ=255)
                                                                                                                                                                                 ttttGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCATCAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGcc                                                >  1:65385/1‑251 (MQ=255)
                                                                                                                                                                                     gg‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGgc                                            >  1:136505/1‑251 (MQ=255)
                                                                                                                                                                                                   gACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGGGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTAt                                >  2:342136/1‑251 (MQ=255)
                                                                                                                                                                                                       cctggATGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTATATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGcc                            <  1:292844/247‑1 (MQ=255)
                                                                                                                                                                                                        cGTGATCGGCGTAAGAGCGCAACCAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTACGTGCTGAACGCAGCTTCGCCAGCCTGTGTTTGCGTGAAGTGGAGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGCCATTTCCGCGACGTAGACGAACTGGGACTGACCATGGTTGCTGAGGGCGGGTTAATGCTAAGCCAACTCATGCACCAGGCGCGTCAGCGTATCGcc                            >  2:480000/1‑250 (MQ=255)
                                                                                                                                                                                                           gtgggggcgaagggggacacaaAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCaaa                         <  1:267212/233‑1 (MQ=255)
                                                                                                                                                                                                             tGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGc                      >  1:238065/1‑251 (MQ=255)
                                                                                                                                                                                                                gagTAAGCGCGCAACAAAAAAAAAAATCCC‑TCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCTGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCggg                   <  1:20264/249‑1 (MQ=255)
                                                                                                                                                                                                                   taagagcaaaaaaaaaaaaaaaaaCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGa                  <  1:354771/237‑1 (MQ=255)
                                                                                                                                                                                                                   tAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCgggga                 >  2:354771/1‑248 (MQ=255)
                                                                                                                                                                                                                    aaGAGCGCAACAAACAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGCGAGCGGTTTAATGCTACGCCAACTCATGCGCCCGGCGCGTCAGCGTATCGCCAAAGGCGGGAgtg               >  2:58734/1‑251 (MQ=255)
                                                                                                                                                                                                                            aacaaAAAGAAAAAACCC‑GCCGTTCGCTGGAGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGc       >  1:387575/1‑251 (MQ=255)
                                                                                                                                                                                                                             acaaAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCCCCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGGCAGCGCATCGCCAAAGGCGGGAgagggagccgcc      >  2:116493/1‑240 (MQ=255)
                                                                                                                                                                                                                                 aaaGAAAAAACCC‑GCCGTTCGCTGGAGGAAGCCGCATTTAACCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTc  >  1:145703/1‑251 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        
CCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGG‑CAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTC  >  NC_000913/4160923‑4161395

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