Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,340,0890GC28.1% 43.2 / 11.5 32intergenic (+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 (15/8);  new base C (5/4);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCACCT  >  NC_002947/4339951‑4340236
                                                                                                                                          |                                                                                                                                                     
ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                             <  5:100915/149‑1 (MQ=25)
 gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgcg                                                                                                                                            >  2:19501/1‑149 (MQ=25)
   tGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCATGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                          >  6:185895/1‑149 (MQ=255)
     gCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGAGGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCATGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                          <  5:185895/147‑1 (MQ=255)
      cGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTACCCCCAATTCGGGGGTGGCACCACCGACGCCTTAGCgggga                                                                                                                                        >  8:152717/1‑146 (MQ=255)
                       gtggtgCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTAACCCGAATGGGGGGGTGTGACACCCCACGCATTGGCGGGGTAACCCCCGCCTaaacggg                                                                                                                     >  3:24506/1‑143 (MQ=255)
                                     cccGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTAACCGCGAAAGCGAGGGTGACACCACCGCCCCATTCGCCGGTAAACCCCCGCCTAAACGGGCCTTGGTGCgttt                                                                                                        >  3:237643/1‑148 (MQ=255)
                                      ccGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                           <  4:103675/113‑1 (MQ=18)
                                      ccGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                           >  3:103675/1‑113 (MQ=18)
                                       cGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTACCGCGCATTGGGGGGGGGCACCCCCCACCCATTAGCGGGGTAAACCGCGCCGACACGGCCCGTGGTGCgttttt                                                                                                      >  5:281416/1‑146 (MQ=255)
                                        gTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTCCCGCGAATGCGACGGTGCCACCACCGACGCCTTAGCGGGTAAACCCGCTCCTACAGGGCCCGTGGTGCGTTTGa                                                                                                      >  7:80310/1‑147 (MQ=37)
                                               ccGGGGCTGGATGGCTGACAACACCCACGCAATTGATTTCGCGGGGCCCGGTAGGAGCGGGTTAACCCGCTACTGCGAGGGGGACACCACCCCCGCATGCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGa                                                                                             <  5:254721/149‑1 (MQ=25)
                                                 ggggATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTCCCGCGCATGCGAGGGTGACACCACCGACGCATTAGCGGGTAAACCCGCGCCTAAAGGGCCCATGGTGGTTTTGACAGGGGGGAGt                                                                                           >  6:70740/1‑149 (MQ=255)
                                                                acaacaaacaACAATTTGATTTCGGGGGCCCGTGTCGGCGGGGGTTACCCCGCTATGGCGAGGGGCACCCCCCCCCCCAATGCGGGGTAAACCCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGc                                                                             <  7:152717/145‑1 (MQ=255)
                                                                     acacaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTAACCGCGAATGCGAGGGGGGCCCCACCCACGCCTTAGCGGGTAAAACCGCGCCTACAGGGCCCATGGTGCGGTTTACAGGGGGGAGTGATTGGCCGGCGCCCCAAAc                                                                      >  2:220658/1‑149 (MQ=18)
                                                                                        cGGGGCCCTGTAGGAGCGGGTTTAACCCCGAATGCGAGGGTGACACCACCGACGCATTAGCGGGTAAACCCGCTACTACCGGGCCCGTGGTGGGTTTGACAGGGCCGAGTGAAGGCCCGTCGCCCCATACCTAGTTCCAAAATCTTgg                                                    >  2:52970/1‑148 (MQ=25)
                                                                                              ccTGTAGGAGCGGGTTTAACCCCGAATTCGAGGGTGACACCCCCGACGCATTCGCCGGTAAACCCCCTCCAACAGGGCCCTTGGTGCTTTTTACAAGGCGGAGTGATGGCCGGTAGCCCCAAAACACATGCCAAAAAATTTGAATaac                                              >  6:148871/1‑146 (MQ=255)
                                                                                                    ggcgggggTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCataata                                        <  5:148871/143‑1 (MQ=255)
                                                                                                      aGCGGTTTACCCCGCTAATGCGAGGGTGACACCACCGACGAATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGa                                     <  1:52970/149‑1 (MQ=39)
                                                                                                               aCCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCaaaa                                                          >  8:256141/1‑119 (MQ=255)
                                                                                                                 ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                          >  5:113660/1‑149 (MQ=255)
                                                                                                                 ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                          >  4:195553/1‑149 (MQ=255)
                                                                                                                       ttGCTAGGGTGACACCACCGACGCATTGGGGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGGCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa                    <  1:220658/148‑1 (MQ=255)
                                                                                                                       ttGCGAGGGGTCCACCACCGACCCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa                    <  4:24506/148‑1 (MQ=255)
                                                                                                                        tGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTCCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCCTACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCACGTAGACCTTGAg                   >  7:120407/1‑149 (MQ=255)
                                                                                                                           gAGGGTGACACCACCGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCGTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGaa                                                  >  4:39395/1‑115 (MQ=255)
                                                                                                                           gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCt                >  2:234232/1‑149 (MQ=255)
                                                                                                                              ggTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGc              <  6:281416/148‑1 (MQ=255)
                                                                                                                                  acacCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGcc         <  5:70740/149‑1 (MQ=255)
                                                                                                                                      caccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGTTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACAATAATTTGAAATCCCGTAGACATTGAGTCTGCAGGCCGGCa     >  8:105641/1‑149 (MQ=255)
                                                                                                                                       accGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCAc    >  6:207240/1‑149 (MQ=255)
                                                                                                                                         cGACGCATTCGGGGGTAAACCCGCTCCTACAGGGCCCGTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCACCt  <  4:237643/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
GGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCACCT  >  NC_002947/4339951‑4340236

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