Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0349431,730,0960GC28.8% 67.6 / 14.1 60E173Q (GAA→CAA) rodZcell shape determination protein RodZ
Reads supporting (aligned to +/- strand):  ref base G (18/24);  new base C (1/16);  total (19/41)
Fisher's exact test for biased strand distribution p-value = 5.80e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACGCAAAGTGAAAGCAAGTACGAAATCCCAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGACAGTTCTTGGATTCGAGTACGCGATGAA  >  NZ_CP034943/1729972‑1730233
                                                                                                                            |                                                                                                                                         
tACGCAATTTTAAATCAATGACGAAATCCCATAATATTAATCGGTAAAATATTATCAATGCAACTGTTGTGAAAAAGAGGCATTAACTAAAATAGAAACGAAAGAAAAAGAAGATTAAATAGAAGAAAGCGACAGTGATa                                                                                                                            <  2:1771702/140‑1 (MQ=255)
     aaaGTGAAAGCAAGTACGAAATCCCAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGAGACAGTGATAAACTg                                                                                                                       >  2:45038/1‑140 (MQ=255)
         tGAAAGCAAGTACGAAATCCCAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGgaga                                                                                                                   >  2:1359225/1‑140 (MQ=255)
         tGAAAGCAAGTACGAAATCCCAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGgaga                                                                                                                   >  2:130714/1‑140 (MQ=255)
         tGAAAGCAAGTACGAAATCCCAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGgaga                                                                                                                   >  2:1182653/1‑140 (MQ=255)
             aGCAAGTACGAAATCCCAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAATATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTaa                                                                                                               >  2:571331/1‑140 (MQ=255)
              gaaaGTAATAAAGCCAAAAAGATTCAACGCTAAAACACAAACAACGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAACAAGAAAGCGACAGTGATAAACTGGAGATTaaa                                                                                                              <  2:14186/138‑1 (MQ=255)
                 aGTACGAAATCCCAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCa                                                                                                           >  2:401280/1‑140 (MQ=255)
                 aGTACGAAATCCCAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCa                                                                                                           >  2:1393152/1‑140 (MQ=255)
                   tACTAAACCCAAAAAGTTAAACCGATAAAACACAATAAAACAAACTGCTCTCAAAAAGAGGCACACACTAAAACCCAACCGAAACAAAAACACCATAAAAACCAACACCGCGACAGTGATAAACTGGAGATTAAAGCAAc                                                                                                         <  1:622857/140‑1 (MQ=255)
                   tACGAAATCCCAAACGATTAAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACAATAAAAAATAAACGAAAGAAAAAGAAGATAAAAAAGAAGAACGCGACAGTGATAAACTGGAGATTAAAGCAAc                                                                                                         <  2:929430/140‑1 (MQ=255)
                   tACGAAATCCCAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAAc                                                                                                         >  1:584101/1‑140 (MQ=255)
                          tCCCAAAAGATTCAACGCTAAAAGAGAATAAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAACAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACt                                                                                                  <  2:612373/140‑1 (MQ=255)
                            caaaaaTATCCAACCATAAAACAAAAAAAAACCAACAGCTCTCAAAAACAGCCACAAAATAAAAAACAACCCAAACAAAAACAACAAAACAAACACCAACGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGa                                                                                                <  1:164614/138‑1 (MQ=255)
                             cAAAATATTCAACGATAAAAGATAATAAAAGAAAAAGATATTAAAAAGAGGCAGACACTAAAAAAGAAACGAAACAAAAACAAGATAAAAAACAACAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGaa                                                                                               <  2:187781/140‑1 (MQ=255)
                             cAAAAGTTTCAACGCTAAAACAGAATAAACGCAACCGCTCTGAAAAAGAGGCAGACACTAAAAAACAAACCAAACAAAAATAATATAAAAAAGAACAACGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGaa                                                                                               <  1:920457/140‑1 (MQ=255)
                             cAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAGGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGTTAAAGTGGAGATTAAAGCAACTGGCACTGac                                                                                               >  2:1880505/1‑139 (MQ=255)
                             cAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGaa                                                                                               >  2:591023/1‑140 (MQ=255)
                             cAAAAGATTAAACGATAAAACACAATCAAAGAAACCGATCTCAAAAAGAGCCACACACAAAAAAACAAACCAAAGAAAAACAACATAAAAAACACGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGaa                                                                                               <  1:1484205/140‑1 (MQ=255)
                             aaaaaGATTCAACGATAAAACACAATAAAAGCAACACATCTCAAAAAGAGGAAGACACTAAAAAACAAACGAAACAAAAACAACATAAAAAAGAACACCGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGaa                                                                                               <  1:90331/139‑1 (MQ=255)
                             aaaaaGATTAAACGCTAAAAGACAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAACATAAAAAAGAAGACAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGaa                                                                                               <  2:60300/139‑1 (MQ=255)
                               aaaGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAATAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGTCACAGAAgt                                                                                             >  2:1865634/1‑139 (MQ=255)
                                aaGATTCAACGATAAAACAAAATCAAAGAAACAGATATTAAAAAGAGGAAGACACTAAAAAAGAAACGAAATAAAAAGAACATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAggg                                                                                            <  2:247919/140‑1 (MQ=255)
                                    ttCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCAt                                                                                        >  2:276344/1‑140 (MQ=255)
                                     ccAACGCTAAAACACAATCAAAGAAACAGCTCTTAAAAAGAGGCAAAAACAAAAAAACAAACGAAACAAACACACGATAAAAAACAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCAtt                                                                                       <  2:1653293/139‑1 (MQ=255)
                                        ccgatAAAAGACACTCAACCCAACAGCTCTCAAACCGAGCCCCACACCAACACACAACCCAAACAAACAGAAGACAAAACCGACCACCGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAAc                                                                                    <  1:565992/136‑1 (MQ=255)
                                        actatAAAAGACAATAAAAGAACCAGCTCTTAAAAAGCGGACCACACTAAAAAAGAAACGAAACAAAATGCAGATAAAAAACAAGAACGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAAc                                                                                    <  1:927035/136‑1 (MQ=255)
                                        aCGCTAAAACATAATCAAACCAACAGCTCTGAAAAAGAGGCACACACTAAAAAATAAACGAAACAAAAAGAACATAAAAAAGAAGAACGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAAc                                                                                    <  2:584101/140‑1 (MQ=255)
                                        aCGATAAAAGACAATAAAAGAAACAGCTCTGAAAAAGAGGCACACAATAAAAAAGAAACGAAACAAAAACAAGATAAAAAACAACAACGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAAc                                                                                    <  1:403805/140‑1 (MQ=255)
                                                  acaaTCAAACCAACCGCTCTGAAAAAGAGGAACACACTAAAAAACAAACCAAACAAAAAGAAGATAAAAACCAACAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGaa                                                                          <  1:977537/138‑1 (MQ=255)
                                                  acaaTCAAACCAACACCTCTCAAAAAGAGGCAGACACTAAAAAACAAACCACACAAAAACAAGATAAAAAACAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGaa                                                                          <  2:1711735/138‑1 (MQ=255)
                                                  acaaTAAAAGCAACACAACTCAAAAAGAGCAACAAAATAAAAAACAAACCAAACAAAAACAAGATAAAAAACAACAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGaa                                                                          <  2:1711526/138‑1 (MQ=255)
                                                    aaTCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAgt                                                                        >  1:1797181/1‑140 (MQ=255)
                                                    aaTCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAAATGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAgt                                                                        >  2:497174/1‑140 (MQ=255)
                                                    aaTCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGAAAAACTGGAGATTAAAGCAACTCGCACAGAAGCGTCATTAACGAATTATCaact                                                                        >  2:1134218/1‑138 (MQ=255)
                                                         aaccAACAGCTCTGAAACAGAGGAAGAACAAAAAAAATAAACGACACAAAAACAAGATAAACCACACGAACGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTTTCCg                                                                   <  2:996905/137‑1 (MQ=255)
                                                              aCAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCt                                                              >  2:1696851/1‑140 (MQ=255)
                                                               cAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTg                                                             >  2:36539/1‑140 (MQ=255)
                                                               cAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAACAAAGCGACAGTGATAAACTGGAGATTAAAGCAAATGGCACTGAAGGGTCATTAACTACTTATGAAGTGTCCGGCGCTg                                                             >  2:1293948/1‑140 (MQ=255)
                                                                           agagGCAGACACTAAAAAATAAACGAAAGAAAAAGAAGATAAAAAATAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGc                                                 <  2:972987/140‑1 (MQ=255)
                                                                             agGCAGACACTAAAAAATAAACGAACTAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTc                                               <  2:575906/140‑1 (MQ=255)
                                                                             agGAAGACAGTAAAAAATAAACGAAACAAAAAGAAGATAAAAAAGAAGAACGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTc                                               <  2:1095907/140‑1 (MQ=255)
                                                                                    caatAAAAAATAAACGAAACAAAAAGACTATAAAAAAGAAGAACGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGa                                        <  1:1360300/137‑1 (MQ=255)
                                                                                    caatAAAAAACAAACCAAACAAAAACACGATAAAACAGAACAACGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGa                                        <  1:764318/137‑1 (MQ=255)
                                                                                    aacTAAAAAAGAAACGAAAGAAAAAGACGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGa                                        <  2:1906634/139‑1 (MQ=255)
                                                                                    aacTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGa                                        <  2:1797181/139‑1 (MQ=255)
                                                                                       tAAAAAATAAACCAAAAAAAAAGAAGATAAAAACCAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAgg                                     <  1:577065/140‑1 (MQ=255)
                                                                                         aaaaaGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCa                                   >  2:592309/1‑140 (MQ=255)
                                                                                              taaaCCAAACAAAAATAAGATAAAAAACAACAACGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGa                              <  1:130714/139‑1 (MQ=255)
                                                                                              gAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGa                              <  1:1359225/140‑1 (MQ=255)
                                                                                              gAAACCAAACAAAAACAAGAAAAAAAAGACTAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGa                              <  1:162247/140‑1 (MQ=255)
                                                                                              caaaCGAAAGAAAAAGAAGATAAAAACGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGa                              <  1:1940142/139‑1 (MQ=255)
                                                                                              caaaCGAAACAAAAACAATCTAAATAAGTAGAAAGCGACAGTGATACCCTGGAGATTAAAGCACCTGGCCCTGTAGGGTCATTAACGACGTATCAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGa                              <  2:354024/139‑1 (MQ=255)
                                                                                              caaaCCAAACAAAATTAATATAAAAAATAACAACGCGACATTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGa                              <  1:1182653/139‑1 (MQ=255)
                                                                                              caaaCCAAACAAAAACAAGAAAAAAAACAACACAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGa                              <  1:1258926/139‑1 (MQ=255)
                                                                                              aaaccgtaagaaaacgaaGATAAATACCAAGAACGAGACAGTGATAACCTGGAGATTAAAGCACCTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGGTGAGAAATTCGAGCTCGAACTGAAGGCATCAGa                              <  2:636038/133‑1 (MQ=255)
                                                                                              aaaaCAAAAGAAAAATAAGCTAAAAATGACGAAAGAGACAGTGATAAACTGGAGATTAAAGCAACTGGCAATGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGa                              <  1:1865634/139‑1 (MQ=255)
                                                                                                    aaataaaatcaaTATAAAAACAAACACCGCGACAGTGATAAACTTGAGATTAAAGCAACTGGCACTGATGGGTCATTAAAGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGTACTGAAGGCATCAGACAGTTc                        <  1:497174/140‑1 (MQ=255)
                                                                                                                  tAAAAAAGAACAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGACAGTTCTTGGATTCGAGTAc          <  1:1450156/140‑1 (MQ=255)
                                                                                                                   aaaaaaGAACAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGACAGTTCTTGGATTCGAGTatg         <  1:101953/140‑3 (MQ=255)
                                                                                                                          aagaaAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGACAGTTCTTGGATTCGAGTACGCGATGaa  <  1:591023/140‑1 (MQ=255)
                                                                                                                            |                                                                                                                                         
TACGCAAAGTGAAAGCAAGTACGAAATCCCAAAAGATTCAACGCTAAAAGAGAATCAAAGCAACAGCTCTGAAAAAGAGGCAGACACTAAAAAAGAAACGAAAGAAAAAGAAGATAAAAAAGAAGAAAGCGACAGTGATAAACTGGAGATTAAAGCAACTGGCACTGAAGGGTCATTAACGACTTATGAAGTGTCCGGCGCTGAGAAATTCGAGCTCGAACTGAAGGCATCAGACAGTTCTTGGATTCGAGTACGCGATGAA  >  NZ_CP034943/1729972‑1730233

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