Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,674,1780AC20.7% 56.9 / 14.1 29E89A (GAG→GCG) ‡tqsApheromone AI‑2 transporter
Reads supporting (aligned to +/- strand):  ref base A (10/13);  new base C (3/3);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTG  >  NC_000913/1674032‑1674320
                                                                                                                                                  |                                                                                                                                               
attattTATTGCTGTTATTCTTAACCCGCTGGTGCAACTCATGGTCCGCTGGCGTGGGCCGCGTGTACTGGCGGTGGCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGt                                                                                                                                               <  4:114785/149‑1 (MQ=255)
          gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGAGCGCGGaa                                                                                                                                    >  5:95300/1‑148 (MQ=255)
           cTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGAGCGCGGaa                                                                                                                                    <  6:95300/148‑2 (MQ=255)
             gTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGagcgcg                                                                                                                                       <  1:55359/143‑6 (MQ=255)
             gTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGagcgcg                                                                                                                                       >  2:55359/1‑138 (MQ=255)
                                                                      gCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCtattat                                                                                                         >  4:92012/1‑116 (MQ=255)
                                                                      gCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCtattat                                                                                                         <  3:92012/116‑1 (MQ=255)
                                                                                      tgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                       >  8:16722/1‑118 (MQ=255)
                                                                                      tgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                       <  7:16722/118‑1 (MQ=255)
                                                                                          catcatCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttg                                                                             >  1:194553/1‑124 (MQ=255)
                                                                                          catcatCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttg                                                                             <  2:194553/124‑1 (MQ=255)
                                                                                                  tGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                <  2:51191/129‑1 (MQ=255)
                                                                                                  tGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                >  1:51191/1‑129 (MQ=255)
                                                                                                              tgtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGa                                                                                                       >  2:60880/1‑78 (MQ=255)
                                                                                                              tgtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGa                                                                                                       <  1:60880/78‑1 (MQ=255)
                                                                                                                   cTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATtata                            <  2:166063/148‑1 (MQ=255)
                                                                                                                   cTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATtata                            >  1:166063/1‑148 (MQ=255)
                                                                                                                   cTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CACGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGTTTGACGTCTCAGTTGACCAGCTGGCGCATtata                            <  2:119703/148‑1 (MQ=255)
                                                                                                                     aTTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCaa                                                                                                                  <  6:186733/60‑1 (MQ=255)
                                                                                                                     aTTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCaa                                                                                                                  >  5:186733/1‑60 (MQ=255)
                                                                                                                      ttAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                        >  1:86442/1‑149 (MQ=255)
                                                                                                                         gCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATcc                     <  1:150341/149‑1 (MQ=255)
                                                                                                                         gCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATcc                     >  2:150341/1‑149 (MQ=255)
                                                                                                                                    ttCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGgc                                    <  2:105681/123‑1 (MQ=255)
                                                                                                                                    ttCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGgc                                    >  1:105681/1‑123 (MQ=255)
                                                                                                                                    gtCCGCG‑TCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                        >  4:144346/2‑135 (MQ=255)
                                                                                                                                    gtCCGCG‑TCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                        <  3:144346/134‑1 (MQ=255)
                                                                                                                                          gcTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGt    <  1:170459/149‑1 (MQ=255)
                                                                                                                                              aaCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTg  <  2:86442/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTG  >  NC_000913/1674032‑1674320

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