Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,058,7310GT24.4% 68.0 / 13.8 41D78E (GAC→GAA) proP_7Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base G (12/19);  new base T (10/0);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 6.59e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGTTTCATCAACGTTTGGGAAAAACAACTTAGCGATGGTCGTAGCAA  >  ACB122/5058605‑5058864
                                                                                                                               |                                                                                                                                      
cAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGGGTTCGCCCCCCTTTTTGACCGA‑TGc                                                                                                                             >  2:297514/1‑137 (MQ=255)
      gACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGGGTTCGCCCCCGTTTGTGACCCC‑TGCGACCaa                                                                                                                       >  6:2886/1‑137 (MQ=255)
           gCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACg                                                                                                                  <  5:100706/137‑1 (MQ=255)
              gAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCGTGCGCCCCCGTTTGTCACCGA‑TGCGCCCAATTACCTaa                                                                                                               >  4:205273/1‑136 (MQ=255)
                aTTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCt                                                                                                             <  7:244547/137‑1 (MQ=255)
                    gAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTAACCGA‑TGCGACCAATTAAGTAGCTAcca                                                                                                         >  2:226003/1‑136 (MQ=255)
                       cACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGc                                                                                                      <  5:35755/137‑1 (MQ=255)
                       cACCAACGCCAATAGTCTCGTATGTTGGAATCAGAACAATACCAATACTGCCGTTTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGc                                                                                                      <  7:94214/137‑1 (MQ=255)
                        aCCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCTCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCg                                                                                                     <  6:42162/137‑1 (MQ=255)
                         ccAACGGCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAATTTATCGTCAGCAGGATCGC‑GGGCTTGCGCCGCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCgg                                                                                                    <  8:237133/137‑1 (MQ=255)
                           aaCGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTc                                                                                                  <  1:226003/137‑1 (MQ=255)
                            aCGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCGTGCGGCCCCGTTTTTGACCGC‑TGCGGCCCAATAAGTAGCTACCAAGCGGTCg                                                                                                 >  1:207834/1‑137 (MQ=255)
                               cctagaGTCTCGTATGTTGGAATCAGACCAATACCACTCCTGCCGATTGCCATGCAATTAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGGCACCGA‑TGCGACCCATTACGTAGCTACCTAGCGGTCGGGc                                                                                              <  8:294643/132‑1 (MQ=255)
                               ccAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGc                                                                                              <  3:205273/137‑1 (MQ=255)
                                                  gAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCGTGGGGCCCCCTTTTTGACCCAGTGCCCCCAATAAC‑TATCTAACAAACCGGCGGGCTATTAAAcccccccccgaa                                                                           >  5:199506/1‑133 (MQ=255)
                                                   aaTCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGGGCCCCCGTTTGTCACCGA‑TGCGCCCAATTACGTAGCTACCAAGCGGTCGGGCTATGAAACCCCCCCCGTAag                                                                          >  8:257335/1‑137 (MQ=255)
                                                   aaTCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGCCCAATTAAGTAGCTACCTAGCGGTCGGGCTATGAGACCCGCCCCGTaac                                                                          >  8:230783/1‑136 (MQ=255)
                                                   aaTCAGACCAATACCAATACTGCCGATTGCCAAGCAAATAATCGTCAGCAGCATCGCGGGGCTGCCGCCCCCTTTTTCCCCCGA‑T‑CGACCAAATAAGTTAATACCAAAAGGTAGGGCGATTAAACCCAcccccgaaa                                                                          >  4:41248/1‑132 (MQ=255)
                                                       aGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCGTGGGCCCCCTTTTTTGACCGC‑TGCGCCCAAATACCTAGGGACCTAGCGGTCGGGCGATGAAAACCGCCCCCTAagcagc                                                                      >  6:98678/1‑137 (MQ=255)
                                                           caataccaatacTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCCCCCCCTTTTTTCACCCC‑TGCGACCCATTACTTAGCTACCTAGCGGTAGGGCTATTAAACCCCCCCCCGAAAAAGCGaaa                                                                  >  8:304725/1‑136 (MQ=255)
                                                            aataccaatacTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGGGTGCGCCCCCATTTTTGAACGC‑TGCGGCCAATTAAGTTGCTACCAATCGGTCGGGGGATGAAACCCCGCCCCTAAACCGCGCAcc                                                                 >  6:33201/1‑136 (MQ=255)
                                                            aataccaatacTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCGTTTTTCAACGC‑TGCGACCAATTAAGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGCAcc                                                                 >  7:166622/1‑136 (MQ=255)
                                                                ccaatacTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGGGCCCCCGTTTGTCACCGA‑TGCGACCAATTTAGTAGCTAACTAGCGGTAGGGCTATGAAACCCGCCCCGTAAGCACCGAACGtgtg                                                             >  7:216897/1‑137 (MQ=255)
                                                                 caatacTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGc                                                            <  1:203838/137‑1 (MQ=255)
                                                                 caatacTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCATTGGCCCCCGTTTTTCAAAGA‑TGCGGCCAATTACGTAGCTACCAAGCCGTCGGGCTATGAAACACGCCCCGGAACCCCCGAACGACTGc                                                            >  6:144504/1‑137 (MQ=255)
                                                                      cTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGGCTGCGCCCCCATTTTTGCCCCC‑TGCGCCCAATTTAGTAGCTACCTAGCCGTCGGGGTATTCAACCCCCCCCGGAAGCAGCGAACCTCTGcaacaa                                                       >  4:39598/1‑137 (MQ=255)
                                                                          cGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGa                                                                                                                                 <  5:335740/60‑1 (MQ=255)
                                                                          cGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGa                                                                                                                                 >  6:335740/1‑60 (MQ=255)
                                                                            aTTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGt                                                 <  3:289127/137‑1 (MQ=255)
                                                                            aTTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCGTGCGCCCCCGTTTGTCACCGA‑TGCGCCCAATTTAGTAGCTACCTAGCCGTCGGGCGATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGt                                                 >  4:126762/1‑137 (MQ=255)
                                                                                                 agcagcATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGTTTCATCAACGTTTGGGaa                            <  4:221961/137‑1 (MQ=255)
                                                                                                 agcagcATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGTTTCATCAACGTTTGGGaa                            <  2:86143/137‑1 (MQ=255)
                                                                                                  gcagcaTCGC‑GGGCTTGCGCCCCCGTTTGTCACCCA‑TGCGCCCAATTACGTAGCTACCTAGCGGGCGGGCTATGAACCCCCCCCCGTAAGCAGCGCACGGCTGCCACAAGGATGTCGTTTCATACACGGTTGGGaaa                           >  2:213418/1‑137 (MQ=255)
                                                                                                                  tGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGTTTCATCAACGTTTGGGAAAAACAACTTAGCGATg            <  6:199506/137‑1 (MQ=255)
                                                                                                                     gcCCCCATTTGTCACCCA‑TGCGCCCAATTACGTAGCTAACTAGCGGTCGGGCTATGAAACCCGCCCCGTAATCAGCGAACGTCTGGAACAAGGATTTCGTTTCATAAACGTTTTGGAAAAACAACATAGCGATGGtt         >  8:328582/1‑136 (MQ=255)
                                                                                                                           ttttGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGTTTCATCAACGTTTGGGAAAAACAACTTAGCGATGGTCGTAGCa   <  6:339631/137‑1 (MQ=255)
                                                                                                                           ttttGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGTTTCATCAACGTTTGGGAAAAACAACTTAGCGATGGTCGTAGCa   <  6:9578/137‑1 (MQ=255)
                                                                                                                           ttttGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGTTTCATCAACGTTTGGGAAAAACAACTTAGCGATGGTCGTAGCa   <  5:33201/137‑1 (MQ=255)
                                                                                                                           ttttGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGTTTCATCAACGTTTGGGAAAAACAACTTAGCGATGGTCGTAGCa   <  4:155567/137‑1 (MQ=255)
                                                                                                                            tttGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGTTTCATCAACGTTTGGGAAAAACAACTTAGCGATGGTCGTAGCaa  >  4:131832/1‑137 (MQ=255)
                                                                                                                            tttGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGGTTCATCTCCGTTTGGGATATGCATCTTAGCGTTGGTCGGAGctt  >  5:179066/1‑135 (MQ=255)
                                                                                                                               |                                                                                                                                      
CAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGC‑GGGCTTGCGCCCCCTTTTGTCACCGA‑TGCGACCAATTACGTAGCTACCTAGCGGTCGGGCTATGAAACCCGCCCCGTAAGCAGCGAACGTCTGCAACAAGGATGTCGTTTCATCAACGTTTGGGAAAAACAACTTAGCGATGGTCGTAGCAA  >  ACB122/5058605‑5058864

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