Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,020,956	A→G	intergenic (‑37/‑182)	sulA ← / → sxy	SOS cell division inhibitor/CRP‑S‑dependent promoter expression factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,020,956	0	A	G	100.0%	145.2 / NA	45	intergenic (‑37/‑182)	sulA/sxy	SOS cell division inhibitor/CRP‑S‑dependent promoter expression factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/20);  total (25/20)

TGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTA  >  NC_000913/1020820‑1021095                                                                                                                                         |                                                                                                                                            tGATAAGCCCGGCTGTAGTGTTTTCCGTAGAAACACGCGCAATTTTACTTGCTGCGGATGAGAAAGACGAAGAACGATGTGCATAGCCTGAAGGGTACATAATCAATCAAGCCCCTGAGAGTTACTGTATGGATgtgcagtac                                                                                                                                       >  2:600259/1‑143 (MQ=255)    tAAGCCTGTTTGTCGTGTCTTCCGCAAGTACACGCGCAATTTTCCTTGCTACGGATGACAACGACAAAGATCGATGTGCATAGCCTGAAGTGTACCTATTCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATc                                                                                                                                    <  2:67942/143‑1 (MQ=255)     atgCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATcc                                                                                                                                   <  2:152036/141‑1 (MQ=255)       gCCCGGCTGTAGTGTTTTCCGGAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAg                                                                                                                                 <  2:22890/143‑1 (MQ=255)        cccGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCact                                                                                                                                >  2:1144145/1‑141 (MQ=255)        cccGGCTGTAGTGTTTGCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGt                                                                                                                                <  1:1059134/143‑1 (MQ=255)              tGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGacaac                                                                                                                          >  1:1015488/1‑143 (MQ=255)               gTAGTGTGTTCCGGAGAGACACGCGCAATTGTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGacaaca                                                                                                                         <  2:212329/143‑1 (MQ=255)                tAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGAcaacaa                                                                                                                        <  1:898194/143‑1 (MQ=255)                       ttCCGTAGAGACACGCGCCATTTTACTTTCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAATTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCaa                                                                                                                 <  2:696173/143‑1 (MQ=255)                            tAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTa                                                                                                            >  1:953500/1‑143 (MQ=255)                                     gcgcAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaa                                                                                                   <  1:795042/143‑1 (MQ=255)                                     gcgcAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGAGAGTTACTGTATCGCTGTGCAGTACATCCAGTGACAACAAAGAACAACCCTAttatagaaa                                                                                                   >  2:1054146/1‑136 (MQ=255)                                         aaTTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAagag                                                                                               <  1:887746/143‑1 (MQ=255)                                         aaTTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAagag                                                                                               >  1:1108104/1‑143 (MQ=255)                                                 ttGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCAGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCaa                                                                                       >  2:473430/1‑143 (MQ=255)                                                     tgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAttt                                                                                   >  1:883158/1‑143 (MQ=255)                                                            gagaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGtt                                                                            <  1:771352/143‑1 (MQ=255)                                                               aaCGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCACCCCTATTTTCGGAAAGCGACTCGCACCTTTCGTCGTTggc                                                                         >  2:1026298/1‑142 (MQ=255)                                                                 cgacgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCGCGCAAATTTTGTCGTTGGTGa                                                                       >  2:635960/1‑143 (MQ=255)                                                                      aagaaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAACGATCCACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGa                                                                  >  1:902679/1‑143 (MQ=255)                                                                       agaaCGATGTGCATAGCCTGAAGTGTACATAATCAAGCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaa                                                                 <  2:635287/143‑1 (MQ=255)                                                                        gtaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaa                                                                <  1:890154/141‑1 (MQ=255)                                                                        gaaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAGCCACGGTTAACCCTAGCTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaa                                                                >  2:820202/1‑143 (MQ=255)                                                                           cGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACa                                                             >  2:670156/1‑143 (MQ=255)                                                                              tgtgCATAGCCTGAAGTGTACATAATCAATCCAGAGCATGCGGGCGACTGTATCGATGTGCAGTACATCCAGTGACAACAAAGATCAACCATATATTCGGAAAGAGCCTCGCAAATTTTGTCGTTCATGACGGGCAAACAtca                                                          >  2:248511/1‑141 (MQ=255)                                                                                   aTAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGCACGTGCCTCGCGTATTTTGTGGTTGGTGACCGGAAAACATAAATTaa                                                     >  1:96283/1‑143 (MQ=255)                                                                                    tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt                                                    <  2:372399/143‑1 (MQ=255)                                                                                    tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt                                                    <  1:701809/143‑1 (MQ=255)                                                                                       ccTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCACACTCTGTCGTTGGTGACGGGAAAACATAAATTAATCtt                                                 <  1:599991/143‑1 (MQ=255)                                                                                       ccTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGGAAACATAAATTAATCtt                                                 >  1:1072278/1‑143 (MQ=255)                                                                                           aaGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGccc                                             >  2:1016416/1‑143 (MQ=255)                                                                                               gtACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTa                                         >  2:32362/1‑143 (MQ=255)                                                                                                tACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa                                        <  1:672227/143‑1 (MQ=255)                                                                                                   aTAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTACGaa                                     >  1:379216/1‑143 (MQ=255)                                                                                                    tAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAAt                                    <  1:736332/143‑1 (MQ=255)                                                                                                      atcaatcCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATaa                                  >  2:75066/1‑143 (MQ=255)                                                                                                      atcaatcCAGCCCCTGTGAGTTACTGTATGGAAGTGCAGTACATCCAGTGAAAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAGATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGTATaa                                  >  2:603510/1‑143 (MQ=255)                                                                                                       tcaatcCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAg                                 <  2:479495/143‑1 (MQ=255)                                                                                                          atcCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTg                              >  2:816918/1‑143 (MQ=255)                                                                                                           tcCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGATAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGc                             <  1:748188/143‑1 (MQ=255)                                                                                                             cAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCt                           >  1:404121/1‑143 (MQ=255)                                                                                                                          ttACTGTATGGACGTGCCGTACATCAAGTGACAACGAAGATCAACCCTGTTTTCGGAAAGAGCCTCGCATTTTTTCTCGTTGATGACGGGATAACATAGATTGATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTa              <  2:764158/143‑1 (MQ=255)                                                                                                                                tATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGAt        >  2:492625/1‑143 (MQ=255)                                                                                                                                tATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTCCCTATGTTCGTAGTTAAccaag        >  2:82990/1‑139 (MQ=255)                                                                                                                                 aTGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATc       >  1:620470/1‑143 (MQ=255)                                                                                                                                      tgtgCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATATATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTa  <  2:724514/143‑1 (MQ=255)                                                                                                                                         |                                                                                                                                            TGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTA  >  NC_000913/1020820‑1021095

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