Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,224,007	0	C	A	44.0%	84.4 / 25.3	50	intergenic (+100/+272)	ycgI/minE	pseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base C (6/22);  new base A (22/0);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 5.24e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGA  >  NC_000913/1223859‑1224157                                                                                                                                                     |                                                                                                                                                       tgtgGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAcc                                                                                                                                                       <  2:818292/150‑1 (MQ=255)      aaTCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCATcc                                                                                                                                                 >  2:633152/1‑151 (MQ=255)                      ggTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATg                                                                                                                                  <  1:637962/150‑1 (MQ=255)                            tGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTATGGGGGGAAAAAGaataa                                                                                                                             >  2:133291/1‑149 (MQ=255)                            tGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCATCCCAATGGAGAACAAAGaataa                                                                                                                             >  1:184001/1‑149 (MQ=255)                              gtgtGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAACCTCATCCTATGGGAGGAAAAAGaataa                                                                                                                             >  2:109242/1‑147 (MQ=255)                                gtgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAAt                                                                                                                       <  1:421021/151‑1 (MQ=255)                                 tgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACTCATCCTTTGGGAGAAAAAAAAATAAAGaaaat                                                                                                                      >  2:279188/1‑149 (MQ=255)                                  gAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTg                                                                                                                     <  1:913977/151‑1 (MQ=255)                                  gAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTg                                                                                                                     <  1:913834/151‑1 (MQ=255)                                    aCTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGCAACCTCCTCCTATGGGGGAACAAAAAATaattgaattgt                                                                                                                    >  2:746452/1‑149 (MQ=255)                                     cTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCCTCCCATGGGAGAACAATGAATAATGAAATTGCgg                                                                                                                  >  1:296585/1‑151 (MQ=255)                                      tGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACTCATCCTATGGGGGAACAAAGAAAAATGAAAATGCggg                                                                                                                 >  1:852262/1‑151 (MQ=255)                                       ggTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCggg                                                                                                                 <  1:566103/150‑1 (MQ=255)                                        gTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCCATTTGGGAAAAAATAATAATTAAATTTCgggg                                                                                                                >  2:797228/1‑150 (MQ=255)                                           gTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTATGGGGGAACAATGAATAATGAAATTGCgggg                                                                                                                >  2:201420/1‑147 (MQ=255)                                           gTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCCATGGGAGAAAAATGAATAATGAAAATGTGGGGtt                                                                                                              >  2:388142/1‑149 (MQ=255)                                               tGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAACCTCCTCCCATGGGGGGAAAAAGAATAATGAAAATGCGGGGGTatcat                                                                                                         >  2:510930/1‑150 (MQ=255)                                                         cGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACCCATCCTATGGGGGAAAAAAGAAAAATGAAAATGTGGGGGTATCCTCTCCCAGta                                                                                                >  2:396266/1‑149 (MQ=255)                                                         cGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGtata                                                                                              <  2:688887/151‑1 (MQ=255)                                                             tAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACCCATCCCATGGGAGGAAAAAGAAAAATTAAAATTCGGGGGTATCATCTCCCAGTatattc                                                                                           >  2:964297/1‑148 (MQ=255)                                                              aaaCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCAt                                                                                         <  1:76133/151‑1 (MQ=255)                                                              aaaCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCCATGGGGGAAAAAAGAAAAAAGAAATTGCGGGGGTATTATCTCCCACTATATCCa                                                                                          >  2:119738/1‑150 (MQ=255)                                                                aCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATAc                                                                                       <  2:107294/151‑1 (MQ=255)                                                                     gcTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACa                                                                                  <  1:551729/151‑1 (MQ=255)                                                                         cGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAACCTCATCCTATGGGGGAAAAAAAAATAATGAAAATTCGGGGGTTTTATCTCCCCGtataaaccaactaaaaaaa                                                                               >  1:87016/1‑135 (MQ=255)                                                                          ggCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAg                                                                             <  2:887036/151‑1 (MQ=255)                                                                            cgcAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGt                                                                           <  2:667818/151‑1 (MQ=255)                                                                              cAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTATTGGGGAACAATGAATAAAGAAAATGCGGGGGTATCATCTCCCCGTATAACCATACTAACAAAAAgggt                                                                          >  2:677849/1‑148 (MQ=255)                                                                                gTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACTCATCCTCTGGGGGAAAAAAGAATAATGAAATTGGGGGGGTATCATCTCCCAGTATATTCATACTAACAAAAAAGTtttt                                                                       >  1:594073/1‑151 (MQ=255)                                                                                  tttCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAACCCCCTCCTATTGGAGAACAAAGAATAAAGAAAATGCGGGGGTATCATCTCCCCGTATATCCATACTAAAAAAAAGGGTAttt                                                                      >  2:142780/1‑150 (MQ=255)                                                                                    tCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTAc                                                                    <  1:488510/150‑1 (MQ=255)                                                                                    tCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTATGGGGGAACAAAGAAAAAATAAAATTCGGGGGTATTATCTCCCCGTATATACATACTAAAAAAAAAGGTAttt                                                                      >  1:681194/1‑148 (MQ=255)                                                                                       tttACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTATGGGGGAAAAATGAATAAAGAAAATGCGGGGGTATTATCTCCCAGTATATCCATAATAAAAAAAAGGGTAttt                                                                      >  2:38473/1‑145 (MQ=255)                                                                                        ttACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCAGCCTCTGGGAGAACAAATAAAGAAAAAATTGCGGGGGTTTCACCTCGCAACTTATCCATACTTACaaaatggaaattttt                                                                    >  1:233929/1‑133 (MQ=255)                                                                                             ggATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAgg                                                           <  2:618917/150‑1 (MQ=255)                                                                                             ggATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGc                                                          <  2:737942/151‑1 (MQ=255)                                                                                              gATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACTTCTACTCTCCGGGGAAAATGAAAAAAAAAAATCTGTGGGTTTCATCTCCCCGTAAATACAAACTAACAAAAAAGATATTTACtaacaccag                                                           >  1:834142/1‑141 (MQ=255)                                                                                                 gTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATaa                                                      <  2:348621/151‑1 (MQ=255)                                                                                                      aaGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACAtt                                                 <  1:443559/151‑1 (MQ=255)                                                                                                        gggAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCATCCTATGGGGGAACAATGAATAATGAAAATGCGGGGGTTTCATCTCCCAGTATATCCATACTAACAATAAAGTTATTTTCTTAACCAGGCATAAACAattt                                               >  2:351747/1‑151 (MQ=255)                                                                                                            aaaTTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTATGGGGGAAAAATGAATAATGAAAATTCGGGGGTATTATTTCCCAGTATATTCATACTAACAAAAAGGGTATTTTCTCAACCCAGCATAAAAAtttttttt                                           >  2:410353/1‑151 (MQ=255)                                                                                                              aTTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTTCTTCCATGGGGGAAAAAAAAATAATTAAATTGGGGGGGTTTTTTTTCCCCGTATATTCATACTTAAAAAAAAGGTATTTTTTCAACCAAGCATAAAAAatttttttt                                          >  1:715847/1‑150 (MQ=255)                                                                                                                   ggggCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTtgtg                                       <  2:612792/148‑1 (MQ=255)                                                                                                                               gCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCtt                        <  2:833947/151‑1 (MQ=255)                                                                                                                                  cGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAg                     <  2:656520/151‑1 (MQ=255)                                                                                                                                    ccAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAgg                    <  1:682153/150‑1 (MQ=255)                                                                                                                                          aCGGCGGGAAAACTCATCCTATGGGAGAAAAAAGAATAATGAAATTGCGGGGGTATCATCTCCCAGGATATCCATACTAACAAAAAAGTTATTTACTCAACCAGGCATAAACAATTTTTTTTTTGCGTGGGAACAGCCTTTAGGGGTaa               >  1:118800/1‑149 (MQ=255)                                                                                                                                            ggcggGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAggg           >  2:887506/1‑151 (MQ=255)                                                                                                                                                 gAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAgg      <  2:861298/151‑1 (MQ=255)                                                                                                                                                  aaaCCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGt     <  1:448994/151‑1 (MQ=255)                                                                                                                                                     ccTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGa  <  2:285798/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       TGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGA  >  NC_000913/1223859‑1224157

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑