Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,492,0280CT42.9% 43.0 / 16.7 37Q90Q (CAG→CAAppiApeptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base C (3/17);  new base T (15/0);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 2.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGGTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCG  >  NC_000913/3491879‑3492178
                                                                                                                                                      |                                                                                                                                                         
ccATGGTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTct                                                                                                                                                          <  2:2520866/151‑1 (MQ=255)
ccATGGTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTc                                                                                                                                                           <  1:2214032/150‑1 (MQ=255)
    ggTCAGGAAAGGCGTTATCGGGGGCAGAGAGAGACACCTTGCTGGTGACGCTGGCGTTGTCGGCGGTACGTGAAAAAAAGATGGTACCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGTTTTTTttgt                                                                                                                                                        >  1:333058/1‑148 (MQ=255)
    ggTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTctgct                                                                                                                                                       <  1:3574690/150‑1 (MQ=255)
       cAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCa                                                                                                                                                    <  1:2614457/150‑1 (MQ=255)
                cGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGTTTTTTTTGCTGCATCTg                                                                                                                                                >  2:2098828/1‑145 (MQ=255)
                   tATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGa                                                                                                                                       <  1:3905773/150‑1 (MQ=255)
                      cGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAAcc                                                                                                                                   <  1:1387627/151‑1 (MQ=255)
                        gCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTTGTTTTTTTTTCTTCCTCTGCTCGGGTGAAccgc                                                                                                                                 >  2:2953116/1‑151 (MQ=255)
                              ttGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTTGGTTTTTTTTCTGCCTCTGCCCGGGGGAAACGCCGccc                                                                                                                            >  2:4061034/1‑150 (MQ=255)
                                 aTAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTTGTTTTTTTTGGTGCATCTGCTCGGGGAAACCGCCGCCCTg                                                                                                                          >  2:912164/1‑149 (MQ=255)
                                                    cgcTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGAc                                                                                                    <  2:2518856/151‑1 (MQ=255)
                                                     gcTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTTTTCTGCATCTGCTCGGGGAAACCGCCGCCCCGAATCATAAAACCAGGaa                                                                                                        >  2:1804032/1‑146 (MQ=255)
                                                     gcTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGAcg                                                                                                   <  2:3363309/151‑1 (MQ=255)
                                                        gTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGAcgc                                                                                                  <  1:3330944/149‑1 (MQ=255)
                                                           tttGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTTGTTTTTTTTTCTTCCACTGCTCCGTGGAAACGCCGCCCCGTAAAATAAAAAAAAGAATGACGCgggg                                                                                             >  2:2338561/1‑148 (MQ=255)
                                                                  gCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGGTTTTTTTTCTTCATCTTCTCGGTGGAACCGCCGCCCCGAATCATAAAAACCGGGAAGACGCgggggaa                                                                                          >  1:625173/1‑141 (MQ=255)
                                                                     gTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGGTTTTTTTGCTGCATCTTCTCGGGGGAACCCCCGCCCCGGATCATAAAACCCGGGAAGACGCGGGGAAAGGGAggg                                                                                   >  2:3926618/1‑148 (MQ=255)
                                                                          tGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTTGTTTTTTTTTGTGGATCTGCTCGGGGAAAACCCCGCCCCGAATCAAAAAACCCGGAAAGACGCCGGGAAAAGGAGTgttgttt                                                                             >  1:205914/1‑150 (MQ=255)
                                                                                 gcgATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAAc                                                                      <  2:2845013/151‑1 (MQ=255)
                                                                                    aTGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGGTTTTTTTTCTGCATCTTCTCGGTGGAAACGCCCCCCTGAATCATAAAACCAGGAAAGACGCGGGGGAAAGGAGGGGTTTTTTTAAAAAc                                                                      >  2:1984923/1‑149 (MQ=255)
                                                                                     tGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAAcc                                                                     <  2:1372281/148‑1 (MQ=255)
                                                                                     tGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGGTTTTTTTTCTCTAAATGCTCGGGGGAAACGCCGCCCCGCATTCTTAAAGCAGGGATGACGGCGTGAAAGGAAGTGGTGTTTTtaaaa                                                                       >  1:1372011/1‑146 (MQ=255)
                                                                                            aCGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGTTTTTTTTGCTGCCTCTGCTCGGTGAAAACGCCGCCCTGAATCCTAAAACCAGGGATGAAGCGGGGGAAGGGGGGGGTGTTATTAAAACCGCTGtt                                                               >  2:508669/1‑147 (MQ=255)
                                                                                                gtgtTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTGGGTTTTTCTGCTGCATCTGCTCGGTGAAAACCCCGCCCCGAATCATAAAACCAGGAATGACGCGGGGGAAGGGAGTGGTGTTTTAAAAACCGCTGTTCACATaa                                                        >  1:2000500/1‑150 (MQ=255)
                                                                                                  gtTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATc                                                      <  2:3062356/150‑1 (MQ=255)
                                                                                                  gtTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCg                                                     <  1:970298/151‑1 (MQ=255)
                                                                                                   ttGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTTGTTTTTTTTGATGCTTCTTCTCGGGGAAAACACCGCCCTTAATCATAAAAACCGGAAAGACGCGGTGGAAAGAAGGGTTGTTATAAAAACCCCCGTTTCCATaa                                                        >  2:1586414/1‑147 (MQ=255)
                                                                                                            ccGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGtttt                                           <  1:2331048/151‑1 (MQ=255)
                                                                                                             cGTTATCGGCTTCATTTGTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGTCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAACCCGCTGTTCACACAACCGTCAAAGtttta                                          <  1:1323405/151‑2 (MQ=255)
                                                                                                                tATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTTTTCTTCATATGCTCTGGGGAAACCCCGACCCGAATCACAAAAACAGGGAAGACGCGGGGAAAGGTAGTGTTGTTATAAAAAACCCTGTTCACATAAACGACaaaaattt                                           >  1:3447936/1‑142 (MQ=255)
                                                                                                                tATCGGCTTCATTTTTGATTGGCGGGGTTGGTTTTTTTTGTTGCATCTGCTCGGTGAAAACGCCGCCCTGCATCATAAAAACAGGGATGACGCGGGGGAAGGGAGTGTTGTTTTAAAAAACCCTGTTTAAAAAAACGACaaaaattttca                                        >  1:3053106/1‑142 (MQ=255)
                                                                                                                 aTCGGCTTCATTTTTGATTGGCGGGGTTTGGTTTTTTTGCTTCATCTTCTCCGGGGAAACGCCCCCCTGAATCATAAAACCAGGGATGAAGCGGGGGAAAGGAGGGTTGTTATAAAAAACGCTGTTTACATAATTGACAAAAATTTTcac                                       >  1:3990526/1‑150 (MQ=255)
                                                                                                                      cTTCATTTTTGATTGGCGGGGTTGGTTTTTTTTGCTCCATCTTCTGCGGGAAAACGCCGCCCCGAAACAAAAAAGCAGGAATGAAGCGGGGTAAGGGAGGGGTGTTATAAAAACCGCAGTTCAAAAAATTTACCAAGTTTTGcaccaaca                                  >  2:1078579/1‑145 (MQ=255)
                                                                                                                          aTTTTTGATTGGCGGGTTTTGATTTTTCTGCTGCGTCTGTCCGGGTAAAAAACCGCCCTCTATTATTAAATCACGGACGGAGCGGGGAAAAGGAGGGGTTCTATAAAAAAAGACGTTCACACAAACGACAAAATTTTGTACAGacaaccg                              >  1:603792/1‑146 (MQ=255)
                                                                                                                            ttttGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGgcg                            <  2:2683465/150‑1 (MQ=255)
                                                                                                                              ttGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCt                          <  2:2784758/150‑1 (MQ=255)
                                                                                                                                   tGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTg                     <  2:2197623/150‑1 (MQ=255)
                                                                                                                                    ggCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGtt                   <  1:1984695/151‑1 (MQ=255)
                                                                                                                                         ggtttggtttTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATc               <  1:3427618/151‑1 (MQ=255)
                                                                                                                                                ttttttCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTcc        >  1:1602856/1‑151 (MQ=255)
                                                                                                                                                      ctgctgCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCg  >  2:430264/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                         
CCATGGTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGTTTGGTTTTTTCTGCTGCATCTGCTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGTGTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCG  >  NC_000913/3491879‑3492178

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: