Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913332,9910GT28.9% 103.9 / 21.3 45L207F (TTG→TTT) pdeLDNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL
Reads supporting (aligned to +/- strand):  ref base G (8/24);  new base T (13/0);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 2.79e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAGCCTTCAT  >  NC_000913/332855‑333129
                                                                                                                                        |                                                                                                                                          
tGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCa                                                                                                                                         >  1:517859/1‑140 (MQ=255)
tGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGATGTTTTTTGGCa                                                                                                                                         >  1:905902/1‑140 (MQ=255)
       gTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGAt                                                                                                                                  <  2:826929/140‑1 (MQ=255)
          aTCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCg                                                                                                                                      >  2:370691/1‑133 (MQ=255)
          aTCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCg                                                                                                                                      <  1:370691/133‑1 (MQ=255)
                tCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAag                                                                                                                         <  2:452925/140‑1 (MQ=255)
                 cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  2:99770/140‑1 (MQ=255)
                 cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  2:7086/140‑1 (MQ=255)
                 cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  1:577321/140‑1 (MQ=255)
                 agtATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  1:701731/138‑1 (MQ=255)
                  ttATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGATTTGAAAAaaag                                                                                                                       >  2:667748/1‑140 (MQ=255)
                           tAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAAc                                                                                                              <  2:819255/140‑1 (MQ=255)
                              tGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTTGAAAAAAGTGTCTGaaaccc                                                                                                           >  2:701731/1‑136 (MQ=255)
                                  ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAgtgt                                                                                                                    >  1:568711/1‑127 (MQ=255)
                                  ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAgtgt                                                                                                                    <  2:568711/127‑1 (MQ=255)
                                         cTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGtt                                                                                                        <  2:195680/132‑1 (MQ=255)
                                         cTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGtt                                                                                                        >  1:195680/1‑132 (MQ=255)
                                            atgaAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTTAAAAAGAGGGGCTGGAACTGGTTAATAAATTAgg                                                                                             >  2:333267/1‑140 (MQ=255)
                                            atgaAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAgg                                                                                             >  2:456475/1‑140 (MQ=255)
                                              gaAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTTAAAAAAAGTGTCTGAACCTGGTTAATAAATTAGGTa                                                                                           >  1:457911/1‑140 (MQ=255)
                                                  cAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTATCTGAACCTGGTTAATAAATTAGGTAACGa                                                                                       <  1:737900/140‑1 (MQ=255)
                                                         cGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAGAGTCTGAACCTGGTTAATAAATTAGGGAACGATaaaaaa                                                                                >  2:759553/1‑138 (MQ=255)
                                                           gATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCaa                                                                              <  1:667748/140‑1 (MQ=255)
                                                             tatTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGc                                                                            <  1:759553/140‑1 (MQ=255)
                                                                   gATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAAAGGGTCTGAACCTGGGTAATAAATTTGGGAACGATAAAATCAAGCTGGTtc                                                                      >  1:151364/1‑140 (MQ=255)
                                                                   gATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTTGAAAAAAGGGTCTGAACCTGGGTAATAAATTAGGTAAAGATAAAAACAAACTGGttt                                                                      >  2:726760/1‑139 (MQ=255)
                                                                   gATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGATTTTAAAAAAAGTGTCTTAACCTGGTTAATAAATTAGGGAACGATAAAATCAAACTgggtt                                                                      >  2:43592/1‑137 (MQ=255)
                                                                      gCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAGAGGGTCTGAAACTGGGTAATAAATTAGGGAACGATAAAAACAAACTGGTtctcc                                                                   >  2:450947/1‑139 (MQ=255)
                                                                      gCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTGGCAGCGGGATTTGAAAAAGAGGGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCg                                                                   >  1:828082/1‑140 (MQ=255)
                                                                        cGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAg                                                                 >  1:55063/1‑140 (MQ=255)
                                                                        cGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAAAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGGTCTCGAg                                                                 >  1:259809/1‑140 (MQ=255)
                                                                            aaaaCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTaa                                                             <  2:55063/140‑1 (MQ=255)
                                                                              aaCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACg                                                           <  2:905902/140‑1 (MQ=255)
                                                                                   ttGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACg                                                      <  1:793810/140‑1 (MQ=255)
                                                                                   ttGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACg                                                      <  1:456475/140‑1 (MQ=255)
                                                                                    tgctgcCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAAAGGGGCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACCGAACGt                                                     >  2:931038/1‑140 (MQ=255)
                                                                                       tgcCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAAc                                                  <  1:43592/140‑1 (MQ=255)
                                                                                       tgcCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAAc                                                  <  2:672387/140‑1 (MQ=255)
                                                                                           ggACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGGTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTa                                              >  1:120462/1‑140 (MQ=255)
                                                                                                   tCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTa                                      <  1:809581/140‑1 (MQ=255)
                                                                                                         ttGGCATCAACGTCTCGGCGGGCTGGTTTTTGGTAGCGGGATTTGAAAAATAAGGGTTGAACCTTGATAATAAATTAAGTAACGATAAAAATAAGGAGGTTCTCTAGCCAAATGAACTTTACCCTATTCCGGTAACGCCa                                >  2:260624/1‑140 (MQ=255)
                                                                                                           ggCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGa                              >  1:775394/1‑140 (MQ=255)
                                                                                                            gCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAAAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGaa                             >  2:691706/1‑140 (MQ=255)
                                                                                                             cATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAAAGTGTCTGAACCTGGTTAATAAATTAGGGAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAAAGCCAGAAg                            >  2:726967/1‑140 (MQ=255)
                                                                                                                   cGTCTCGGCGGGGTGTTTTTTTGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCagag                      >  2:237968/1‑140 (MQ=255)
                                                                                                                          gcggGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATAtt               <  1:450947/140‑1 (MQ=255)
                                                                                                                             ggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa            <  2:120462/140‑1 (MQ=255)
                                                                                                                             ggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa            <  2:775394/140‑1 (MQ=255)
                                                                                                                                tgttTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAg         <  1:237968/140‑1 (MQ=255)
                                                                                                                                       tGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAGCCTtcat  >  1:160305/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAGCCTTCAT  >  NC_000913/332855‑333129

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