Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,753,986	0	A	T	23.1%	97.0 / 27.1	52	intergenic (+40/+22)	bamE/ratB	lipoprotein component of BamABCDE OM biogenesis complex/UPF0125 family protein
Reads supporting (aligned to +/- strand):  ref base A (22/18);  new base T (6/6);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCC  >  NC_000913/2753838‑2754122                                                                                                                                                     |                                                                                                                                         ccAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCaa                                                                                                                                          <  5:28197/149‑1 (MQ=255)            gTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCttt                                                                                                                               <  3:82907/149‑1 (MQ=255)                 gAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtct                                                                                                                          >  6:153386/1‑149 (MQ=255)                   aGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtctgt                                                                                                                        <  1:110574/149‑1 (MQ=255)                      tgtAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttttaacgtct                                                                                                                     <  5:1014/149‑8 (MQ=255)                         aaCTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCtttttt                                                                                                                            >  3:74151/1‑139 (MQ=255)                         aaCTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCtttttt                                                                                                                            <  4:71828/139‑1 (MQ=255)                         aaCTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCtttttt                                                                                                                            <  4:74151/139‑1 (MQ=255)                         aaCTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGATCATTGAGCACCtttttt                                                                                                                            >  3:71828/1‑139 (MQ=255)                                  aaCGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCt                                                                                                                                 >  1:93252/1‑125 (MQ=255)                                  aaCGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCt                                                                                                                                 <  2:93252/125‑1 (MQ=255)                                       tgacgctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGAtt                                                                                                    >  4:11765/1‑149 (MQ=255)                                             tgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAACAGGTGCTCAATGAGCACCttttttaacgt                                                                                                                       >  2:102226/1‑119 (MQ=255)                                             tgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAACAGGTGCTCAATGAGCACCttttttaacgt                                                                                                                       <  1:102226/124‑6 (MQ=255)                                               acCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTg                                                                                              <  5:153386/147‑1 (MQ=255)                                                   ttAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGtt                                                                                        >  1:138041/1‑149 (MQ=255)                                                          tAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGc                                                                                  >  5:27717/1‑148 (MQ=255)                                                                    ttGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                       >  6:134433/1‑149 (MQ=255)                                                                      gACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                         >  1:177939/1‑81 (MQ=255)                                                                      gACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                         <  2:177939/81‑1 (MQ=255)                                                                        ccAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTAGTTATTCGCTGATTTTTCTGCTCGTTGCCTg                                                                                   >  8:29900/1‑133 (MQ=255)                                                                        acAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCACAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGAt                                                                   <  6:27717/148‑1 (MQ=255)                                                                             atTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGc                                                               >  7:97308/1‑148 (MQ=255)                                                                                gATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAAt                                                            >  5:34667/1‑148 (MQ=255)                                                                                gATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAAt                                                            <  2:138041/148‑1 (MQ=255)                                                                                  tAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                         >  2:187257/1‑69 (MQ=255)                                                                                  tAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                         <  1:187257/69‑1 (MQ=255)                                                                                     cAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAgg                                                      <  8:97308/149‑1 (MQ=255)                                                                                                    tGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATgaga                                                        <  4:115107/132‑1 (MQ=255)                                                                                                    tGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATgaga                                                        >  3:115107/1‑132 (MQ=255)                                                                                                    tGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCAc                                       >  7:92565/1‑149 (MQ=255)                                                                                                    tGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCAc                                       <  6:34667/149‑1 (MQ=255)                                                                                                    tGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCAc                                       <  3:11765/149‑1 (MQ=255)                                                                                                        aaCTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGa                                                  <  5:165782/134‑1 (MQ=255)                                                                                                        aaCTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGa                                                  >  6:165782/1‑134 (MQ=255)                                                                                                           taataataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttt                                                                                                                             <  2:54771/56‑1 (MQ=255)                                                                                                           taataataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttt                                                                                                                             >  1:54771/1‑56 (MQ=255)                                                                                                           taataataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCg                                >  3:130620/1‑149 (MQ=255)                                                                                                                 taaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATg                                                           <  2:71332/116‑1 (MQ=255)                                                                                                                 taaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATg                                                           >  1:71332/1‑116 (MQ=255)                                                                                                                     gTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCt                                                                                               >  8:56458/1‑76 (MQ=255)                                                                                                                     gTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGca                      <  4:130620/149‑1 (MQ=255)                                                                                                                        gCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtctgtct                                                                                                                      <  8:55677/50‑1 (MQ=255)                                                                                                                        gCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtctgtct                                                                                                                      >  7:55677/1‑50 (MQ=255)                                                                                                                              aaaGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGc             >  2:21094/1‑149 (MQ=255)                                                                                                                                     ttAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGc                                                                                             >  5:163061/1‑62 (MQ=255)                                                                                                                                     ttAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGc                                                                                             <  6:163061/62‑1 (MQ=255)                                                                                                                                       aaaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCca                             <  4:108148/124‑1 (MQ=255)                                                                                                                                       aaaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCca                             >  3:108148/1‑124 (MQ=255)                                                                                                                                         aaaaGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGcc  >  2:40555/1‑149 (MQ=255)                                                                                                                                               tGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCCGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                       >  2:155945/1‑74 (MQ=255)                                                                                                                                                gCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCt                                                                                    <  2:180484/60‑1 (MQ=255)                                                                                                                                                gCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCt                                                                                    >  1:180484/1‑60 (MQ=255)                                                                                                                                                     |                                                                                                                                         CCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCC  >  NC_000913/2753838‑2754122

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