Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0129671,132,8510AG20.6% 526.8 / 30.4 240N65N (AAT→AACyceIYceI family protein
Reads supporting (aligned to +/- strand):  ref base A (50/139);  new base G (49/0);  total (101/139)
Fisher's exact test for biased strand distribution p-value = 2.20e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 268 total aligned reads displayed.

TTATCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGT  >  NC_012967/1132780‑1132921
                                                                       |                                                                      
ttATCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGa                                                                        <  1:4165072/72‑1 (MQ=255)
 tATCGCGTTCGGCGTGATTAGTATCGACGCNGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGAt                                                                       <  2:620589/72‑1 (MQ=255)
  aTCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGAtt                                                                      <  2:3106656/72‑1 (MQ=255)
    cgcgTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGtttt                                                                    >  1:4926679/1‑72 (MQ=255)
       gTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGAttttttt                                                                 <  2:9157659/72‑1 (MQ=255)
        ttCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTc                                                                <  1:7254358/72‑1 (MQ=255)
         tCGGCGTGAGTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGGCGGGctttttttg                                                               >  1:908256/1‑70 (MQ=255)
          cGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGt                                                              >  1:8288531/1‑72 (MQ=255)
          cGGCGTGATTAGTATCGACGCCGGTGGTGTTAATTGTAACATTCACTTTATCGGCAGGCGGGTGTTTTTTGt                                                              >  2:3499183/1‑72 (MQ=255)
            gCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGCTTTTTTTTGTct                                                            >  1:6739669/1‑71 (MQ=255)
              gTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCaaa                                                          <  2:2237187/72‑1 (MQ=255)
               tGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAg                                                         <  2:1339372/72‑1 (MQ=255)
                gATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTCGTCCAAgg                                                        >  2:7514700/1‑72 (MQ=255)
                gATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTCGTCAAAgg                                                        >  1:2468690/1‑72 (MQ=255)
                gATTAGTATCGACGCTGGCGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAgg                                                        <  1:3323168/72‑1 (MQ=255)
                 aTTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGt                                                       <  1:3256342/72‑1 (MQ=255)
                 aTTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGt                                                       <  2:4156392/72‑1 (MQ=255)
                  ttAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTa                                                      <  2:6681208/72‑1 (MQ=255)
                   tAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGATTTTTTCGTCAAAGGTaa                                                     >  2:7879869/1‑72 (MQ=255)
                   tAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTaa                                                     <  2:7840244/72‑1 (MQ=255)
                    aGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGCTTTTTTTGTTAAAAGTaaa                                                    >  2:3893954/1‑72 (MQ=255)
                    aGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTaaa                                                    <  2:6859044/72‑1 (MQ=255)
                     gTATCGACCCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGGTTTTTTCGTCTAAAGTaaaa                                                   >  2:131335/1‑72 (MQ=255)
                      tATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAg                                                  <  2:4026091/72‑1 (MQ=255)
                       aTCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGt                                                 <  1:3665395/72‑1 (MQ=255)
                        tCGACGCTGGTGTTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTa                                                <  2:7723559/72‑1 (MQ=255)
                          gACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTAcc                                              <  1:1390213/72‑1 (MQ=255)
                           aCGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCACAGGTGAAAAGACCg                                             >  2:2799445/1‑72 (MQ=255)
                           aCGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAAGGAAAAATAACg                                             >  2:2202710/1‑72 (MQ=255)
                           aCGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCCGGTTTTTTTTGTCAAAGGTAAAAGTaacc                                             >  1:4482570/1‑69 (MQ=255)
                           aCGCTGGTGGTGTTAATTGTCACATTCAATTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCg                                             <  1:2704290/72‑1 (MQ=255)
                            cGCTGGTGGTGTTAATTGTCACATTCACTTTATCTGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGt                                            <  2:8075932/72‑1 (MQ=255)
                            cGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGt                                            <  1:6083123/72‑1 (MQ=255)
                             gCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAGGTAAAAGTACCGTc                                           >  1:2832904/1‑72 (MQ=255)
                             gCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAGGTAAAAGTAACGTc                                           >  2:8896933/1‑72 (MQ=255)
                             gCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTCGTCAAAGGGAAAAATAACGTc                                           >  1:5491494/1‑72 (MQ=255)
                               tggtggTGTTAATTGTCACATTCACTTTATCGGCAGCCGGCCTTTTTTTTTCCAAGGTAAAAATAACGTcca                                         >  1:2390677/1‑70 (MQ=255)
                               tggtggTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGa                                         <  2:3551807/72‑1 (MQ=255)
                                ggtggtGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTCGTTAAAGGTAAAAATACCTTCgga                                        >  1:5997329/1‑70 (MQ=255)
                                 gtggtgTTAATTGTCACATTCACTTTATCGGCAGCCGGGGTTTTTTTGTCAAAGGGTAAAGTACCGTCGaaa                                       >  1:8693465/1‑72 (MQ=255)
                                  tggtgTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAAt                                      <  2:2732011/72‑1 (MQ=255)
                                     tgtTAATTGTCTCATTCACTTTATCGGCAGCCGGGCTTTTTTTGTCAAAGGTAAAAATACCGTCGAAATCtt                                   >  2:7346810/1‑72 (MQ=38)
                                     tgtTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCtt                                   <  2:1428360/72‑1 (MQ=255)
                                      gtTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAGGTAAAAGTAACGTCGAAATCttt                                  >  2:191097/1‑72 (MQ=255)
                                      gtTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCttt                                  <  2:2453837/72‑1 (MQ=255)
                                       ttAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTa                                 <  1:7707773/72‑1 (MQ=255)
                                       ttAATGGTCACAGTCACTGTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTa                                 <  1:5819691/72‑1 (MQ=255)
                                        tAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAGGGAAAAGTACCGTCGAAAACTTTaa                                >  2:5600393/1‑72 (MQ=255)
                                        tAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTaa                                <  2:7447031/72‑1 (MQ=255)
                                         aaTTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTaaa                               <  1:9650006/72‑1 (MQ=255)
                                          atagTCACATTCACTTTATGGGCAGCCGGATTTTTTTCGTCANAGGTAAAAGTACCGTCGAAATCTTTAAAg                              <  2:841526/69‑1 (MQ=255)
                                          aTTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAg                              <  1:9116982/72‑1 (MQ=255)
                                          aTTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAg                              <  1:6947773/72‑1 (MQ=255)
                                          aTTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAg                              <  1:81655/72‑1 (MQ=255)
                                            tGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTTTTAAAAGTAAAAATAACGTCGGAAACTTTTAAAGt                            >  2:3261738/1‑72 (MQ=255)
                                            tGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGt                            <  1:7106563/72‑1 (MQ=255)
                                            tGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGt                            <  2:7534236/72‑1 (MQ=255)
                                             gTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTg                           >  1:9905947/1‑72 (MQ=255)
                                               cacaTTCACTTTATCGGCAGTCGGGTTTTTTTTGTCAAAGGGAAAAAGACCGGCGGAATCTTTAAAGGTGcc                         >  1:8876667/1‑72 (MQ=255)
                                               cacaTTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGcc                         <  2:1137637/72‑1 (MQ=255)
                                                acaTTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCg                        <  1:1880864/72‑1 (MQ=255)
                                                 caTTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAAGTAAAAATAACGTCGAAAACTTTTAAAGTGCCGt                       >  2:3360290/1‑72 (MQ=255)
                                                  aTTCACTTTATCGGCAGCCGGGTTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGta                      >  2:9652867/1‑72 (MQ=255)
                                                    tCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAGGGAAAAGTACCGTCGAAAACTTTTAAGGTGCCGtata                    >  2:342821/1‑72 (MQ=255)
                                                    tCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGtata                    <  2:8836830/72‑1 (MQ=255)
                                                    tCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGtata                    <  2:8566249/72‑1 (MQ=255)
                                                     cACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATaa                   <  2:3794646/72‑1 (MQ=255)
                                                     cACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATaa                   <  2:4103833/72‑1 (MQ=255)
                                                     cACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGAATaa                   <  1:7490966/72‑1 (MQ=255)
                                                      aCTTTATCGGCAGCCGGGTTTTTGTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTTAAAGTGCCGTATAAc                  >  1:6324888/1‑72 (MQ=255)
                                                        tttATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCa                <  2:5445140/72‑1 (MQ=255)
                                                          tATCGGCAGCCGTATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGc              >  1:3524487/1‑72 (MQ=255)
                                                          tATCGGCAGCCGGGTTTTTTTCGTCNAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGc              >  1:1095702/1‑72 (MQ=255)
                                                          tATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGc              >  1:9264482/1‑72 (MQ=255)
                                                           aTCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAATTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCt             <  1:5423232/72‑1 (MQ=255)
                                                           aTCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCCGCt             >  1:2347225/1‑72 (MQ=255)
                                                           aTCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCt             <  1:9575639/72‑1 (MQ=255)
                                                           aTCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCt             <  1:2764687/72‑1 (MQ=255)
                                                           aTCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGGACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCt             >  2:9456604/1‑72 (MQ=255)
                                                            tCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCta            >  2:7314322/1‑72 (MQ=255)
                                                            tCGGCAGCCGGATTTTTTTAGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCta            <  1:318946/72‑1 (MQ=255)
                                                             cGGCAGCCGGGTTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGGGCCGTATAACCAGCtat           >  1:4069786/1‑72 (MQ=255)
                                                              ggCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCtata          >  1:856968/1‑72 (MQ=255)
                                                              ggCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCtata          >  1:4306356/1‑72 (MQ=255)
                                                              ggCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCtata          <  1:3852118/72‑1 (MQ=255)
                                                               gaggCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCGTTAAAGGTGCCGTATAACCAGCTATAg         <  1:1801728/69‑1 (MQ=255)
                                                                cAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATNTTTAAAGGTGCCGTATAACCAGCTATAGc        <  2:484919/72‑1 (MQ=255)
                                                                cAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGc        >  2:4779059/1‑72 (MQ=255)
                                                                cAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGc        >  2:101611/1‑72 (MQ=255)
                                                                cAGCCGGACTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGc        >  2:6425340/1‑72 (MQ=255)
                                                                 aGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGcc       <  2:4488341/72‑1 (MQ=255)
                                                                  gCCGGATTTTTTTCGTCAAAGGTAAAAGTNCCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCa      >  2:493037/1‑72 (MQ=255)
                                                                   ccGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCaa     >  1:4706354/1‑72 (MQ=255)
                                                                    cGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCCGCTATAGCCAAg    >  1:7915390/1‑72 (MQ=255)
                                                                    cGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAg    <  1:5321948/72‑1 (MQ=255)
                                                                    cGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAg    <  1:1308206/72‑1 (MQ=255)
                                                                    cGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAg    >  2:1629390/1‑72 (MQ=255)
                                                                    cGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAACCTTTAAAGGTGCCGTATAACACGCTATAGCCAAg    >  1:294910/1‑72 (MQ=255)
                                                                      gATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGt  <  2:3701115/72‑1 (MQ=255)
                                                                      gATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGt  <  2:6076838/72‑1 (MQ=255)
                                                                       |                                                                      
TTATCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGT  >  NC_012967/1132780‑1132921

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: