Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913354,2820TG51.8% 67.9 / 84.7 54V526G (GTT→GGT) prpEpropionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (5/21);  new base G (28/0);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 1.26e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-10

TGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGC  >  NC_000913/354140‑354430
                                                                                                                                              |                                                                                                                                                     
tGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTgggg                                                                                                                                                  >  1:213843/1‑148 (MQ=255)
 ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGa                                                                                                                                               <  1:27835/150‑1 (MQ=255)
  gggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGaaa                                                                                                                                             <  1:372733/151‑1 (MQ=255)
   ggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTggggggga                                                                                                                                              >  1:704479/1‑145 (MQ=255)
     cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt                                                                                                                                          <  2:356407/151‑1 (MQ=255)
      aTCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATg                                                                                                                                         <  1:743190/151‑1 (MQ=255)
      aTCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTGGGGGGGAAGGTg                                                                                                                                         >  1:248965/1‑151 (MQ=255)
           tGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTggggggttaagtggggtg                                                                                                                                    >  2:549341/1‑137 (MQ=255)
                 tGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGGAAggg                                                                                                                             >  1:203961/1‑151 (MQ=255)
                       ttATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTGGGGGGGAAAATGTGCTGAAAGGGGAggtg                                                                                                                       >  2:516146/1‑151 (MQ=255)
                        tATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAggtgg                                                                                                                      >  2:752641/1‑151 (MQ=255)
                         aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtgg                                                                                                                      <  1:44322/150‑1 (MQ=255)
                             cTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCCgggggcgggg                                                                                                                 >  2:90534/1‑149 (MQ=255)
                              tttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCCGGAAggggggggggggg                                                                                                                   >  1:558168/1‑138 (MQ=255)
                                tATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGGAgggggcggggggg                                                                                                              >  2:931179/1‑149 (MQ=255)
                                      cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGGAAGGGGGGGGGGGGGGGGCGgtt                                                                                                           >  1:866622/1‑148 (MQ=255)
                                       gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAAAGGGCTGGAAGGGCAggggggggggggggtt                                                                                                           >  1:22024/1‑142 (MQ=255)
                                           gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTcc                                                                                                   <  1:219208/151‑1 (MQ=255)
                                              cTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAAGGGAGGTGGCGGGGGCGGTTGTCATTTCGa                                                                                                 >  2:314567/1‑150 (MQ=255)
                                                 aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaa                                                                                                <  1:431921/148‑1 (MQ=255)
                                                             aCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCt                                                                                 <  2:31351/151‑1 (MQ=255)
                                                             aCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGGGGCGGGGGGGTTTGTCCTTTCGAAAAAGAGCGGAAGTCt                                                                                 >  1:965381/1‑151 (MQ=255)
                                                               gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTg                                                                                <  2:354823/150‑1 (MQ=255)
                                                                ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTGGGGGGGGAAGAAGCGCTGGAAAGGGCGGGGGGGGGGGGGATTGGCATTTCGGAagagcgggagcggaccgg                                                                              >  1:234858/1‑137 (MQ=255)
                                                                 tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGa                                                                              <  2:409943/150‑1 (MQ=255)
                                                                  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGaa                                                                             <  1:610810/150‑1 (MQ=255)
                                                                  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg                                                                            <  2:808001/151‑1 (MQ=255)
                                                                  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAAAAGAGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGGGCGACAGTCTgg                                                                               >  2:818730/1‑148 (MQ=255)
                                                                        cATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGTGGCGGGGGGGGTTTTCATTCCGAAAGAGAGCGACAGTCTGGGAGATCggg                                                                      >  2:782351/1‑149 (MQ=255)
                                                                         aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGGAGGGGGGGGGGGGGGTTGTTATTTCGGAAGAGAGCGACCGTCTGGAAGATCg                                                                        >  2:70003/1‑148 (MQ=255)
                                                                         aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGCTGGAAAGGGAGGGGGGGGGGGGGGTTGTTATTCCGAAAGAGAGGGACAGGCTGGGAGGTc                                                                         >  1:855569/1‑147 (MQ=255)
                                                                             gCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAAGCGCTGAAAGGGGCGGGGGGGGGGGGGGTTTTCATTTCCAAAAAGAGCGACAGGCTGGAAGGTCGTGAtg                                                                   >  1:554865/1‑149 (MQ=255)
                                                                                    aCGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGGTTGTCATTCCGAAAGAGAGCGACAGGCTGGGAGATCGTGATGTGGCGCAc                                                           >  2:260692/1‑150 (MQ=255)
                                                                                     cgcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTc                                                         <  2:279023/151‑1 (MQ=255)
                                                                                         cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCCGAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGGGCACTCGc                                                       >  1:602529/1‑149 (MQ=255)
                                                                                         cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCCGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGGGCACTccc                                                       >  1:389955/1‑147 (MQ=255)
                                                                                         cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTTGAAGATCGTGATGTGGCGCACTc                                                         >  2:688523/1‑147 (MQ=255)
                                                                                          gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGGAGGTGGCGGGGGCGTTTGGCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATTTGGCGCACTCGc                                                       >  2:119816/1‑148 (MQ=255)
                                                                                           agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAaga                                                   <  2:902200/151‑1 (MQ=255)
                                                                                           agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTTCGAAAGAGAGCGACAGGCTGGGAGATCGTGATGTGGCGCCCTc                                                         >  1:160050/1‑145 (MQ=255)
                                                                                             aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCCGGAAGGGCAGGGGGCGGGGGGGTTTGTCATTCCGGAAGAGAGCGACAGTCTGGAAGATCGGGGTGTGGCGCACTCGcca                                                     >  1:961356/1‑145 (MQ=255)
                                                                                                         aTCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATg                                      <  2:736242/150‑1 (MQ=255)
                                                                                                           ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAGATGCGCTGGAAAGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGGAAATCGTGGTGTGGCGCACTCGCAAGAGAAGGGGATTATGgcg                                   >  2:336972/1‑151 (MQ=255)
                                                                                                              cAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtg                                <  2:380212/151‑1 (MQ=255)
                                                                                                              cAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtg                                <  2:334444/151‑1 (MQ=255)
                                                                                                                 ccATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGGTTGTCCTTCCGGAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGCGCCCTCGCAAGAGAAGGCGGTTATGGCGCtgggg                             >  1:495411/1‑149 (MQ=255)
                                                                                                                  cATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGTGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtgg                            >  1:769472/1‑151 (MQ=255)
                                                                                                                     ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGAc                          <  2:129892/150‑1 (MQ=255)
                                                                                                                      cGGGCGTTGCCGAAGTGGCGGAGGGTGGGGTGGAAAATGCGCTGAAAAGGGAGGGGGCGGTGGGGGTTGTCGTTTCGAAAAAAAACGACAGTCTGGGAGAACGGGGTGTGGGGGACTCGCGAGAGAAAGCGAATATGGCGCTGGGGGAc                          >  1:96040/1‑149 (MQ=255)
                                                                                                                                  aaGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCa             >  1:137012/1‑150 (MQ=255)
                                                                                                                                  aaGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCa             >  2:791016/1‑150 (MQ=255)
                                                                                                                                       gCGGTGGTTGGGGTGAAAGATGCTCTCAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGCGGCCCACTCGCATGAGAAGGCGATTATGGCCCTGGTGGACAGCCTGATCGGCAACtttc       <  1:549333/151‑2 (MQ=255)
                                                                                                                                       gCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTg       <  2:635342/151‑1 (MQ=255)
                                                                                                                                        cGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTgg      <  1:90530/151‑1 (MQ=255)
                                                                                                                                            ggtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccgc  >  1:121810/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
TGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGC  >  NC_000913/354140‑354430

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