Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,462,8810C.22.2% 62.7 / 17.3 27intergenic (+235/‑131)fadL/yfdFlong‑chain fatty acid outer membrane transporter/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (9/12);  new base . (3/3);  total (12/15)
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.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCAT  >  NC_000913/2462746‑2463021
                                                                                                                                       |                                                                                                                                            
cAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACGGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCt                                                                                                                                 <  2:51057/149‑1 (MQ=255)
   aCAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTAtt                                                                                                                              <  5:155318/149‑1 (MQ=255)
              gCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCAt                                                                                                                    <  1:169538/148‑1 (MQ=255)
                       ccTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTc                                                                                                                                              >  1:209084/1‑113 (MQ=255)
                       ccTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTc                                                                                                                                              <  2:209084/113‑1 (MQ=255)
                        cTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATATTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAAtt                                                                                                         >  6:151854/1‑148 (MQ=255)
                             aaaaCAATAAGCGGTATACACAGACAAAACTGAAGGGAAAGGCATATATTCAGTCGTTCTGCACGACCTTCTCAAAGAGTTTACTTTTCTGCATATCCAGGATACTCCACCCCCCAGCGTATTGTACGCTCAAACAATCAAATTAAAg                                                                                                     >  6:151400/1‑148 (MQ=255)
                                                acaGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGc                                                                                >  7:94754/1‑149 (MQ=255)
                                                      aaaCCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTc                                                                                                                      <  1:213146/105‑1 (MQ=255)
                                                      aaaCCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTc                                                                                                                      >  2:213146/1‑105 (MQ=255)
                                                          cTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGAt                                                                      <  8:94754/149‑1 (MQ=255)
                                                           tGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATtg                                                                                                                             >  2:48299/1‑94 (MQ=255)
                                                           tGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATtg                                                                                                                             <  1:48299/94‑1 (MQ=255)
                                                                  aaaGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCGGGCTATTGTGCGCTCATCCACTCAAATTAAAGATAGGTTCTGAATAAATGCGCGTTTTTTGATAGTCTAt                                                               >  1:59850/1‑149 (MQ=255)
                                                                               cAGGCCTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATAATCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTaa                                                   <  3:98495/148‑1 (MQ=255)
                                                                                          cTCGTCCTTCTCAAAGAGTTTACTTTTCTTCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATATGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGt                                        <  5:151400/148‑1 (MQ=255)
                                                                                                   ctcAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATCGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAaatcaa                              <  5:151854/148‑1 (MQ=255)
                                                                                                        agagTTTTCTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGGTATTGTGCGCGCATACACTCAAATTAAAGATAGGTTCTAAATAAATGCGCGTTTTTTGATAGGCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCaaa                         <  2:59850/149‑1 (MQ=255)
                                                                                                        agagTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTaa                                                                                                       >  6:16988/1‑71 (MQ=255)
                                                                                                        agagTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTaa                                                                                                       <  5:16988/71‑1 (MQ=255)
                                                                                                            tttACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATCCACTAAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTGTATTAAATTAGGTAATATAAATTTGTAAAAAATCAATCAAAATg                      >  4:117141/1‑148 (MQ=255)
                                                                                                                     cTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCGCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAAt             >  1:127313/1‑148 (MQ=255)
                                                                                                                      tGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAATATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATc            <  3:117141/148‑1 (MQ=255)
                                                                                                                              cAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTaaa                                                                                                      >  5:85211/1‑50 (MQ=255)
                                                                                                                              cAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTaaa                                                                                                      <  6:85211/50‑1 (MQ=255)
                                                                                                                              cAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTAcc    <  2:127313/148‑1 (MQ=255)
                                                                                                                               aGGATACTCCCCCCCCCTGGCGATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATAATGCTACCAt  >  1:86453/1‑149 (MQ=255)
                                                                                                                                 gATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCAtt                                                         >  1:46814/1‑92 (MQ=255)
                                                                                                                                 gATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCAtt                                                         <  2:46814/92‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                            
CAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCAT  >  NC_000913/2462746‑2463021

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: