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.

GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913_3_vch_tpiA/4295942‑4296236
                                                                                                                                                     |                                                                                                                                                  
gTCTGATAAGGCGTGTATGCCGTATCCGACGTCAATGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAc                                                                                                                                                    <  2:158494/149‑1 (MQ=11)
    gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGAGACGCTTAACGCGTCTTCCCGGTCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                                >  1:112432/1‑149 (MQ=11)
         ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGTATACGGCGTTTACGCCGCATCCGACGTCAACGCCTGATGc                                                                                                                                           >  2:125569/1‑149 (MQ=255)
                   ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCgg                                                                                                                                         >  2:511748/1‑140 (MQ=17)
                   ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                         <  1:511748/141‑1 (MQ=32)
                    cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcg                                                                                                                                >  2:149902/1‑149 (MQ=35)
                     gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGCATAAGACGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                               >  2:478077/1‑149 (MQ=25)
                          cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGCTGCGACGCTTGCCGCGTCtt                                                                                                                          >  1:383481/1‑149 (MQ=255)
                          cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                          <  1:149902/149‑1 (MQ=39)
                            aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCAGCGTCTTAt                                                                                                                        >  2:552671/1‑149 (MQ=25)
                                cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCgg                                                                                                                    >  2:125674/1‑149 (MQ=21)
                                 aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGCCCTACGCCAGACAGCGCAATCGCCTGATCTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                   <  1:697900/149‑1 (MQ=255)
                                  aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTGAGCCATCAGGTCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                  >  1:538747/1‑149 (MQ=11)
                                   cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                 >  2:309970/1‑149 (MQ=25)
                                    gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                <  1:626196/149‑1 (MQ=21)
                                      cTGATTCGACGCTTAACTCGTCGTTTTAGGCCTAGGCCAGACAGCGCAATAGCCTGATTTAGCGGGATTGTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACGTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                              >  2:693861/1‑147 (MQ=11)
                                      cTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                 <  1:183498/114‑1 (MQ=12)
                                      cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATGAATGCCTGCTGCGACGCTGGCCGCGTCTTCTCAGGCCtgtg                                                                                                              >  2:357896/1‑146 (MQ=11)
                                          tGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAg                                                                                                          >  2:603442/1‑149 (MQ=37)
                                                gcttgccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                                                                                                    <  1:491259/143‑1 (MQ=21)
                                                gCTTAACGCGTCCTATCCGGCCTCCGCCATACAGCGCTCTCGACAATAAACTCTGGATTTTGTAGGTCGGATAAGGCTTGTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                                                                                                    <  1:572304/149‑1 (MQ=2)
                                                          tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                               >  2:216076/1‑112 (MQ=39)
                                                          tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                               <  1:216076/112‑1 (MQ=39)
                                                          tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTg                                                                                          >  2:626196/1‑149 (MQ=37)
                                                               tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                     >  1:727695/1‑149 (MQ=38)
                                                               tCAGGCCTACGCCAGACAGAGAAATAGCCTGATTCAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                     <  2:524356/149‑1 (MQ=255)
                                                                cAGGGCTAAGCCGGGCAGCGAAATTGTGGGATTTAGCGTGATTTTGTAGGTCGGGTAAGGGGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAg                                                                                    <  2:37894/149‑1 (MQ=11)
                                                                 aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGc                                                                                   >  1:50095/1‑149 (MQ=38)
                                                                 aGGCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                   >  1:324051/1‑149 (MQ=255)
                                                                         gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                             >  1:429132/1‑83 (MQ=25)
                                                                         gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                             <  2:429132/83‑1 (MQ=25)
                                                                            aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                        <  1:478077/149‑1 (MQ=38)
                                                                              aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                 <  2:55269/74‑1 (MQ=25)
                                                                                aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACGTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                    <  1:770129/149‑1 (MQ=38)
                                                                                                 tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGc                                                   <  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)
                                                                                                                                                     |                                                                                                                                                  
GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  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: