Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,3410AG20.6% 62.7 / 13.2 34Y283H (TAC→CAC) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base G (5/2);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCTCCGGAAACACATTGTCACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACGCCACCGTAGACCGCCGTCAGGAAGATGGCTTTCGGGAT  >  NC_000913/1358204‑1358476
                                                                                                                                         |                                                                                                                                       
gcgcTCCGGAAACACATTGTCACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACa                                                                                                                                       >  2:845525/1‑140 (MQ=255)
gcgcTCCGGAAACACATTGTCACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACa                                                                                                                                       >  1:757988/1‑140 (MQ=255)
   cTCCGGAAACACATTGTCACGCCCCATCACATACAGCAGACGCGAAACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAgcg                                                                                                                                    <  2:609045/140‑1 (MQ=255)
       ggAAACACATTGTCACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAAt                                                                                                                                >  1:556371/1‑140 (MQ=255)
               aTTGTCACGCCCCAGCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCa                                                                                                                        >  1:87967/1‑140 (MQ=255)
                ttGTCACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAg                                                                                                                       >  1:914697/1‑140 (MQ=255)
                  gTCACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTg                                                                                                                     >  1:254821/1‑140 (MQ=255)
                    cACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTGCAGCGCAATTTCAGGCAGTGCg                                                                                                                   <  2:815026/140‑1 (MQ=255)
                    cACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCg                                                                                                                   <  2:914697/140‑1 (MQ=255)
                           cagcacaTACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCg                                                                                                            <  2:87967/137‑1 (MQ=255)
                            aTCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTGCAGCGCAATTTCAGGCAGTGCGGCATCCgg                                                                                                           >  1:386429/1‑140 (MQ=255)
                                     agcagACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTaaa                                                                                                  >  2:830735/1‑140 (MQ=255)
                                      gcagACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCgg                                                                                                           >  2:384546/1‑130 (MQ=255)
                                      gcagACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCgg                                                                                                           <  1:384546/130‑1 (MQ=255)
                                          aCGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGTTTGCCGCCCACGTACAGCGCAATTTCAAGCAGTGCGGCATCCGGATCTTTAACGCGGc                                                                                             >  2:551125/1‑140 (MQ=255)
                                               acacGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGata                                                                                        >  2:268676/1‑140 (MQ=255)
                                                 acGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTGCAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATc                                                                                      >  2:753131/1‑140 (MQ=255)
                                                            tgAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACa                                                                           <  1:268676/140‑1 (MQ=255)
                                                             gAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAg                                                                          <  1:830735/140‑1 (MQ=255)
                                                             gAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAg                                                                          <  2:541290/140‑1 (MQ=255)
                                                               aGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCt                                                                        >  2:678620/1‑140 (MQ=255)
                                                                    ggccagAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATAGACTGGAAGAGCTTGCCGCCGACGTGCAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATa                                                                   >  1:219422/1‑140 (MQ=255)
                                                                           ggCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAacg                                                            >  2:98759/1‑140 (MQ=255)
                                                                               aaCGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGAcgc                                                        >  1:780441/1‑140 (MQ=255)
                                                                                aCGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGAcgcc                                                       <  1:678620/140‑1 (MQ=255)
                                                                                        acaAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATaaa                                               >  2:239388/1‑140 (MQ=255)
                                                                                          aaaCGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGa                                             >  1:51173/1‑140 (MQ=255)
                                                                                          aaaCGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGa                                             >  2:357780/1‑140 (MQ=255)
                                                                                           aaCGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGAt                                            <  2:757988/140‑1 (MQ=255)
                                                                                               tggtgCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGTAATTTCAGGCAGTGCGACATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCAGAGATAAAGATAACg                                        <  2:661522/140‑1 (MQ=255)
                                                                                               tggtgCAGAGGAAAATAGACTGGAAGAGCTTGCCGCCGACGTGCAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACg                                        <  2:219422/140‑1 (MQ=255)
                                                                                                              tCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACGCCACCGTAGAccgcc                         >  2:972241/1‑140 (MQ=255)
                                                                                                                 aCTGGAAGAGCTTGCCGCCGACGTGCAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACGCCACCGTAGACCGCCGTc                      >  1:65415/1‑140 (MQ=255)
                                                                                                                   tGGAAGAGCTTGCCGCCGACGTGCAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACGCCACCGTAGACCGCCGTCAg                    >  2:71998/1‑140 (MQ=255)
                                                                                                                          gCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACGCCACCGTAGACCGCCGTCAGGAAGATg             <  2:780441/140‑1 (MQ=255)
                                                                                                                                     aCGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACGCCACCGTAGACCGCCGTCAGGAAGATGGCTTTCGGGAt  <  2:51173/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
GCGCTCCGGAAACACATTGTCACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACGCCACCGTAGACCGCCGTCAGGAAGATGGCTTTCGGGAT  >  NC_000913/1358204‑1358476

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