Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA4,296,0900CT21.8% 86.1 / 24.4 56intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (24/19);  new base T (6/6);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913_3_vch_tpiA/4295942‑4296236
                                                                                                                                                     |                                                                                                                                                  
gTCTGATAAGGCGTGTATGCCGTATCCGACGTCAATGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAc                                                                                                                                                    <  2:158494/149‑1 (MQ=11)
    gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGAGACGCTTAACGCGTCTTCCCGGTCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                                >  1:112432/1‑149 (MQ=11)
         ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCTTGATTTTGTAGGTCGTATACGGCGTTTACGCCGCATCCGACGTCAACGCCTGATGc                                                                                                                                           >  2:125569/1‑149 (MQ=255)
                   ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCgg                                                                                                                                         >  2:511748/1‑140 (MQ=17)
                   ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                         <  1:511748/141‑1 (MQ=32)
                    cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcg                                                                                                                                >  2:149902/1‑149 (MQ=35)
                     gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGCATAAGACGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                               >  2:478077/1‑149 (MQ=25)
                          cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGAAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGCTGCGACGCTTGCCGCGTCtt                                                                                                                          >  1:383481/1‑149 (MQ=255)
                          cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                          <  1:149902/149‑1 (MQ=39)
                            aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCAGCGTCTTAt                                                                                                                        >  2:552671/1‑149 (MQ=25)
                                cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCgg                                                                                                                    >  2:125674/1‑149 (MQ=21)
                                 aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGCCCTACGCCAGACAGCGCAATCGCCTGATCTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                   <  1:697900/149‑1 (MQ=255)
                                  aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTGAGCCATCAGGTCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                  >  1:538747/1‑149 (MQ=11)
                                   cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                 >  2:309970/1‑149 (MQ=25)
                                    gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                <  1:626196/149‑1 (MQ=21)
                                      cTGATTCGACGCTTAACTCGTCGTTTTAGGCCTAGGCCAGACAGCGCAATAGCCTGATTTA‑GCGGGATTGTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACGTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                              >  2:693861/1‑147 (MQ=11)
                                      cTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                 <  1:183498/114‑1 (MQ=12)
                                      cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTGTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATGAATGCCTGCTGCGACGCTGGCCGCGTCTTCTCAGGCCtgtg                                                                                                              >  2:357896/1‑146 (MQ=11)
                                          tGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAg                                                                                                          >  2:603442/1‑149 (MQ=37)
                                                gcttgccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                                                                                                    <  1:491259/143‑1 (MQ=21)
                                                gCTTAACGCGTCCTATCCGGCCTCCGCCATACAGCGCTCTCGAC‑AATAAACTCTGGATTTTGTAGGTCGGATAAGGCTTGTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                                                                                                    <  1:572304/149‑1 (MQ=2)
                                                          tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                               >  2:216076/1‑112 (MQ=39)
                                                          tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                               <  1:216076/112‑1 (MQ=39)
                                                          tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTg                                                                                          >  2:626196/1‑149 (MQ=37)
                                                               tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                     >  1:727695/1‑149 (MQ=38)
                                                               tCAGGCCTACGCCAGACAGAGAAATAGCCTGATTCA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                     <  2:524356/149‑1 (MQ=255)
                                                                cAGGGCTAAGCCGGGCAGCGAAATTGTGGGATTTA‑GCGTGATTTTGTAGGTCGGGTAAGGGGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAg                                                                                    <  2:37894/149‑1 (MQ=11)
                                                                 aGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGc                                                                                   >  1:50095/1‑149 (MQ=38)
                                                                 aGGCCTACGCCAGACAGAGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                   >  1:324051/1‑149 (MQ=255)
                                                                         gCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                             >  1:429132/1‑83 (MQ=25)
                                                                         gCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                             <  2:429132/83‑1 (MQ=25)
                                                                            aGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                        <  1:478077/149‑1 (MQ=38)
                                                                              aCAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                 <  2:55269/74‑1 (MQ=25)
                                                                                aGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACGTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                    <  1:770129/149‑1 (MQ=38)
                                                                                                 tA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGc                                                   <  1:489828/149‑1 (MQ=39)
                                                                                                       tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGGCATCAACGCCTGATGCGACGCTTGCCCCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAt                                              >  1:346816/1‑149 (MQ=255)
                                                                                                        gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATc                                             >  2:340263/1‑149 (MQ=39)
                                                                                                           tttGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGa                                          <  1:100302/149‑1 (MQ=38)
                                                                                                            ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAAATCGACGTTTGATGCGACGCTTGCCGCGTCTTATCAGGCCGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAc                                         <  1:610362/149‑1 (MQ=11)
                                                                                                             tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATg                                                                                                                                            <  1:553059/49‑1 (MQ=255)
                                                                                                             tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCATTAGCCTGATTTAGCTTGATTTTGTAGGCCGGATAAGGCGTTTAAGCCGCATCCGCCa                                        >  2:61729/1‑149 (MQ=11)
                                                                                                                 ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTCGCGTGATTTTCTCGATCGGATAAGTCGTTTATGCCCCATCCGACCTCaa                                    >  2:11555/1‑149 (MQ=14)
                                                                                                                  ggcGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAc                                   <  1:231976/147‑1 (MQ=38)
                                                                                                                   tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                >  2:513813/1‑71 (MQ=19)
                                                                                                                   tCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGc                                                                                                                                      <  2:631263/49‑1 (MQ=255)
                                                                                                                   tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                <  1:513813/71‑1 (MQ=255)
                                                                                                                     ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCACCGcc                                >  1:480908/1‑149 (MQ=255)
                                                                                                                       aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTTCTGCCGCATCCGACATCAACGCCTg                              >  2:619927/1‑149 (MQ=255)
                                                                                                                         agggCGTTTCCGCCCCGCCAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAt                            <  2:633964/147‑1 (MQ=11)
                                                                                                                           ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAATGCGTTTATGCCGCATCCGACATCTAGCCCtgctgc                          >  2:304142/1‑149 (MQ=14)
                                                                                                                            gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCTAGCCAGCGCATTAGCCTGCTTTAGCGTGACTTTGTAGGTCGGATAAGGCGTTTATACCGCAAGCGGCGTCAACGCCGCATGCg                         >  2:151708/1‑149 (MQ=14)
                                                                                                                                    gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGACACCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTATGCCGCATCCGACATCGACGCCTGCTGCGGAGCTTac                 >  1:165697/1‑148 (MQ=11)
                                                                                                                                     ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                   <  1:181232/82‑1 (MQ=40)
                                                                                                                                     ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                   >  2:181232/1‑82 (MQ=40)
                                                                                                                                          tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                       <  1:125588/41‑1 (MQ=40)
                                                                                                                                          tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                       >  2:125588/1‑41 (MQ=40)
                                                                                                                                                   aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc  <  1:552975/149‑1 (MQ=32)
                                                                                                                                                     |                                                                                                                                                  
GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTA‑GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913_3_vch_tpiA/4295942‑4296236

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