Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,9170AT29.2% 27.3 / 16.3 24intergenic (‑192/+18)uof/fldAryhB‑regulated fur leader peptide/flavodoxin 1
Reads supporting (aligned to +/- strand):  ref base A (8/9);  new base T (3/4);  total (11/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 = 4.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGG  >  NC_000913/710778‑711055
                                                                                                                                           |                                                                                                                                          
aataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgc                                                                                                                                    <  4:65241/148‑1 (MQ=255)
  taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgccg                                                                                                                                  >  3:62750/1‑148 (MQ=255)
   aaaaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                                                                                                                >  6:100105/1‑149 (MQ=255)
     aaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTACGCCGtat                                                                                                                               >  3:140035/1‑148 (MQ=255)
     aaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGt                                                                                                                                          <  6:131101/137‑1 (MQ=255)
         gTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTACGCCGtattaata                                                                                                                          <  4:140035/149‑6 (MQ=255)
                     aCATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTga                                                                                                                <  3:75570/147‑1 (MQ=255)
                     aCATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTga                                                                                                                <  5:100105/147‑1 (MQ=255)
                                   aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATg                                                                                                >  4:31268/1‑149 (MQ=255)
                                        cACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgc                                                                                                                                    <  4:124064/108‑1 (MQ=255)
                                        cACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgc                                                                                                                                    >  3:124064/1‑108 (MQ=255)
                                              aaCACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCa                                                                                      <  3:31268/148‑1 (MQ=255)
                                                          ttGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTaa                                                                         <  4:62750/149‑1 (MQ=255)
                                                                      aacaGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTa                                                             <  1:133562/149‑1 (MQ=255)
                                                                         aGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTacac                                                          >  3:118929/1‑149 (MQ=255)
                                                                                           gTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTcc                                         <  4:118929/148‑1 (MQ=255)
                                                                                                          aTTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtcgt                         >  2:70259/1‑149 (MQ=255)
                                                                                                                     cTATTAAAACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCAATGAGAATTTCGTAGAGATGCAACTCTTCTGAAATCTGATTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGa              >  8:63524/1‑149 (MQ=255)
                                                                                                                            tgcGGCGTAGACTCATGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGa                                   >  5:39600/3‑121 (MQ=255)
                                                                                                                            tgcGGCGTAGACTCATGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGa                                   <  6:39600/119‑1 (MQ=255)
                                                                                                                                        aTGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACg                                 >  5:46190/1‑111 (MQ=255)
                                                                                                                                        aTGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACg                                 <  6:46190/111‑1 (MQ=255)
                                                                                                                                          gAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTgg  >  2:178121/1‑140 (MQ=255)
                                                                                                                                          gAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTgg  <  1:178121/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                          
AATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGG  >  NC_000913/710778‑711055

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: