Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,590,3730GC27.3% 32.6 / 14.3 22G284A (GGC→GCC) nicP‑Iporin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (6/10);  new base C (3/3);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTACCAGGTAC  >  NC_002947/3590228‑3590513
                                                                                                                                                 |                                                                                                                                            
tttGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCggg                                                                                                                                          <  2:215233/150‑1 (MQ=255)
       gCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCt                                                                                                                                   <  1:353994/150‑1 (MQ=255)
       gCCTGAAGCACAGCAAACCCGTCCGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGCCATGCCGGGcagtgaa                                                                                                                                   <  1:158315/150‑6 (MQ=255)
                      aaCCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTAcc                                                                                                                    >  1:156976/1‑150 (MQ=255)
                        cccGTCGGCGACTGGGTACTCAGTGCCCACCTGGGTTATTTCCCTAAACCCTGCCAATGCTATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAg                                                                                                                  <  2:116551/150‑1 (MQ=255)
                                       gTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGCCATGCCGGGca                                                                                                                                        >  2:127733/1‑113 (MQ=255)
                                       gTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGCCATGCCGGGca                                                                                                                                        <  1:127733/113‑1 (MQ=255)
                                        tACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATa                                                                                                  <  2:23916/150‑1 (MQ=255)
                                                                                       gcgcTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAg                                                   <  2:156976/150‑1 (MQ=255)
                                                                                                   aaCATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTAt                                       >  2:222895/1‑150 (MQ=255)
                                                                                                                 ggCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTacac                         <  1:222895/150‑1 (MQ=255)
                                                                                                                    gTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGcc                                              >  1:379465/1‑126 (MQ=255)
                                                                                                                    gTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGcc                                              <  2:379465/126‑1 (MQ=255)
                                                                                                                          ctcGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGt                                                              <  2:181005/104‑1 (MQ=255)
                                                                                                                          ctcGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGt                                                              >  1:181005/1‑104 (MQ=255)
                                                                                                                                    tgtGCCCGCCATGCCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGc                               <  1:201403/123‑1 (MQ=255)
                                                                                                                                    tgtGCCCGCCATGCCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGc                               >  2:201403/3‑125 (MQ=255)
                                                                                                                                      tGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTACCAGGt    >  1:267941/1‑150 (MQ=255)
                                                                                                                                      tGCCCGCCATGCCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTACCAGGt    >  2:285817/1‑150 (MQ=255)
                                                                                                                                        cccGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTACCAGGTAc  <  2:369642/150‑1 (MQ=255)
                                                                                                                                         ccGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGAc                                           >  1:187779/1‑108 (MQ=255)
                                                                                                                                         ccGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGAc                                           <  2:187779/108‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
TTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTACCAGGTAC  >  NC_002947/3590228‑3590513

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