Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,143,140	T→C	T410A (ACT→GCT)	rne ←	ribonuclease E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,143,140	0	T	C	97.8%	155.4 / NA	46	T410A (ACT→GCT)	rne	ribonuclease E
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (22/23);  total (22/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 = 1.00e+00

CAGCAGGTAAGAAGCGATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCA  >  NC_000913/1143012‑1143255                                                                                                                                 |                                                                                                                    cagcagGTAAGAAGCGATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGaacaac                                                                                                          <  1:774768/140‑1 (MQ=255)  agcagGTAAGAAGCGATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAAcg                                                                                                         >  1:161411/1‑140 (MQ=255)       gTAAGAAGCGATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGAc                                                                                                    >  1:359885/1‑140 (MQ=255)        tAAGAAGCGATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACa                                                                                                   <  1:230642/140‑1 (MQ=255)          agaagCGATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACaaa                                                                                                 <  2:268597/140‑1 (MQ=255)                 aTTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGtgatg                                                                                          >  2:542558/1‑140 (MQ=255)                 aTTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGtgatg                                                                                          >  2:788930/1‑140 (MQ=255)                  ttGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGtgatga                                                                                         <  1:32935/140‑1 (MQ=255)                   tGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGAc                                                                                        <  2:129750/140‑1 (MQ=255)                       acaGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGa                                                                                    <  2:161411/140‑1 (MQ=255)                             aCAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCAcc                                                                              <  1:327695/140‑1 (MQ=255)                             aCAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCAcc                                                                              <  1:638988/140‑1 (MQ=255)                              cAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAGACGTGATGACTGGATTCAccc                                                                             >  2:896277/1‑140 (MQ=255)                                     gTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTc                                                                                 >  2:511374/1‑129 (MQ=255)                                     gTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTc                                                                                 <  1:511374/129‑1 (MQ=255)                                      tGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGAt                                                                                   >  1:168113/1‑126 (MQ=255)                                      tGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGAt                                                                                   <  2:168113/126‑1 (MQ=255)                                             ccTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCGAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAgg                                                              >  2:861040/1‑140 (MQ=255)                                                 ggTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCt                                                          <  1:788930/140‑1 (MQ=255)                                                  gtgtTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTg                                                         >  1:343895/1‑140 (MQ=255)                                                   tgtTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGTGCTGa                                                        <  2:421717/140‑1 (MQ=255)                                                   tgtTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGa                                                        >  2:562371/1‑140 (MQ=255)                                                   tgtTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGa                                                        >  1:142544/1‑140 (MQ=255)                                                    gtTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGAc                                                       <  1:896277/140‑1 (MQ=255)                                                     ttCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACg                                                      >  2:574673/1‑140 (MQ=255)                                                       ctctTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACggg                                                    <  1:861040/140‑1 (MQ=255)                                                       ctctTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACggg                                                    <  1:625471/140‑1 (MQ=255)                                                               gcgcTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACggg                                                    <  2:489069/132‑1 (MQ=255)                                                               gcgcTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACggg                                                    >  1:489069/1‑132 (MQ=255)                                                               gcgcTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTcc                                            >  1:194035/1‑140 (MQ=255)                                                                 gcTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcag                                          <  2:359885/140‑1 (MQ=255)                                                                 gcTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcag                                          <  2:454691/140‑1 (MQ=255)                                                                  cttcttcttCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCaagc                                         >  2:937242/1‑140 (MQ=255)                                                                    tcttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcagcag                                       >  2:214248/1‑140 (MQ=255)                                                                    tcttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcagcag                                       >  2:89225/1‑140 (MQ=255)                                                                       tcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGcc                                    <  1:214248/140‑1 (MQ=255)                                                                                aGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCaa                                  >  2:967154/1‑133 (MQ=255)                                                                                aGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCaa                                  <  1:967154/133‑1 (MQ=255)                                                                                      aGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATg                     <  1:959755/140‑1 (MQ=255)                                                                                                gCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGa           >  1:462744/1‑140 (MQ=255)                                                                                                gCGACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGa           >  2:370968/1‑140 (MQ=255)                                                                                                  gACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAAt         >  2:220459/1‑140 (MQ=255)                                                                                                  gACAGCGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAAt         >  1:14160/1‑140 (MQ=255)                                                                                                       cGATTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATAcgcg    <  2:343895/140‑1 (MQ=255)                                                                                                         aTTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCa  <  1:220459/140‑1 (MQ=255)                                                                                                         aTTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCa  <  2:194035/140‑1 (MQ=255)                                                                                                         aTTCGTTGTCACGCACGGTGCCAGCACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGAGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCa  >  2:895849/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                    CAGCAGGTAAGAAGCGATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCCAAAGCGAGAAATATGGCTGATTTGAATACGCGCA  >  NC_000913/1143012‑1143255

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