Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,665,6980TC24.0% 104.1 / 12.4 50D38G (GAT→GGT) gadXacid resistance regulon transcriptional activator; autoactivator
Reads supporting (aligned to +/- strand):  ref base T (30/8);  new base C (0/12);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 1.04e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGA  >  NC_000913/3665493‑3665920
                                                                                                                                                                                                              |                                                                                                                                                                                                                              
aaGATTATTATTCTGAGATAGTTGCGCAACTCCCTAAAAAACAAGAATATTTCTACAACAGTGATGAACGAGTAATGAAAAAACTTTTTTATGTGCATGAAGATTTAATCACCCCCCCTCGAGCATCGGGAAAAATAATAAAAGCGTGCCCCTTCACACAAAAAAAAAATTTTCAACCCCCTTATCTCCTCAATTTTGGTTTCCATCCAAAAAACCCAGGTCACCGCCATTAAAATAg                                                                                                                                                                                                <  2:699139/238‑1 (MQ=255)
 agtttCTTATTCTGCGATAGTTGCGCAACTTCCTAAAAAACAGGAATATCTCTACTAAAGTGATGAACGAGTAATGAAGAAATTTTTTTATGTGCATAAAGTTTTAATGCCCCCTCCTTGAGCATGGGGAAAAAAAATAAAAGCGTGTCCCTTGACACGAATAAAAATTTTTAAACACCCTTATCTACCCGATTTGGGCTTGCACCCGCAAAAACAAGGTCACCGCCATTAAAATAGCg                                                                                                                                                                                              <  1:726002/236‑1 (MQ=255)
               cGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCAt                                                                                                                                                                                                                                <  1:797793/191‑1 (MQ=255)
               cGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCAt                                                                                                                                                                                                                                >  2:797793/1‑191 (MQ=255)
               cGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTa                                                                                                                                                                                 >  1:164901/1‑238 (MQ=255)
               cGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTa                                                                                                                                                                                 >  1:523855/1‑238 (MQ=255)
               cGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTa                                                                                                                                                                                 >  1:36385/1‑238 (MQ=255)
                          aaacttctaaaaaaaaAGAAATATCTCTAATACAGTGATGAACGAGTAATAAAAAAACTTTTTTATGTGAATAAATATTTAATGCCCCCTCCTTGAACATCGGGAAAAATAATGAAAGCGTGTCCCTTGACACGAAAAAAAAATTTTAAACACACTTATCTACTAAATTTTGTCTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTgaga                                                                                                                                                                     <  2:36385/229‑1 (MQ=255)
                           atctcccagAAAAACAGAAATAGCTCAACTACAGGAATGAAGGAGTAATGAAAAATCTTTTTTGGGAGCAGAAATATATAATCCCTCCCCCTTGAGCAGCGGGAAAAAAAATAAAAGGGTGTCCTTTGACACGAACACAAAATTTACAACACCCTTATCTCCACAAATTTGGTTTGCTCCCACAAAAACCGGGCCACCGCCATTAAAATAGCGATATTCCCCATTAACCATGGTGAGaa                                                                                                                                                                    <  2:625262/231‑1 (MQ=255)
                                                   tatactacAGGAATGAACAAGTAAGGAAAAAACTTTTTTATGTGAATAAAGTTTAAATGCCCCCTCCTTGAGCAGCGGAAAAAAAAATAAAAGCGTGTCCCTTGCAACAAAAACAAAATTTCCAACACACTTATCTCCTAAATTTGGCTTTGCACCCGAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                                                                                              <  2:77978/203‑1 (MQ=255)
                                                     tacaaaaGTGATGAACGATTAATGAAAAATTTTTTTTATTTGAATAAAGTTTTAATGCCTCCTCCTAGAGAATGGGAAAAAATAATGAAGGCGTTCCCCTTAACACAAAAAAAAAATTTCCAACCCACTTATCTACACAATTTTGGTTTCCACCCGAAAAACCCGGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGt                                                                                                                                                                         <  2:375563/202‑1 (MQ=255)
                                                          cAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAc                                                                                                                                                                               >  1:339747/1‑197 (MQ=255)
                                                           aGTGAGAAACGAGTAATAAAGAATCTTTTTTATGTGAATAAAGATTTACCCCCTCCCCCTGGAGCATGGGGAAAAATAATGAAAGCGTGTCCCTTAACACAAATAAAAAATTTTCAACACACTTATCTACTCGATTTAGGCTTTCATCCGAAAAAACCAGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACaa                                                                                                                                    <  1:683810/239‑1 (MQ=255)
                                                             tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAACTTTTTAAAACACTTATCTACTCGAAATTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTAAGTCTTGCATACGCAATTATACAAtt                                                                                                                                  >  2:40005/1‑239 (MQ=255)
                                                             tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                                                                                                                                                    >  1:77978/1‑205 (MQ=255)
                                                             tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAAt                                                                                                                                   >  1:851732/1‑238 (MQ=255)
                                                               atgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatatat                                                                                                                                                               >  1:375563/1‑208 (MQ=255)
                                                               atgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTc                                                                                                                                 >  1:603730/1‑238 (MQ=255)
                                                               aaaaaCGATTAATGAAAAAATTTTTTTATGTGCATGAATTTTTAATGCCTCCTCCTGGAGAATCGGGAAAAAAAATGAAAGCGTGCCCCTAAACACAAATACAAAATTTTAAACACATTTATCTTCTTAATTTTGTCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTcc                                                                                                                                <  2:523855/236‑1 (MQ=255)
                                                                         aaTGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAACTTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATg                                                                                                                            >  1:706080/1‑233 (MQ=255)
                                                                                aattttttttATGTAAATAAATTTTTAATCCCTCCTCTGTAGAAATGGGAAAAAATAATAAAATGGTTTCCTTGAACAGGAAAAAAAAATTTCAAACACTTTTATCTACTAGAATTTGGCTTCCACCCGAAAAACCCGGGACACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATa                                                                                                               <  2:851732/235‑1 (MQ=255)
                                                                                aaaCTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATg                                                                                                                                                          >  1:792733/1‑196 (MQ=255)
                                                                                            gtgtaagAAGATTTAATCCCCCCCCCGTAAAAAGCGGGAAAAATAATGAAAGCTTGTCCCTTGACACGAATACAAAATTTCCACCACACTTATCTACTCAATTTTGGCTTCCATCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAata                                                                                                   <  1:89628/233‑1 (MQ=255)
                                                                                                          aaTGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGt                                                                                                                                                                         >  1:476392/1‑155 (MQ=255)
                                                                                                            agCCCCCTCCTGGAGGTTGGGGAAAAATAAGAAAGGCGTGTCCCTCAACCCGAATAAAAAATTTTCAACACTTTTCTCTAGTCGAATTTGTTTGTCATCCCAAAAAACCAGGTTCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTa                                                                                   <  1:601752/238‑1 (MQ=255)
                                                                                                                tcctcctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATAtt                                                                                    >  1:566101/1‑234 (MQ=255)
                                                                                                                  ctcctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGa                                                                                                                                       >  1:267537/1‑181 (MQ=255)
                                                                                                                  ctcctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTa                                                                                   >  1:206665/1‑233 (MQ=255)
                                                                                                                   cctttGAGCTGCGGAAAAAATAATGAAAGCGTGTCCCTTGACAAAAATAAAAAATTTCAAACACATTTATCTCCTCAATTTTGGTTTCCACCCGAAAAACCCGGTTAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                                                                                                                                                    <  2:790955/148‑1 (MQ=255)
                                                                                                                            aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTa                                                                    >  2:153692/1‑238 (MQ=255)
                                                                                                                             tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                                                                                                                                                                   >  1:641413/1‑110 (MQ=255)
                                                                                                                             tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATg                                                                                                                                                          >  1:790955/1‑151 (MQ=255)
                                                                                                                                   aaaaaaaTGAAAGCTTCTCCTTTGCCCCAAATACAAAATTTAAACCACCCTTATCGACACGATTTTGGCTGGCACCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCAcgccg                                                            <  1:291660/238‑1 (MQ=255)
                                                                                                                                          tGAAAGCTTGTCCTTTGACACAAATAAAAAATTTTCAACCCCCTTATCTCCTAGATTTTGTCTTCCACCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGa                                                                                                                      <  2:231129/174‑1 (MQ=255)
                                                                                                                                              aGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGTTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTa                                                                                   <  1:410103/205‑1 (MQ=255)
                                                                                                                                              aGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCAt                                                                                                                                                                                                                                <  1:594204/64‑1 (MQ=255)
                                                                                                                                              aGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCAt                                                                                                                                                                                                                                >  2:594204/1‑64 (MQ=255)
                                                                                                                                              aGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCAAGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTa                                                                                   >  2:410103/1‑205 (MQ=255)
                                                                                                                                                    tCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACTATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGa                                                                                                                                       >  1:620310/1‑147 (MQ=255)
                                                                                                                                                    tCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGtt                                                                                                            >  1:231129/1‑174 (MQ=255)
                                                                                                                                                                attACAAAATTTAAAACCCCTTTTTCTACTGGAATTTGGTTTCCACCCAAAAAACCCAGGCCGCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTga                                                                                                                                                                       <  1:4157/101‑1 (MQ=255)
                                                                                                                                                                          tttCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                                                                                                                                                    >  1:556821/1‑96 (MQ=255)
                                                                                                                                                                          tttCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAttt                                                                                                                                                             >  2:4157/1‑103 (MQ=255)
                                                                                                                                                                                     ttATCTACTCTATTTTGGCTTGCATCCACAAAACACAGCTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCTTGAGaa          <  2:566101/239‑1 (MQ=255)
                                                                                                                                                                                       aTCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGt                                                 <  1:174727/198‑1 (MQ=255)
                                                                                                                                                                                       aTCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGt                                                 >  2:174727/1‑198 (MQ=255)
                                                                                                                                                                                             tCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCg                                                                                                                                                                                              >  1:381012/1‑51 (MQ=255)
                                                                                                                                                                                             tCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACa                                                                                              >  1:732768/1‑147 (MQ=255)
                                                                                                                                                                                             tCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGgatga  >  1:213080/1‑239 (MQ=255)
                                                                                                                                                                                               gattttGGCTTGCACCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                                                                                              <  1:765483/62‑1 (MQ=255)
                                                                                                                                                                                               gAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                                                                                              >  2:765483/1‑65 (MQ=255)
                                                                                                                                                                                                 aTTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCa                                                                                                                                                                                                         >  2:685524/1‑36 (MQ=255)
                                                                                                                                                                                                 aTTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCa                                                                                                                                                                                                         <  1:685524/36‑1 (MQ=255)
                                                                                                                                                                                                    tGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTga                                                                                                                                                                       >  1:370147/1‑67 (MQ=255)
                                                                                                                                                                                                              |                                                                                                                                                                                                                              
AAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGA  >  NC_000913/3665493‑3665920

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