Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,492,0280CT39.5% 54.9 / 15.3 38Q90Q (CAG→CAA) ppiApeptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base C (4/19);  new base T (15/0);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 5.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTCC  >  NC_000913/3491883‑3492172
                                                                                                                                                  |                                                                                                                                                  
ggTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTctgctg                                                                                                                                               <  1:369888/151‑1 (MQ=255)
   cAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCa                                                                                                                                             <  2:555400/150‑1 (MQ=255)
     aggaaggCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCt                                                                                                                                          <  1:484631/151‑1 (MQ=255)
            cGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGGGTTTTTTTCTGCctcttccccg                                                                                                                                    >  2:653703/1‑144 (MQ=255)
                aTCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGaa                                                                                                                               <  1:347713/150‑1 (MQ=255)
                  cGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGaa                                                                                                                               <  1:229586/148‑1 (MQ=255)
                    gCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAAccgc                                                                                                                          <  1:100171/151‑1 (MQ=255)
                    gCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGCACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGGTTTTTTTTCTGCCTCTT‑CTCGGGGGAAcc                                                                                                                            >  1:626385/1‑149 (MQ=255)
                               aaaGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAAt                                                                                                                <  1:307493/150‑1 (MQ=255)
                                  gAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGGTTTTTTTTCCGCATCTG‑TTTGGGGGAAACCCCGccc                                                                                                                     >  1:662256/1‑142 (MQ=255)
                                   aaCTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTTTGCTGCGTCTG‑CTCGGTGGAACCGCCGCCCTGaa                                                                                                                 >  1:22001/1‑145 (MQ=255)
                                   aaCTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATaa                                                                                                           <  1:550688/151‑1 (MQ=255)
                                     ctggctggTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTTGTTTTTTTTTCTTCATCTG‑CTCCGGGACACCGCCCCCCTCCATCATaaa                                                                                                          >  2:601341/1‑150 (MQ=255)
                                            gtggCGCTGTCTGTGGCAGCGGGCCGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCTCCGCCCTGAATCATAAAGCCATGaa                                                                                                  <  2:548496/151‑1 (MQ=255)
                                                           tCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGCCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTaaaa                                                                                     <  2:161277/149‑1 (MQ=255)
                                                             aGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTTGTTTTTTCTTCTGCCTCTT‑CTCGGTGGAAACGCCGCCCCGAAACATAAAACCAGGGATGAcggcggggaa                                                                                     >  2:165104/1‑139 (MQ=255)
                                                             aGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGGTTTTTTTTCTGCATCTG‑CTCGGGGAAACCGCCGCCCTGCATCATAAAGACAGGGATAACGCggggga                                                                                      >  2:256038/1‑142 (MQ=255)
                                                                        ccATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTTGTTTTTTTTGCTGCATCTC‑CTCCGTGGAAACCCCGCCCCGTATCCTAAAACCCGGGATGACCCGGTGAAAGGTAGA‑GGTGGtat                                                                      >  1:90928/1‑150 (MQ=255)
                                                                             gcgATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAAc                                                                <  2:624388/151‑1 (MQ=255)
                                                                                    tGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTTGGTTTTTTTGCTGCATCTG‑CTCGGGGAAACCGCCCGCCCGAATCATAAAAACCGGAATGACGCGGGGAAAAGGAGG‑GTTGTTATTAAAAACGCTGtt                                                         >  1:256360/1‑151 (MQ=255)
                                                                                    tGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGTTTTTTTTTCTGCATCCG‑CTCGGTGGAACCGCCGCCCCGAATCATAAAACCCGGAATGACGCGGGGAAAGGGAGT‑GGTTTTTTAAAAAACCCTggt                                                         >  1:129896/1‑149 (MQ=255)
                                                                                    tGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTCGC‑GGGTTTGGGGTTTTTTCCTACAACTGCCTCGGGCAAAGCG‑CGCCCCCTATAATAAAGGCCGGGCTAACGCGGCGGGAGGGAGA‑GGTTTTATATAAAcc                                                               >  2:645708/1‑145 (MQ=255)
                                                                                       cACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGGTTTTTTTTCTTCCTCTG‑CTCCGGGGAACCGCCCCCCTGAATCCTAAAACCAGGAAAGAAGCGGGGGAAGGGAGG‑GTTGTTTTTAAAACCGCCGtttcc                                                      >  1:270945/1‑148 (MQ=255)
                                                                                            gtgtTGCGCAGGCCGTTATCGGCTTCATGTTTGATGGGC‑GGGTTTGGTGTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTcaca                                                     <  2:647142/147‑1 (MQ=255)
                                                                                                tGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGGTTTTTTTGCTGCATCTG‑CTCGGTGGAACCGCCGCCCTGTATCATAAAAACAGGGATGACGCGGGGAAAGGGAGG‑GTTGTTATAAAAAACGCCGTTtaca                                                     >  2:694095/1‑143 (MQ=255)
                                                                                                 gcgcAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGAc                                             <  1:761642/150‑1 (MQ=255)
                                                                                                  cgcAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTTTGCTGCATCTG‑CTCGGGGAAACCGCCGCCCTGAATCATAAAACCAGGAATGACGCGGGGAAAAGTAGT‑GGTGTTATAAAAAACGCTGTTCACATAATCGaaa                                            >  1:559507/1‑148 (MQ=255)
                                                                                                  cgcAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACaa                                           <  2:514259/151‑1 (MQ=255)
                                                                                                          gTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGc                                   <  1:169523/151‑1 (MQ=255)
                                                                                                          gTTATCGGCTTCATTTTTGATTGGC‑GGGGTTTGTTTTTTTTGCTGCATCTT‑CTCTGTGGAAACGCCGCCCTGAATTATAAATCCCGGAATGACGCGGGGAAAAGGAGG‑GTTGTTTTAAAAAACGCTGTTCACATAATCGACAAAAtttt                                     >  2:100426/1‑149 (MQ=255)
                                                                                                               cGGCTTCATTTTTGATTGGC‑GGGGTTGGTTTTTTCTTCTGCCTCTG‑CTCCGTGAAAACGCCGCCCTGAATCATAAAAACCGGAATGGCGCGGGGGAAGGGAGT‑GTTGTTTTAAAAAACGCTGTTTACATAATCGACAAAGattt                                     >  2:506576/1‑144 (MQ=255)
                                                                                                                 gCTTCATTTTTGATTGGCGGGGTTTTGTTTTTTTTCCT‑CATCTT‑CTCTGTGAAACCGCCGCCCTGGATTATAAAAACCGAAAATAAGCGGGGGAACGGAGTGGTTGTTA‑AAAAACCCCTCGTCACATCAACGACAAAGtttttcccc                                >  2:614254/1‑142 (MQ=255)
                                                                                                                      aTTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGgc                       <  2:127930/151‑1 (MQ=255)
                                                                                                                                ggC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGt              <  2:426912/150‑1 (MQ=255)
                                                                                                                                  c‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGttt            >  2:276013/1‑150 (MQ=255)
                                                                                                                                  c‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTa           <  2:396156/151‑1 (MQ=255)
                                                                                                                                    gggTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTa           <  1:241764/150‑1 (MQ=255)
                                                                                                                                            ttttttCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTcc  >  2:605260/1‑151 (MQ=255)
                                                                                                                                             tttttCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTcc  <  1:576564/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGTCAAGGAAGGCGTTATCGGCAACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTG‑CTCGGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATGACGCGGTGAAAGGTAGT‑GTTGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTGCACAGACACTGGCGCTTTTTGTTTATCCAGCTCC  >  NC_000913/3491883‑3492172

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