Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0020941,004,3600TG42.1% 85.6 / 36.9 77G22G (GGT→GGGtrpAtryptophan synthase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (9/35);  new base G (32/0);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 2.30e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGTCTTGAAGCAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCAGACCC  >  CP002094/1004264‑1004458
                                                                                                |                                                                                                  
aaCGTCTTGAAGCAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGa                                                                                                  <  2:205507/99‑1 (MQ=255)
 aCGTCTTGAAGCAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGag                                                                                                 >  1:164687/1‑98 (MQ=255)
 aCGTCTTGAAGCAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGAc                                                                                                 >  2:1015926/1‑99 (MQ=255)
     cTTGAAGCAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCCATATTATGGCAGGTGACCATg                                                                                             <  2:980789/99‑1 (MQ=255)
         aaGCAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGaaaa                                                                                         >  1:1462013/1‑98 (MQ=255)
         aaGCAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCATGAAag                                                                                         >  2:1391790/1‑99 (MQ=255)
           gCAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGGCCATGAAaaag                                                                                       >  1:688540/1‑99 (MQ=255)
            cAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAgg                                                                                      <  1:1183547/99‑1 (MQ=255)
            cAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCATGAAAAAgg                                                                                      >  2:616486/1‑99 (MQ=255)
              gAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGAt                                                                                    <  1:74241/99‑1 (MQ=255)
              gAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCCTGAAAAAGGAt                                                                                    >  2:238392/1‑99 (MQ=255)
              gAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATATTGTTCCCTATATTATGGCAGGGGACCATGAAAGAGGAt                                                                                    >  2:540837/1‑99 (MQ=255)
                aaaaGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTa                                                                                  <  1:1593471/99‑1 (MQ=255)
                 aaaGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGGCCATGAAAGAGGATTAg                                                                                 >  1:760282/1‑99 (MQ=255)
                  aaGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGa                                                                                >  2:68566/1‑99 (MQ=255)
                   agagGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGaa                                                                               >  2:144830/1‑99 (MQ=255)
                   agagGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGGCCATGAAAGAGGATTAGaa                                                                               >  1:66558/1‑99 (MQ=255)
                     agGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAgg                                                                             <  2:1336965/99‑1 (MQ=255)
                        tGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCt                                                                          <  1:1017606/99‑1 (MQ=255)
                        tGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCCTGAAAAAGGATTAGAAAGGCt                                                                          >  2:1596312/1‑99 (MQ=255)
                        tGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCATGAAAGAGGATTAGAAggggt                                                                          >  2:1169997/1‑97 (MQ=255)
                          aaaTAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCATGAAAAAGGGTTAGAAgggggca                                                                        >  2:995337/1‑95 (MQ=255)
                               atGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGaa                                                                   <  1:106316/99‑1 (MQ=255)
                                tGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGaaa                                                                  <  2:700350/99‑1 (MQ=255)
                                tGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGGCCATGAAAAAGGATTAGAAGGGGGCAATGaaa                                                                  >  1:1530654/1‑99 (MQ=255)
                                  aaaaCTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCCTGAAAGAGGATTAGAAGGGGTCAATGAAAcc                                                                >  2:1528867/1‑99 (MQ=255)
                                  aaaaCTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCATGAAAAAGGATTAGAAGGGGTCAATGAAAcc                                                                >  2:138611/1‑99 (MQ=255)
                                    aaCTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCAt                                                              <  1:1304452/99‑1 (MQ=255)
                                    aaCTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCAt                                                              <  2:316649/99‑1 (MQ=255)
                                     aCTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGGTCAATGAAAACAtt                                                             >  2:1086827/1‑99 (MQ=255)
                                          aCAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAAct                                                        <  2:129046/99‑1 (MQ=255)
                                           cAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCATTAAAAAGGATTAGAAGGGGTCCATGAAACAATTCAActt                                                       >  2:778186/1‑99 (MQ=255)
                                           cAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCATGAAAGAGGATTAGAAGGGGTCAATGAAAACATTAAAatt                                                       >  2:139468/1‑99 (MQ=255)
                                           cAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGAACATGAAAAAGGATTAGAAAGGCTCAATGAAACCATTCAAatt                                                       >  2:803074/1‑99 (MQ=255)
                                                  aaTTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGGCCCCTAAAAACGGATATAAAGGCCCCAAGGAACCAATCAACTTCTTGTAc                                                >  2:350454/1‑99 (MQ=255)
                                                  aaTTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAAATTCTTGAAc                                                >  2:143489/1‑99 (MQ=255)
                                                    ttGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGCCCCGCCAAAAGGGTTTAAAAGGCCCAATGACAACATTCAAATTTTTTAaaaa                                              >  2:1284335/1‑99 (MQ=255)
                                                    ttGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGGACATGAAAAAGGATTAGAAGGGGTACATGAAAACAATCAACTTTTTTAACaa                                              >  1:1429417/1‑99 (MQ=255)
                                                     tGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGACCATGAAAGAGGAGTAGAAAGGGTACATGAAAACATTCAACCTCTTGAACAAg                                             >  2:1090632/1‑99 (MQ=255)
                                                      gAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGc                                            <  2:13221/99‑1 (MQ=255)
                                                       aGCAATAAAAAAAATAATTTTTTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGcc                                           <  2:448764/99‑1 (MQ=255)
                                                         cAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGGGAAGATGGAAGAGACATAGAAGGGCGCAATAAAACCATTCAAATTCTTTAACAAAcccg                                         >  2:1388807/1‑97 (MQ=255)
                                                           aTAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGa                                       <  2:327843/99‑1 (MQ=255)
                                                            tAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAg                                      <  1:680156/99‑1 (MQ=255)
                                                              aaaaaaaTAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGtt                                    <  1:872789/99‑1 (MQ=255)
                                                              aaaaaaaTAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGtt                                    <  1:582207/99‑1 (MQ=255)
                                                               aaaaaaTAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCTTTCAACTTCTTGAACAAGCCGGGAGTTc                                   <  1:479533/99‑1 (MQ=255)
                                                               aaaaaaTAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTc                                   <  1:69834/99‑1 (MQ=255)
                                                                aaaaaTAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCt                                  <  2:267213/99‑1 (MQ=255)
                                                                  aaaTAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGc                                <  1:312313/99‑1 (MQ=255)
                                                                       atTTTGTTCCCTATATTATGGCAGGGGGCCCTGAAAAAGGATTAGAAGGGGTCAAAGAAAACATTCAACTTCTTTTAAAAACCCGGGGGTCTGCCCTTg                           >  1:1340707/1‑99 (MQ=255)
                                                                       atTTTGTTCCCTATATTATGGCAGGGGACCCTGAAAGAGGGGTTGAAGGGCTCAATGAAACCATTCAACTTCTTTAACAAGCCGGGAGGTCTGCCATTg                           >  2:641713/1‑99 (MQ=255)
                                                                            gTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAACGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAAtt                      <  2:528867/99‑1 (MQ=255)
                                                                                ccTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGgtg                  <  2:1405918/99‑1 (MQ=255)
                                                                                ccTATATTATGGCAGGGGACCATGAAAAAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGgtg                  >  2:177500/1‑99 (MQ=255)
                                                                                 cTATATTATGGCAGGGGACCATGAAAGAGGATTAGAAGGGGTCAATGAAACCATTCAAATTCTTTAACAAGCCGGGAGTTCTGCCATTGAAATTGgtgg                 >  2:302523/1‑98 (MQ=255)
                                                                                  tataTTATGGCAGGGGACCATGAAAGAGGATTAGAAGGGGTCAATGAAAACATTCAAATTCTTGAACAAGCCGGGAGGTCTGCCCTTGAAATTGGTGTa                >  2:166204/1‑99 (MQ=255)
                                                                                   atatTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAACCCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTAc               <  2:533344/99‑1 (MQ=255)
                                                                                   atatTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTAc               <  2:397732/99‑1 (MQ=255)
                                                                                   atatTATGGCAGGGGACCATGAAAGAGGATTAGAAGGGGTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTAc               >  2:578722/1‑99 (MQ=255)
                                                                                    tattatGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTAcc              <  2:2630/99‑1 (MQ=255)
                                                                                     attatGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTAccc             <  1:16729/99‑1 (MQ=255)
                                                                                     attatGGCAGGGGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTAccc             >  2:1324317/1‑99 (MQ=255)
                                                                                        atGGCAGGGGACCCTGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGGTCTGCCATTGAAATTGGGGTACCCTtc          >  2:76098/1‑99 (MQ=255)
                                                                                         tGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTtct         <  2:993709/99‑1 (MQ=255)
                                                                                         tGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTtct         <  1:1582171/99‑1 (MQ=255)
                                                                                         tGGCAGGGGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGGTCTGCCATTGAAATTGGGGTACCCTtct         >  1:42963/1‑99 (MQ=255)
                                                                                          ggCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTtctc        <  2:835832/99‑1 (MQ=255)
                                                                                           gCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCa       >  1:676308/1‑99 (MQ=255)
                                                                                           gCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCa       >  1:989470/1‑99 (MQ=255)
                                                                                           gCAGGAGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCa       >  2:1465759/1‑99 (MQ=255)
                                                                                             aGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCAGa     <  1:914741/99‑1 (MQ=255)
                                                                                             aGGTGACCAGAAAAGAGTATTAGAAGGTCTCAATGAAACCTTTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCAGa     <  1:956162/99‑1 (MQ=255)
                                                                                              ggTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCAGAc    >  1:399649/1‑99 (MQ=255)
                                                                                              ggTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCAGAc    >  2:425745/1‑99 (MQ=255)
                                                                                               gTGACCATGAATGAGGATTACAGGGGTTGAATGAAACAATTAGCCTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCATTCTCAGGcc   <  2:718135/99‑1 (MQ=255)
                                                                                               gTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCAGAcc   <  2:464178/99‑1 (MQ=255)
                                                                                                tGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCAGAccc  <  2:22031/99‑1 (MQ=255)
                                                                                                |                                                                                                  
AACGTCTTGAAGCAGAAAAAGAGGTGAAATAATGAAAACTTTACAAGCAAAATTGAGCAATAAAAAAAATAATTTTGTTCCCTATATTATGGCAGGTGACCATGAAAGAGGATTAGAAGGGCTCAATGAAACCATTCAACTTCTTGAACAAGCCGGGAGTTCTGCCATTGAAATTGGTGTACCCTTCTCAGACCC  >  CP002094/1004264‑1004458

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