Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,340,0890GC32.0% 30.7 / 14.4 26intergenic (+74/‑257)ribAB‑II/PP_38143,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (8/9);  new base C (4/4);  total (12/14)
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.79e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCA  >  NC_002947/4339952‑4340233
                                                                                                                                          |                                                                                                                                                
gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCTCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgcg                                                                                                                                       <  2:372468/149‑1 (MQ=25)
 cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCATTTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                      <  2:102783/149‑1 (MQ=255)
 cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  6:372548/149‑1 (MQ=25)
 cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  2:171481/149‑1 (MQ=25)
 cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                      <  1:15655/149‑1 (MQ=255)
                            cAGAGCATCCCGTTCGACCCGGGGATGGATGGCTCAGAACACACAAGCAATTGATTTAGCGGGGCCCTTTAGGAGCGGGTTTTACCCCCAATGGCTGGGTTACACCCCCCACGCCATACCGGGTAAACCCCCTACTACAGGGCCCGTgg                                                                                                           >  7:365440/1‑149 (MQ=255)
                                                      ggatggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGGAGCGGGTTTAACCCCGAATGCGAGGGTGACACCCCCGACGCATTAGCGGGTAAACCCCCACCTACAGGGCCCTTGGGGGGTTTGACAGGGGGGAGTGAttg                                                                                  >  1:88081/1‑147 (MQ=255)
                                                                caaCCCACAAGCAATTGCTTTAGCGGGGCCCTGGAGGAGCGGGGTTACCCGGTAATGGCCGGGTGACACCACCGACGCATGCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGccc                                                                        <  4:252547/147‑1 (MQ=37)
                                                                caaCACACAAGAAATTTATTTCGCGGGGCCCGGTAGGCGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGcccc                                                                       <  2:88081/148‑1 (MQ=37)
                                                                                          ggCCCTGTAGGCGCGGGTTTACCCGCTAATGCGAGGGTTACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATa                                             <  6:358649/149‑1 (MQ=255)
                                                                                               tGTAGGAGCGGGTTTACCCCCGAATGCGAGGGTGACACCCCCCCCGCATTAGCCGGTAAACCCCCTCCAACAGGGCCCATGGTGGGTTTTACAGGGCGGAGTGATGGCCGGGAGCCCCAGACCACGTGCCAAAAACTTGGAAAACCaa                                         >  2:93578/1‑147 (MQ=19)
                                                                                                                cccGCTAATGCGAGGGTGACACCCCCGACGCATTCGCGGGAAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa                         <  8:102086/148‑1 (MQ=255)
                                                                                                                cccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAAt                                              >  8:163114/1‑127 (MQ=255)
                                                                                                                cccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAAt                                              <  7:163114/127‑1 (MQ=255)
                                                                                                                 ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTTGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                  <  5:166128/122‑1 (MQ=18)
                                                                                                                 ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                       >  5:67840/1‑149 (MQ=25)
                                                                                                                 ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGATTTTGACAGGGAGGAGTGATGGCCGGTCGCCCCATACCTATTGCCAAAATATTg                                                  >  6:166128/1‑122 (MQ=17)
                                                                                                                 ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                       >  2:347073/1‑149 (MQ=255)
                                                                                                                    cgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                    >  8:249433/1‑149 (MQ=255)
                                                                                                                           gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt                   <  8:72112/143‑1 (MQ=255)
                                                                                                                           gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt                   >  7:72112/1‑143 (MQ=255)
                                                                                                                             ggctGACACCACCTCCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCGTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGc           <  8:365440/146‑1 (MQ=255)
                                                                                                                                    accaccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCgg    >  6:22870/1‑149 (MQ=255)
                                                                                                                                      caccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa                 <  2:284549/134‑1 (MQ=255)
                                                                                                                                      caccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa                 >  1:284549/1‑134 (MQ=255)
                                                                                                                                      caccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCa  >  2:364652/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                
GCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCA  >  NC_002947/4339952‑4340233

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