Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,340,0760GC30.4% 25.2 / 14.3 23intergenic (+61/‑270)ribAB‑II/PP_38143,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (7/9);  new base C (3/4);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTT  >  NC_002947/4339928‑4340216
                                                                                                                                                     |                                                                                                                                            
ggggCCGCCGCTGAACTCCGCGGGGGGGGCGGGCTATGCGCTGGAGGGGGGGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACCCAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAg                                                                                                                                              <  8:115205/149‑1 (MQ=255)
 gggCCGCCGCTGAAATACGGGGGGCTGGCGGGCTATTCGCTGGCGGTGGTGCAGAGCATCCCTTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAg                                                                                                                                              <  8:120840/148‑1 (MQ=255)
  ggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAggg                                                                                                                                            <  7:235330/149‑1 (MQ=255)
       ccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGaca                                                                                                                                       <  1:91935/149‑1 (MQ=255)
            gAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGCGGATGGATGGCTGAGAACACACAAGCAAGTGATTTAGCGGGGCCCTGTAGGAGCGGCTTTACCCCCGAAAGCGACGGTGACcccacc                                                                                                                                  >  3:252547/1‑149 (MQ=255)
                        gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCTCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgcg                                                                                                                      <  2:372468/149‑1 (MQ=25)
                         cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCATTTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                     <  2:102783/149‑1 (MQ=255)
                         cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                     <  2:171481/149‑1 (MQ=25)
                         cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                     <  6:372548/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                             <  7:163114/127‑1 (MQ=255)
                                                                                                                                        cccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAAt                             >  8:163114/1‑127 (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  >  7:72112/1‑143 (MQ=255)
                                                                                                                                                   gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt  <  8:72112/143‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GGGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTT  >  NC_002947/4339928‑4340216

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: