Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,1480GA28.6% 71.9 / 12.7 42E28K (GAA→AAA) fabRDNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base G (28/2);  new base A (1/11);  total (29/13)
Fisher's exact test for biased strand distribution p-value = 2.06e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCT  >  NC_000913/4160916‑4161394
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
aCTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGTCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTACATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCAAGTATGTTCA‑TTCTCTGGT‑TTAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAACAAGAAAAAACCCGCCGTTcc                                                                                                                                                                                                                                       >  2:172253/1‑249 (MQ=255)
               ttttttTGGTTTTTTACCGACGATACTCTCCGTTAAAGCGGCAGATTTCCTCTGTACGTAAAAGACCCGGCCAAAGATTTGAAGAAATTATGTTTTATGTCTTTCCCTTC‑CATTCAAAATACTGGAGAAACCCAGTATTTTCA‑TTCCCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGAAAGCGTGAGGGGCTTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCa                                                                                                                                                                                                                        <  1:299648/250‑1 (MQ=255)
               ttttATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACACGCCGTTCGATGGTGGAAGCCGCAt                                                                                                                                                                                                                       >  1:349349/1‑251 (MQ=255)
               ttttATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGAGATGGGCGTAAGAGCGAAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt                                                                                                                                                                                                                       >  1:236467/1‑251 (MQ=255)
                              accgacgaatctctCCATTTAAGGGGCAGGTTTCTGCTGAACGAAAAAGAACCGACAAAAGATTTGCAGTAAATATTTTTTATTGTTTTACCTTT‑CATTCAAAATACGGGAGCAATCTAGTAT‑TTCAATCCTCTGGT‑ATAGTCCCAGCAGTACTTTTGAAAAGGA‑TTCAGACATGGTGATGGGCGAAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt                                                                                                                                                                                                         <  2:15997/240‑1 (MQ=255)
                              aTCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACACAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTg                                                                                                                                                                                                        >  1:152225/1‑251 (MQ=255)
                                 gacgaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CCTTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCCAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGATGGGGGCAGCCGCATTTAGCCAATTAAATGCTg                                                                                                                                                                                                     >  2:272513/1‑251 (MQ=255)
                                              gTTTAAGCGGCAGGATTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGAAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGc                                                                                                                                                                                        >  1:187548/1‑251 (MQ=255)
                                                  aaGCGGCAGGTTTCCGCTGTACGTAAAAGACCCGGCAAAAGAATTGCAGTAAATATTTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAACCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTTCTTTTGGCAAGGA‑TTCAGACATCGGGATGGGCGTAAGAGCGAAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGc                                                                                                                                                                                    <  1:5282/251‑1 (MQ=255)
                                                  aaGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGc                                                                                                                                                                                    >  1:2844/1‑251 (MQ=255)
                                                    ggggggggTTTCCTGTGTACGAAAAAGAACCGGCCAAAGATTTGCAGAAAATATGTTTTATTGGGTTGCCGTT‑CATTCAAAATACGGGAGCAACCCATTAGGTTCA‑TCCTCTGGA‑ATAGTGCAACAATTACTTT‑GGGAAGGAATTCAGACATCGGGGTGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCt                                                                                                                                                                                  <  2:125223/245‑1 (MQ=255)
                                                                 gCTGTACGTAAAAGAACCGGCCAACGAATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGAGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCATTTAAGTGCTGCACGCAGCTTCGCCAGCCTGGGTTTGCGTgc                                                                                                                                                                      >  2:265892/1‑249 (MQ=255)
                                                                    gTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGA‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGCCATCGTGATGAGCGAAAGCACGCAACAAAAAAACAAAACCCGCCGTTCGCTAGTGGAAGCCGCATTAAGCCAATTAAGTGCTGAACGCAGCTGCGCCAGCCTGAGTTTGCATGAAggg                                                                                                                                                                  >  2:249338/1‑249 (MQ=255)
                                                                         aaaaaaacccGGCAAAAAAATTGAAGTAAATATTTTTTATGGCGTTACCGTT‑CATTCAAAAAACTGGAGCAATCAAGTATTTTCA‑TTCTTTGGT‑ATAGTGCCAGCAGTCCTTTGGACAAGGA‑TTAAGACACCGTGAGGGGCGAAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAGTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGcgcg                                                                                                                                                             <  1:397937/243‑1 (MQ=255)
                                                                                           aaTTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGTAGTGGCGCGTGAAGCGGTCATTGCTcc                                                                                                                                           >  2:315062/1‑251 (MQ=255)
                                                                                            aTTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTccc                                                                                                                                          >  1:355443/1‑251 (MQ=255)
                                                                                                    aaaTATGTTTTATTGC‑TTTCCGTTTCATCCACAATACGGGAGCAATCCAGTTGGTTCA‑TTCTCTGGT‑ATAGTGCCAGCATTACTTTTGGCAGGGA‑TACAGATATCGGGAGGGGGTTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCttt                                                                                                                                  <  1:272513/251‑1 (MQ=255)
                                                                                                       tatGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCACTCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAACAAGAAAAAACCCGCCGTTCGCTGGTGGAAGcc                                                                                                                                                                                                                           >  2:211626/1‑159 (MQ=255)
                                                                                                              tATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAt                                                                                                                        >  1:182817/1‑251 (MQ=255)
                                                                                                              tATTGCGTTACAGTT‑AATTGACAATACTGGTGCAATCCAGTATGTTCA‑GTCTCTGGT‑ATAGTACCAGCAGTACTTTTGGCAATGA‑TTCAGGCATCGTGATGGGCGTAAGAGCGCAACAAAAAGAACAAACCCGCCGTTCGCTGGTGGCAGCCGCATTTAGCCAATTCAGTGCTGATCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCACACCTCTTTTTATCGGCAt                                                                                                                        >  2:305665/1‑251 (MQ=255)
                                                                                                                ttGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTc                                                                                                                                                                                                                                        >  1:60267/1‑137 (MQ=255)
                                                                                                                    gTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGCGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcgc                                                                                                                  >  1:279208/1‑251 (MQ=255)
                                                                                                                       cccGTT‑CATTCAAATTCCTGGAGCAATCAAGTTTGTTCA‑TTCTCGGGT‑ATGGTGCCAGCGTTATTTTGGGAAAGAA‑TTAAGACATCGGGATGGGCGTAAGGGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGAc                                                                                                               <  1:431443/250‑1 (MQ=255)
                                                                                                                         cGTT‑CATTCCCACTACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGAACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTAATCCACCCTCTTTTTATAGGCATTTCCGCGAcat                                                                                                             >  2:378729/1‑249 (MQ=255)
                                                                                                                             tcttccacaatcatGGAACAATCCATTATGTTCA‑TCCT‑TGGTTATAGTGCCAGCAGTTCTTTTGGCAAGGA‑TCCAGCCATCGTGATGGGCTTAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAg                                                                                                           <  1:265892/237‑1 (MQ=255)
                                                                                                                                   acaATACGGGACCAATCCAGTATGTTCA‑TTCTCGGGT‑ATGGGGCCAGCAGTACTTTTGGAAAGGA‑TTCAGACATCGGGATGGGCGTAAGAGCGAAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                                                                                    <  1:320489/251‑1 (MQ=255)
                                                                                                                                                                ttCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGGTGAGTGCggg                                                                        >  2:342682/1‑250 (MQ=255)
                                                                                                                                                                          ataGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTCGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGCAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAATTCCCAGACGGAGACGAACTGTGACTGAACATGGTTGGTGGGAGCGGTTTactggt                                                                >  2:154670/1‑245 (MQ=255)
                                                                                                                                                                           aaGTGCCAGCAGTTCTTTTGGCAAGAA‑TTCAGACATCGGGATGGGCGAAAGAGCGAAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTAc                                                              <  1:342682/250‑1 (MQ=255)
                                                                                                                                                                             gTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGc                                                            >  1:62188/1‑251 (MQ=255)
                                                                                                                                                                                ccAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCAGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGGGGAGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGGAGACGAACTGGGTCTGACCATGGTTGGTGAGTGCGGTTTAATGCTACGCCac                                                         >  2:185982/1‑250 (MQ=255)
                                                                                                                                                                                           tttGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGATTGCGTGAAGTGGCGCGTGATGCGGACATTGCTCCCACCTCTTTTTATCGGCATTTCCACGACGTAGACGAGCTGAGTCTGACCCTTGTTGTTGAGAGCGGTTTAATGCTACGCCAACTCATGCGccc                                              >  2:255523/1‑250 (MQ=255)
                                                                                                                                                                                              ggCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGgc                                           >  1:381604/1‑251 (MQ=255)
                                                                                                                                                                                                 aaGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGcgcg                                        >  1:290334/1‑251 (MQ=255)
                                                                                                                                                                                                  aGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCAACTCTTTTTATCGACATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGGTTAATGCTACGCCACCTCATGCGcccggagcat                                       >  2:399823/1‑243 (MQ=255)
                                                                                                                                                                                                           aacatcgtgtggggcttaagacgaaaacaaAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTAt                               <  1:315062/226‑1 (MQ=255)
                                                                                                                                                                                                            aCATCGTGATGGGCGTACGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCAACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGGCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCCACTCATGCGCCAGGCGCGGCCGCGTATc                              >  2:198942/1‑251 (MQ=255)
                                                                                                                                                                                                            aCATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGTTTAGAGCGGTTTACTGATACGCCAACTCATGCGCCAGGCGCGTCAGCGTATc                              >  2:163051/1‑251 (MQ=255)
                                                                                                                                                                                                                   gATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCac                         >  1:362348/1‑248 (MQ=255)
                                                                                                                                                                                                                    atgtgcataagagagtaaaaaaaaGAAATAACCCGCCGTTCGCTGGTGGAATCCTCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGTCATTTCCGCTACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACGCATGCGCCAGGCGCGTCAGCGTATCGCCAAAgg                      <  1:34055/232‑1 (MQ=255)
                                                                                                                                                                                                                             agagCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCTCCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTTGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAgtgt             <  1:282879/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGTTGAGAGCGGTTTAATGCTACGCCAACTCCTGCGCCCGGAGCGTCAGCGTATCGCCAAAGGCGGGAGGGTGCTCCGCAccc  >  2:87349/1‑250 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
ACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCA‑TTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCT  >  NC_000913/4160916‑4161394

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