Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,231,4070GA76.7% 66.6 / 23.0 43E18K (GAA→AAA) prpCpossible PP2C protein phosphatase
Reads supporting (aligned to +/- strand):  ref base G (3/7);  new base A (17/16);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGG  >  CP000730/1231270‑1231542
                                                                                                                                         |                                                                                                                                       
tACGTAGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGGAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAg                                                                                                                                      >  2:22384/1‑141 (MQ=255)
                  cGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTAtttttt                                                                                                                    <  2:65821/141‑1 (MQ=255)
                       atTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAAt                                                                                                               <  1:263162/141‑1 (MQ=255)
                        ttGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATc                                                                                                              >  2:200914/1‑141 (MQ=255)
                        ttGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATc                                                                                                              >  2:343675/1‑141 (MQ=255)
                         tGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCa                                                                                                             >  1:151087/1‑141 (MQ=255)
                            cGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAAc                                                                                                          >  1:162755/1‑141 (MQ=255)
                              cAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTa                                                                                                        <  1:235946/141‑1 (MQ=255)
                                    gaggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaac                                                                                                  <  2:421380/139‑1 (MQ=255)
                                      ggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaacaa                                                                                                >  1:309750/1‑141 (MQ=255)
                                      ggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaacaa                                                                                                <  1:226496/141‑1 (MQ=255)
                                          aaTTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTATTGTTACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACttt                                                                                            >  2:331961/1‑141 (MQ=255)
                                              aaGAGCAAAGGAACGACAAGTAGAAACGAGGTAACGACAAATGCTAGAGGCACGATATTTCAATGAAACCGCCCACCATATCGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGt                                                                                        <  1:153930/141‑1 (MQ=255)
                                              aaGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGt                                                                                        <  2:445505/141‑1 (MQ=255)
                                                gagCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTc                                                                                      >  1:231050/1‑141 (MQ=255)
                                                   cAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCtgt                                                                                   >  1:326188/1‑141 (MQ=255)
                                                     aaGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACttt                                                                                            >  2:171495/1‑130 (MQ=255)
                                                     aaGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACttt                                                                                            <  1:171495/130‑1 (MQ=255)
                                                       ggAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGa                                                                               <  1:266757/141‑1 (MQ=255)
                                                          aCGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATgg                                                                            <  1:408859/141‑1 (MQ=255)
                                                          aCGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATgg                                                                            >  1:166574/1‑141 (MQ=255)
                                                               aaGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATgg                                                                       <  2:389769/141‑1 (MQ=255)
                                                                aGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATggg                                                                      <  2:223878/141‑1 (MQ=255)
                                                                    gAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGc                                                                  <  2:103823/141‑1 (MQ=255)
                                                                     aaaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGcc                                                                 >  1:128451/1‑141 (MQ=255)
                                                                     aaaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGcc                                                                 <  1:27868/141‑1 (MQ=255)
                                                                     aaaCGAGGTAAAGAAAAATGATAGAGGAAAAATTTTTTACTGATAATGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTACTCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGcc                                                                 <  2:36669/141‑1 (MQ=255)
                                                                      aaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCa                                                                <  2:323488/141‑1 (MQ=255)
                                                                              aaaGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAg                                                        >  2:337119/1‑141 (MQ=255)
                                                                                aGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGga                                                      >  1:157600/1‑141 (MQ=255)
                                                                                       tGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTg                                               >  2:27486/1‑141 (MQ=255)
                                                                                         cTAGAGGCACAATTTTTTACTGATACTGGACGAGATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                             <  2:169573/141‑1 (MQ=255)
                                                                                          tAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCaa                                            <  2:252875/141‑1 (MQ=255)
                                                                                           agagGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAg                                           <  2:326188/141‑1 (MQ=255)
                                                                                                 acaATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAAtt                                     >  2:354027/1‑141 (MQ=255)
                                                                                                      tttttACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTa                                >  1:136501/1‑141 (MQ=255)
                                                                                                      tttttACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTa                                >  1:195748/1‑141 (MQ=255)
                                                                                                        tttACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACa                              <  2:405666/141‑1 (MQ=255)
                                                                                                           aCTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGAt                           <  1:344736/141‑1 (MQ=255)
                                                                                                                  ctagACAACATAGAGATAAGAAAAAAGATGCGGGTGGTATTTTTTATAATAAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                    <  2:434636/138‑1 (MQ=255)
                                                                                                                            tAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGttt          <  1:23665/141‑1 (MQ=255)
                                                                                                                              gagaTAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTAGTTACAGATGAGTTGAAATCCCGtttgg        >  2:141512/1‑139 (MQ=255)
                                                                                                                                    aGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCgg  >  1:89679/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGG  >  CP000730/1231270‑1231542

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