Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,6490AG61.5% 26.6 / 39.3 39W107R (TGG→CGG) hdfRDNA‑binding transcriptional dual regulator HdfR
Reads supporting (aligned to +/- strand):  ref base A (9/6);  new base G (13/11);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTGAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACT  >  NC_000913/3947406‑3947875
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                  
cccACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGAtt                                                                                                                                                                                                                              >  2:276147/1‑250 (MQ=255)
     acggcactgtggtAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCGGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCAt                                                                                                                                                                                                                        <  1:25132/239‑1 (MQ=255)
                    ttAAGATCTCCCTTTAGTTTTGCAGGGGCACTGGTATAACGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACgg                                                                                                                                                                                                         >  2:308/1‑250 (MQ=255)
                      aaGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAg                                                                                                                                                                                                        <  2:131973/250‑1 (MQ=255)
                       aGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGATAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGgc                                                                                                                                                                                                      >  1:458753/1‑251 (MQ=255)
                        gATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGgcg                                                                                                                                                                                                     >  1:257785/1‑251 (MQ=255)
                            tcCCTTTAGTTTTGAAGGGGCACTGGTATCAAGCGCTAAAGTGAACTATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTgg                                                                                                                                                                                                 >  2:466164/1‑251 (MQ=255)
                                      tttGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAg                                                                                                                                                                                        <  1:347160/250‑1 (MQ=255)
                                                    ggTATAAAGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAACTTCGTCCATTTTGGGCGCTTCAGTGGTCATAAGAAGATCAAGCTGGCGATCATGCCGCTGTTTTACCCGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTc                                                                                                                                                                         >  1:243038/1‑251 (MQ=255)
                                                                tAAAGTGAAATATCCCAGAAAATGACTACTAAATTTGTCCATTTTGCGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGCGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTACGCATACATTCCCACAACGAGGCGCTGGCCCCGATAGAAAACTCGTTATGTCGTGAGGTATgcg                                                                                                                                                             >  2:438872/1‑251 (MQ=255)
                                                                 aaaGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCCTTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGCTTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGAATgcgc                                                                                                                                                            >  2:234109/1‑251 (MQ=255)
                                                                        aaTATCCCAGCAACTGACTCCTCAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCAcctcc                                                                                                                                                     >  1:173037/1‑251 (MQ=255)
                                                                                   aaCTGACTACTAAATTCGTCCATTTTGGGCGCTTCCGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGAATGCGCCACCTCCTTACGggcggc                                                                                                                                          >  2:414376/1‑251 (MQ=255)
                                                                                     cTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCt                                                                                                                                        >  1:270639/1‑251 (MQ=255)
                                                                                      tGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTg                                                                                                                                       >  1:247686/1‑251 (MQ=255)
                                                                                        actactAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGcc                                                                                                                                     >  1:146056/1‑251 (MQ=255)
                                                                                               aaTTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCt                                                                                                                              >  1:268730/1‑251 (MQ=255)
                                                                                                    gTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTGTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGa                                                                                                                         <  1:169227/251‑1 (MQ=255)
                                                                                                       cATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGTGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTTCTGGCCGCCTGCCAGGTGCTCATGAGCg                                                                                                                      <  1:415921/251‑1 (MQ=255)
                                                                                                           ttGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTGTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCTGCATGCCAGGTGCTCATGAGCGttt                                                                                                                   <  1:272070/250‑1 (MQ=255)
                                                                                                               gcgcTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCa                                                                                                              <  2:135203/251‑1 (MQ=255)
                                                                                                                                 gaagaTCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTc                                                                                             >  1:408347/1‑250 (MQ=255)
                                                                                                                                   agaTCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAATCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGTCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGTGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCAcc                                                                                          <  1:276835/251‑1 (MQ=255)
                                                                                                                                   agaTCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCAcc                                                                                          <  1:133815/251‑1 (MQ=255)
                                                                                                                                         aGCTGGCGTTCATGCAGCTGTTTTGCCAGAGACTGCCGTTGGGCAAGTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGGTATGTCGTGAGGTATGCGCCACCTCCGTACGGGCGGGCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGtttt                                                                                               <  1:242355/240‑1 (MQ=255)
                                                                                                                                         aGCTGGCGTTCATGCAGCTGTTTTCCCAGAGACTGCCGTTGGGCAATTGGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGCTACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGCGCGTTTCTGCATAAGGCAGTAGttttagatcggcagg                                                                                    >  2:242355/1‑240 (MQ=255)
                                                                                                                                              gCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTa                                                                               <  1:184525/251‑1 (MQ=255)
                                                                                                                                                gTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTaa                                                                              >  2:405958/1‑250 (MQ=255)
                                                                                                                                                      gCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGat                                                                        >  1:364642/1‑250 (MQ=255)
                                                                                                                                                           tGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGcc                                                                                                                                         >  1:379578/1‑180 (MQ=255)
                                                                                                                                                           tGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGcc                                                                                                                                         <  2:379578/180‑1 (MQ=255)
                                                                                                                                                           tGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGtt                                                                  >  1:14667/1‑251 (MQ=255)
                                                                                                                                                            gTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGttt                                                                 <  2:152052/251‑1 (MQ=255)
                                                                                                                                                                  ccAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTGGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATg                                                                                                                          <  1:206775/188‑1 (MQ=255)
                                                                                                                                                                  ccAGAGACTCACGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGCTTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATg                                                                                                                          >  2:206775/1‑188 (MQ=255)
                                                                                                                                                                                      aTTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGGGAGGTAGGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTGCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTGAAAAGGTTCACAcc                                       <  1:466164/251‑1 (MQ=255)
                                                                                                                                                                                         cgcgCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTGAAAAGGTTCACACCCAg                                    <  1:276147/251‑1 (MQ=255)
                                                                                                                                                                                                                          tGATTTTGATACCAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGGGAAAGGGTTACACCCAGTTGGTTTTTCtgttgtttggtttgaaaagtccc   >  1:122279/1‑234 (MQ=255)
                                                                                                                                                                                                                           gATTTTGATACAAGCGTCCCAGCCGCTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTGAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCAct  <  2:122279/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                  
CCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTGAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACT  >  NC_000913/3947406‑3947875

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