Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0860AC53.2% 46.8 / 73.0 47Y580D (TAC→GAC) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (21/1);  new base C (0/25);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.75e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-10

TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGT  >  NC_000913/3584938‑3585227
                                                                                                                                                       |                                                                                                                                             
tattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcc                                                                                                                                             >  2:344739/1‑151 (MQ=255)
        aaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt                                                                                                                                     >  2:79569/1‑151 (MQ=255)
            cAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAat                                                                                                                                 >  1:36759/1‑151 (MQ=255)
              gtggtaAAATTCATTCTCGGGGAACGGCAAACCCGCTGAAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTAACCTTAGTCCCCCGCCGTTAAatc                                                                                                                                <  2:791535/150‑1 (MQ=255)
                         cattcTCTGGGAACGGCAAAACCGTTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGATACCTTCGGCTTGCCCTGACAAAATACCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAACCATCCACCGAgtg                                                                                                                    >  1:648159/1‑149 (MQ=255)
                          attcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAgcgc                                                                                                                   >  1:983242/1‑151 (MQ=255)
                           ttcCCTGGGAGCGCAAACCCCGCGGGAAGAATATCTACAGGCTTACAAAAAAAAAGAGGCTACCTCCGTCTCCCCCTAACAAAATCCCCCTCTCCCCGCAAGGGGGCCCGTAAACTTAAGCCCCCCGCCGTTAAATCATCCACCGAgc                                                                                                                     <  1:205738/148‑1 (MQ=255)
                               ctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTc                                                                                                              >  2:423025/1‑151 (MQ=255)
                                cgggAAAGGAAAAACCCGCGGAAAAAATATCTAAAGAATAAAAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCggg                                                                                                                <  1:573405/148‑1 (MQ=255)
                                 gggAACGGAAAACCCGCTGGAAGAATATCTAAAGGATTAAAAAAAAAAAAAGGCTTCCTCCGGCTTGCCCTGACAAAGTACCCCTCTCCCCCGGAGGGGGCCCGCTACCCTTTGTCCCCCGCCGTTAAATCATCCACCGAgcgc                                                                                                                   <  2:154880/144‑1 (MQ=255)
                                 gggAACGACAAACCCGCTGAAAGAATATCTACAGGATTACAAAATAAAAAAGGCTCCCTTCGGCTTGCCCTGACAAAATGCCCCTCTTCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                            <  2:901473/151‑1 (MQ=255)
                                  gggaaggaaaaaCCCGTGGAAAAATTATATACAGGTTAACAAAAAAAAAAGAGCTCCTTTCGCCTTGCCCTAACAAAATCGCCCTCTCCCCGGGAAGAGGGCCGTTACCCTTATACCCCCGCCGTTAAATCATCCACCGAGCGCGGGTc                                                                                                              <  1:1026047/141‑1 (MQ=255)
                                      cgccAAAACCGCGGAAAGAAAATCTACAGGATTACAAAAAAAAAAAGGCTCCCTCCGGCTTCCCCTAACAAAATACCCCTCCCCCCCCGAAGGGGGCCGCTACCCTTATTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGa                                                                                                           <  1:781448/144‑1 (MQ=255)
                                      cGGAAAACCCGCTGAAAAAATATCTACAGGATTAAAAAAAAAAAAAGCCTCCCTTCGGCTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                        <  2:676377/150‑1 (MQ=255)
                                             aCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTAc                                                                                                 >  2:1125968/1‑150 (MQ=255)
                                               ccGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                >  2:335131/1‑150 (MQ=255)
                                                      cagaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTACCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                        >  1:702448/2‑151 (MQ=255)
                                                      aaaaaTATCTACAGGATAAAAAAAAAAAAAGGGCTCCCTGCGGCTTCCCCTGAAAAAATACCCCTCCCCCCCGAAAGAGGGCCGCAACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                           <  1:937674/145‑1 (MQ=255)
                                                      aaaaaTATCAACAGAATTAAAAAAAAAAAAGGGCTCCCTCCTGCTTGCCCGGAAAAAATGCCCTTCTCCCCAAGAAGGGGGCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                <  1:901229/140‑1 (MQ=255)
                                                       aaaaTATCAACAGGTTAACAAAAAAAAAAAGCCTCCCTCCGCCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAGGGGGGCCGTTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGt                                                                                       <  1:79546/149‑1 (MQ=255)
                                                         aaTATCTACAGGATTACAAAAAAAAAAAGGCTCCCTCCGGCTTCCCCTGACAAAATACCCCTCTCCCCACGAGGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                     <  2:932724/151‑1 (MQ=255)
                                                            tatCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCgg                                                                                   >  2:926733/1‑151 (MQ=255)
                                                               cTACAGTATTACGAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTAAAAAAATAGCCCTCTCCCCACGAAGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                    <  2:568320/147‑1 (MQ=255)
                                                               cTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                >  2:30149/1‑151 (MQ=255)
                                                                tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                >  1:681382/1‑150 (MQ=255)
                                                                           aaaaaaaaaaaGCCTCCCTTCGGCTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                    <  1:504519/150‑1 (MQ=255)
                                                                            aaaaaaaaaaGGCTCCCTTCGGCTTCCCCTAAAAAAATACCCTTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATa                                                                         <  2:635816/145‑1 (MQ=255)
                                                                                 aaaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                                              >  1:161554/1‑151 (MQ=255)
                                                                                             ttCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                  >  1:1029377/1‑151 (MQ=255)
                                                                                                     gCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                           >  1:1188594/1‑150 (MQ=255)
                                                                                                        cTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgc                                        >  1:1124977/1‑150 (MQ=255)
                                                                                                        cTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgc                                        >  1:529412/1‑150 (MQ=255)
                                                                                                           aaaaaaTAGCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGcc                                     <  2:1118832/148‑1 (MQ=255)
                                                                                                             aaaaTCCCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACt                                  <  2:137389/151‑1 (MQ=255)
                                                                                                              aaaTAGCCCTCTTCCCACGAGGAGGCCCGCTACCCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACtt                                 <  2:1188858/151‑1 (MQ=255)
                                                                                                               aaTACCCCTCTCCCCACGAAGGGGGCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                <  2:102475/151‑1 (MQ=255)
                                                                                                                atcccccTCCCCCCACGAAGAGGGCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTc                               <  1:161513/147‑1 (MQ=255)
                                                                                                                  aGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                              >  1:458616/1‑150 (MQ=255)
                                                                                                                    cccTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAc                           >  1:375724/1‑151 (MQ=255)
                                                                                                                    cccTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTg                             <  2:1157148/149‑1 (MQ=255)
                                                                                                                          ccccACGAGGGGGCCCCCAAACTTTACTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtttt                      <  1:885455/149‑1 (MQ=255)
                                                                                                                          ccccACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtt                        <  2:75598/147‑1 (MQ=255)
                                                                                                                           cccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGc                    >  2:765371/1‑151 (MQ=255)
                                                                                                                                 agggggccggctccccTGAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCca               <  1:963425/136‑1 (MQ=255)
                                                                                                                                         cGCTACCTTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa      <  1:1173500/151‑1 (MQ=255)
                                                                                                                                         cGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa      <  2:693301/151‑1 (MQ=255)
                                                                                                                                              aCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt  >  1:242825/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGT  >  NC_000913/3584938‑3585227

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