Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,340,0750AT31.7% 38.5 / 25.8 45intergenic (+60/‑271)ribAB‑II/PP_38143,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (12/16);  new base T (7/6);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTG  >  NC_002947/4339929‑4340223
                                                                                                                                                   |                                                                                                                                                       
gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAgg                                                                                                                                                       <  6:34117/149‑1 (MQ=255)
gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAgg                                                                                                                                                       <  5:161188/149‑1 (MQ=255)
gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAgg                                                                                                                                                       <  5:147786/149‑1 (MQ=255)
    cgccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGa                                                                                                                                                   <  6:360564/149‑1 (MQ=255)
       cgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTggtgt                                                                                                                                                >  7:167559/1‑145 (MQ=255)
       cgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGacac                                                                                                                                                <  4:229230/149‑1 (MQ=255)
          tGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTggtgtcac                                                                                                                                             >  8:320788/1‑142 (MQ=255)
               tACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACg                                                                                                                                        <  3:232480/149‑1 (MQ=255)
                      ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGTTATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                 <  1:313695/149‑1 (MQ=25)
                      ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                 <  8:167559/149‑1 (MQ=25)
                                cTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                 >  6:300433/1‑139 (MQ=25)
                                cTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                 <  5:300433/139‑1 (MQ=25)
                                   tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                 <  6:246887/136‑1 (MQ=25)
                                   tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                 >  5:246887/1‑136 (MQ=25)
                                   tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCCGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgcg                                                                                                                               <  3:10259/137‑1 (MQ=255)
                                   tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCCGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgcg                                                                                                                               >  4:10259/1‑137 (MQ=255)
                                     aGCTGGCGGGGGGGCAGAGCATCCCGTGCGACCCGGGTTTGGATGGCTGATAACACACAAGCAATTGATTTAGGGGGGCCCTGTTGGAGGGGGTTTACCCGCGAATGCGAGGGTGACACCACCTACGCATTCGCGGGT‑‑AAACCCGCTc                                                                                                                <  2:88243/148‑1 (MQ=37)
                                                        cATCCCGTGCGACCGGGGGATGTATGGCTGAGAAAACACACGAAATTTATTTTGCGGGGCCCTGTAGGAGGGGGTTTACCCGCTAAGGCGAGGGTGACACCACCGACGCTTTCGCGGTT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCtt                                                                                            <  2:291259/149‑1 (MQ=25)
                                                          tCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTAGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCtttt                                                                                          >  1:300460/1‑149 (MQ=255)
                                                            ccGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGTCACCCCCGCCGCATTAGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGa                                                                                        >  1:24356/1‑149 (MQ=255)
                                                                 cGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCCCGAATGCGAGGGGGACACCACCCCCGCCTTAGCGGGT‑‑AAACCCGCTCCTACACGGCCCGAGGTGCGTGTGACAggg                                                                                   >  8:283493/1‑149 (MQ=25)
                                                                  gACCGGGGTATGGATGGCTGATAACACACAAGCAATTTTTTTCGCGGGGCCCTGTAGGAGCGGTTTTACCCGCTAATGCGAGGGTTACACCACCGACGCATTCGCGGTT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAggg                                                                                   <  2:24356/148‑1 (MQ=37)
                                                                   aCCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTACCCCCAAAGCGAGGGTGTGCCCCCCGCCCCAATAGCGGGG‑‑AAACCCCCGGCGAAACGGCCCTTGGGGGTATTTACACGGCg                                                                                 >  6:160442/1‑149 (MQ=255)
                                                                             ggatggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTCCCCCGAAATCGAGGGTGACAACACCGCCCCATTAGCGGGGT‑‑AAACCCCACCTACAGGGCCCATGGTGCTTTTGACAgg                                                                                    >  1:73706/1‑136 (MQ=255)
                                                                              gatggCTGACAACACACACCCAATTGATTTCGCGGGGCCCTGTCGGCGGGGGTTTCCCCGCTAATGCGCGGGTTACACCACCGACGCATTCGGGGTT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGc                                                                       <  2:300460/148‑1 (MQ=37)
                                                                                      acaacaaacaAGCATTTTTTTTCGCGGCGCCGTGTGGGGGGGGTTTACCCGGCTATGGCGCGGTTCACCCCCCCGCCGATTGCGGGTTAAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGccc                                                              <  8:173545/147‑1 (MQ=255)
                                                                                           acacaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATAc                                                         >  8:283039/1‑149 (MQ=255)
                                                                                               aaGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGc                                                                                                                  >  5:357792/1‑88 (MQ=255)
                                                                                               aaGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGc                                                                                                                  <  6:357792/88‑1 (MQ=255)
                                                                                                 gCAATTTATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCTAATGCTAGGGTGACACCACCTACGCATGCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTg                                                   <  7:299668/149‑1 (MQ=38)
                                                                                                        aTTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTg                                                                            <  5:362005/117‑1 (MQ=255)
                                                                                                             gCGGGGCCCTGTAGGCGCGGGTTAACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                        <  7:283493/148‑1 (MQ=255)
                                                                                                             gCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                        <  7:283039/148‑1 (MQ=255)
                                                                                                             gCGGGCCCCTGTTGGTGGGGTTATACCGGGTAATGCGCGGGTGACACCACCGCCTCTTTCGGGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                        <  1:311827/148‑1 (MQ=255)
                                                                                                             gCGGGCCCCTGTCGGTGCGGTTTTACCCGCTAATCCGCGGGTGACACCACCGACGCATTCGCGGTT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                        <  5:12948/148‑1 (MQ=38)
                                                                                                                         cggcgCGGGTTTACCCGCTAATGCGTCGGTGGTGCCACCCTCGCATTCGCGGTT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGaa                                     <  2:298490/135‑1 (MQ=21)
                                                                                                                         aGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGGTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGaa                                     >  1:298490/1‑139 (MQ=18)
                                                                                                                                 gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATcc                   >  3:241188/1‑149 (MQ=255)
                                                                                                                                       ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa             >  3:339078/1‑149 (MQ=25)
                                                                                                                                       ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa             >  1:76081/1‑149 (MQ=255)
                                                                                                                                         cgtaattGCGTGGGTCACACCACCGACGCATTGGGGTGT‑‑AACCCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc            <  5:207819/143‑1 (MQ=255)
                                                                                                                                          gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc           >  8:213613/1‑149 (MQ=25)
                                                                                                                                          gcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc           >  5:349898/1‑149 (MQ=255)
                                                                                                                                          gcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc           >  2:148170/1‑149 (MQ=255)
                                                                                                                                           cgAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt          <  2:76081/149‑1 (MQ=255)
                                                                                                                                                   aGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTg  >  3:351310/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTG  >  NC_002947/4339929‑4340223

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: