Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1525,1320TA36.7% 46.8 / 11.8 30intergenic (‑141/+160)unknown/unknownInositol transport system sugar‑binding protein/Myo‑inositol 2‑dehydrogenase 2 (EC 1.1.1.18)
Reads supporting (aligned to +/- strand):  ref base T (16/3);  new base A (0/11);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 6.66e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGTCTATCCAACGCTTTTGTTCACATGAAAACGCAATAGCATCCAACACGCGGGAAACCTT  >  1/524926‑525369
                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
aaagagccaggCTTGAGACTTTAATAGGATCCCCCCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGTTAAAACAAATAACAACAAAAATTAAAATACTGCTATCCCCCACCACAACAATCATATGCAAATATAATTTTCACCAATAAACATTTT‑AAAACAAACCTTCCGTT‑GAGTATCTATATCATAAAAAGAAGTTTTCCTGAAACCGCCATTTTTTGAAAGTGTGagat                                                                                                                                                                                                   <  1:288336/241‑1 (MQ=255)
  aaaTACCAGTTTGAGCCGTTAATAGGTTTCCCCCAAAAAGTAACGAAAGTAGTGTTTTTTTCATATTGCTCCTCCTTTTTATCGTTTAAACAAATAACAACAAAAATTAAAATACTCCTATCCCCCCCCAAAACAATCATATGCGAATAAAAATTTCCCCAATATAAAATTT‑AAAACAAATCGTCCGTT‑GATTATCTATATAATAAAAAAAAGTTTTCCTGCAACCGCCATTTTTTGAAAGTGTGagat                                                                                                                                                                                                   <  1:868719/248‑1 (MQ=255)
          gCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTAtt                                                                                                                                                                                         >  1:222559/1‑251 (MQ=255)
                               cccccacAAAGAAACGCAAGTATTTTCTTTTTAATAATGCTCCTCCCTTTAATCGATTAAACAAAAAAAAACAAAAATTAAAATTCTGCTATACCCCGCCACAACAATCATATGCGAAAATAAATTTCACCAATATA‑AATTTGGAAAAAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAAATTTCCTACAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAact                                                                                                                                                                    <  1:9863/248‑1 (MQ=255)
                               ccaccaAAAAGAAACGAAAGTAGTTTTTTTTTCATAATGCTCCCCCTTTTTATCGGTTAAACGAAAAAAAAAAAAAATAAAAATACTGCTATACCCCCCCACAACAATCATATGAGAATATAAATTCCACAAAAAAACAATTT‑GAAACAAATCGTCCGTT‑GAGTATTTATATAATAAAAAAAAGATTTCCGGAAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAact                                                                                                                                                                    <  1:540316/251‑1 (MQ=255)
                                       aaGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCaaa                                                                                                                                                            >  1:522467/1‑251 (MQ=255)
                                       aaGTAACCAAAGTGGTGTCTTTTTCATATTGCTACTCCCTTTTATCGTTTAAACGAAAAACAAAAAAAATAAAAATTCTGCTATCCCCCCCCACAACAATCATAGGCGAATATAAATTTCACAAATATACAATTA‑AAAACAAATGGCCCGGT‑GATTATCAATATCATAAAAAAAAGTTTCCCGGCAACCGCCATTTTTTGAAAGTGTGagat                                                                                                                                                                                                   <  1:187939/212‑1 (MQ=255)
                                         gTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTAttt                                                                                                                                                                                        >  1:755433/1‑221 (MQ=255)
                                                 agtagtGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTaa                                                                                                                                                  >  1:144002/1‑251 (MQ=255)
                                                         cTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGcc                                                                                                                                                                                                                       <  1:582233/174‑1 (MQ=255)
                                                         cTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAAtatttta                                                                                                                                          >  1:78656/1‑251 (MQ=255)
                                                                                             gaaaaacaaaAAAAATAAAAATACTGCAACCCCCCCCAACAACAAACATAGGCAAATATAATTTTCCCAAATATACATTTA‑GAAAAAAATC‑TCCCTTAGATTATCTATATAATAAAAAGAATATTCCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTa                                                                                                                                                                       <  1:539900/182‑1 (MQ=255)
                                                                                                       aaaaaTTAAAATACTGCTATCCCCCACCACAACAATAATATGCAAAAATAAATTCCACAAATATAAAATTT‑AAAAAAAATCGTCCGTT‑GAGTATCTATATCATAAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTc                                                                                            <  1:373157/251‑1 (MQ=255)
                                                                                                          aaTTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTAttt                                                                                                                                                                                        >  1:674773/1‑156 (MQ=255)
                                                                                                               aaaTACTTCTATCCCCCACCACAAAAATCATAGGCAAATATAAATTTCACCAATATACAATTT‑GAAAAAAATCGTCCTTT‑GAGTATCTATATCATAAAAAGAAGTTTTCCTGAAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGaa                                                                                    <  1:812807/251‑1 (MQ=255)
                                                                                                                    cTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGagatag                                                                                                                                                                                                 >  1:232114/1‑137 (MQ=255)
                                                                                                                    cTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATc                                                                                                                                                                                              >  1:719460/1‑140 (MQ=255)
                                                                                                                      gCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTa                                                                             >  1:867489/1‑251 (MQ=255)
                                                                                                                      gCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTa                                                                             >  1:846799/1‑251 (MQ=255)
                                                                                                                                    caacaaTCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTa                                                               >  1:798291/1‑251 (MQ=255)
                                                                                                                                    caacaaCCATATGCAAATATAAATTTCCCCAATATACAATTT‑AAAAAAAATCGTCCGTT‑GATTATCTATATAATAAAAAGAAGATTCCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTc                                                                                            <  1:357227/222‑1 (MQ=255)
                                                                                                                                       aaaTCATATGAGAATATAAATTTCACAAATATACAATTT‑GAAACAAATCGTCCGTG‑GAGTATCTATATCATAAAAAGAAGTTTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGt                                                            <  1:792644/250‑1 (MQ=255)
                                                                                                                                                              cACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGTCTATCCAACGc                                                 >  1:753846/1‑239 (MQ=255)
                                                                                                                                                                ccAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGTCTATcc                                                      >  1:129333/1‑232 (MQ=255)
                                                                                                                                                                ccAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGTCTATCCAACGCTTTTGTTCACATGa                                   >  1:32871/1‑251 (MQ=255)
                                                                                                                                                                caaaTATACAATTT‑GAAACAAATCGTCCTTT‑GAGTATCTATATCATAAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTc                                                                                                  <  1:788322/186‑1 (MQ=255)
                                                                                                                                                                          tttt‑GAAACAAATCGTCCGTT‑GAGTATCTATATAATTAAAAGAAGTTTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTAttt                                                                                                                                                                                        <  1:639233/91‑1 (MQ=255)
                                                                                                                                                                               aaaaaaaaTCGTCCGTT‑GAGTATATATATCATAAAAAAAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGTCTATCCAACGCTTTTGTTCACATGAAAACGCAATAGCAt                     <  1:501722/246‑1 (MQ=255)
                                                                                                                                                                                                   gTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGTCTATCCAACGCTTTTGTTCACATGAAAACGCAATAGCATCCAACACGCGGGAAACCtt  >  1:427275/1‑251 (MQ=255)
                                                                                                                                                                                                     aTCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAAt                                                                                                                                                         >  1:245199/1‑98 (MQ=255)
                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
AGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTT‑GAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGTCTATCCAACGCTTTTGTTCACATGAAAACGCAATAGCATCCAACACGCGGGAAACCTT  >  1/524926‑525369

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