Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,613,0530CT70.3% 44.9 / 17.4 38A385V (GCC→GTC) hyfRDNA‑binding transcriptional activator HyfR
Reads supporting (aligned to +/- strand):  ref base C (8/3);  new base T (22/4);  total (31/7)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGCATCGCTATTGCCGTTGGCAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGAGCGGG  >  NC_000913/2612814‑2613298
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
tcgcatcgcTATTGCCGTTGGCAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATcc                                                                                                                                                                                                                                             >  1:393163/1‑251 (MQ=255)
     tcgcTATTGCCGTTGGCAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGCAAACCACCAGATTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGAGATCCATCa                                                                                                                                                                                                                                         >  1:615621/1‑250 (MQ=255)
      cgcTATTGCCGTTGGCAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAAc                                                                                                                                                                                                                                       >  1:14282/1‑251 (MQ=255)
                       aTGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCTATCCATCAACTTAGCCCgcgacgcgac                                                                                                                                                                                                                     <  2:572098/251‑1 (MQ=255)
                            gATGCCTGGCGTAGCATGACCGATTTGCAGGAACGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGcc                                                                                                                                                                                                                >  1:403698/1‑251 (MQ=255)
                            gATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGTTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGcc                                                                                                                                                                                                                <  1:247062/251‑1 (MQ=255)
                                      gTAGCATGACCGATTTGCAGGAAAGTTTGCAGCCAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGACTCTGGGCATCGGTGACATTATCTATCAAAGACAGGCAATGGAAGACCAGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGCTTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATAGCCAGAGCGATCCATCCACTTAGCCCGCGACGCGACACGCCGCTGGTCaa                                                                                                                                                                                                       >  1:367313/1‑250 (MQ=255)
                                       tAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGACAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGGGGTGAAACCGGAACTGGCAAAGAGGGGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCaca                                                                                                                                                                                                      >  2:238622/1‑248 (MQ=255)
                                           aTGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGCGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGCGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACCAGCCGCTGGTCAAAATCa                                                                                                                                                                                                  >  2:529088/1‑250 (MQ=255)
                                                 gATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTgcgc                                                                                                                                                                                           >  1:81389/1‑251 (MQ=255)
                                                    ttGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTg                                                                                                                                                                                         <  1:300680/250‑1 (MQ=255)
                                                             aGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGCTTTGTGGTGAACCCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGc                                                                                                                                                                               >  1:434304/1‑251 (MQ=255)
                                                                     gcaAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGtcttc                                                                                                                                                                       >  1:23923/1‑251 (MQ=255)
                                                                           aaaCCACCCGCTTAGCGATCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACATGCTACAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAGACCGGAACTGGCAAAGTGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGCCAAGCCGCTGGGCTAAATCAACTGCGCTGCCAGCCCCGCCCGTCTTCTggcc                                                                                                                                                                  >  2:143477/1‑248 (MQ=255)
                                                                                 ccaGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGt                                                                                                                                                           >  2:192035/1‑251 (MQ=255)
                                                                                  caGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTCGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGtt                                                                                                                                                          >  2:470838/1‑251 (MQ=255)
                                                                                        aGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCCATGGAAGACCTGCTCCAGCAGGCAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCACAATCAACTGCGCTGCCATCCCCGCCCGTCTTCTGGAAAGTGAGTTATTCgg                                                                                                                                                    >  2:305526/1‑251 (MQ=255)
                                                                                                  tCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAACCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGAGAGTTATTCGGTCATGACac                                                                                                                                           >  1:107474/1‑249 (MQ=255)
                                                                                                                     cGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCg                                                                                                                        >  1:477605/1‑250 (MQ=255)
                                                                                                                         gACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTa                                                                                                                    >  2:97118/1‑250 (MQ=255)
                                                                                                                         aatattataaaaaaaaGACAGGAAATGGAAGACCTTCCCCAGCAGATAGATATTGGGGCAAAGAGCGACAGTGGGGTGTTGATTGGGGGTAAACCCGGAAGTGAAAAAGAGGTGATCGTCAGAGCACTCCATCAACTTGCCCCGCACCGCGACATGCCGCGGGTCAAAATCAACGGCGCTGCCACCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTaa                                                                                                                   <  1:589187/238‑1 (MQ=255)
                                                                                                                          aCATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACGGGTGCGATTAAt                                                                                                                  >  1:162444/1‑251 (MQ=255)
                                                                                                                             ttATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATAcc                                                                                                               >  2:18140/1‑251 (MQ=255)
                                                                                                                                        aGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCg                                                                                                    >  1:322314/1‑251 (MQ=255)
                                                                                                                                            aGGCAATGGAAGACCTGCTCCAGAAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAACGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGtttt                                                                                                >  2:86923/1‑251 (MQ=255)
                                                                                                                                                aaTGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGaa                                                                                             >  1:103546/1‑250 (MQ=255)
                                                                                                                                                  tGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGCGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTCCTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAAtt                                                                                          >  2:594620/1‑251 (MQ=255)
                                                                                                                                                     aacaccgtccacagcagGGAGATTTTGGGGCGAAGAGCGACAGGGTGGTGTTGATTTGTGGAAAACCCGGAACGGGAAAAGGGGTGATCGTCGGCGCGATCCACAAACTTCCCCCGCGCCGCAAAACGCCGCTGGTCAAAATCACCTGCGCTGCCACCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGcc                                                                                       <  1:201365/240‑1 (MQ=255)
                                                                                                                                                      gaccctgctcaggcgggTAGATTTTGTGGCAAAGAGCAACAGTGGGTTTTTAATTTGTGGAAAAACCGAAACGGGCAAGGAGGTGATCGCAAGAGCAACCCAAAAACTAAGCCCGCGACGCACCAAGCCGCTGGTAAAAAACAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCg                                                                                      <  1:104752/236‑1 (MQ=255)
                                                                                                                                                               ccAGCAGGGAGATATTGTGGAGAAGAGCGACAGTACGATGTTGATTTGTGGTGAAAACGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGCTTAATACCCATCGTGGTCGCTTTGAAATTGCCGATggcggc                                                                              >  2:106767/1‑250 (MQ=255)
                                                                                                                                                                      gTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACgttgtt                                                                      >  1:283941/1‑251 (MQ=255)
                                                                                                                                                                                                 cgggtggtGATTGGTGGTAAACCCGAACCTGGAAAAGGGGTGATCGCCAGAGCGACCCACCACTTAAGCCCGCGACGCGACAAGCCGCTGGTAAAAATCAACGGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGtt                                           <  1:475537/249‑1 (MQ=255)
                                                                                                                                                                                                         gATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGACGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATGGGCGACCTGCCGTTAGACCTTc                                   >  2:468572/1‑251 (MQ=255)
                                                                                                                                                                                                          aTTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGATTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTc                                   >  2:283587/1‑250 (MQ=255)
                                                                                                                                                                                                                                   aaaGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCCACTGCGCTGCCCTCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGCAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCAGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGc           >  1:253715/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aaGAGGTGGTCGACAGATCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCCTCCCCGCCAGTCTTCTGGCAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATggcggaa                                                                             >  2:330892/1‑181 (MQ=255)
                                                                                                                                                                                                                                    aaGAGGTGATCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGc           >  2:497954/1‑249 (MQ=255)
                                                                                                                                                                                                                                            aTCGTCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGAGCggg  >  1:423265/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
TCGCATCGCTATTGCCGTTGGCAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGAGCGGG  >  NC_000913/2612814‑2613298

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