Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,965,557	(T)8→7	coding (74/1266 nt)	rhlB ←	ATP‑dependent RNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,965,550	0	T	.	100.0%	190.2 / NA	47	coding (81/1266 nt)	rhlB	ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (22/25);  total (22/25)

AGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTG  >  NC_000913/3965416‑3965680                                                                                                                                       |                                                                                                                                   agaagaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttttt                                                                                                                                >  1:298946/1‑139 (MQ=255)         aTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAggg                                                                                                                       <  2:260500/139‑1 (MQ=255)         aTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAggg                                                                                                                       <  2:298946/139‑1 (MQ=255)         aTGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAggg                                                                                                                       <  1:334740/139‑1 (MQ=255)          tGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGc                                                                                                                      >  1:215154/1‑139 (MQ=255)                 ttGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACa                                                                                                               >  2:1122733/1‑139 (MQ=255)                    aCGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTa                                                                                                                 >  1:1659701/1‑134 (MQ=255)                    aCGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTa                                                                                                                 <  2:1659701/134‑1 (MQ=255)                          gAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGa                                                                                                      >  2:1106264/1‑139 (MQ=255)                              cGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCa                                                                                                  >  1:105473/1‑139 (MQ=255)                                   tCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCg                                                                                             <  2:1392516/139‑1 (MQ=255)                                    cGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGa                                                                                            <  2:105473/139‑1 (MQ=255)                                            cGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                                                                                    <  2:215154/139‑1 (MQ=255)                                                    gTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctg                                                                            >  1:160528/1‑139 (MQ=255)                                                     tttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                           <  2:712374/139‑1 (MQ=255)                                                         cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttcaag                                                                                                                          <  2:1442243/89‑5 (MQ=255)                                                         cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttcaag                                                                                                                          >  1:1442243/1‑85 (MQ=255)                                                         cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCt                                                                                                                  >  1:119881/1‑96 (MQ=255)                                                         cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCt                                                                                                                  <  2:119881/96‑1 (MQ=255)                                                         cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCg                                                                                                        >  1:250106/1‑106 (MQ=255)                                                         cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCg                                                                                                        <  2:250106/106‑1 (MQ=255)                                                          gcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTtctgttctgt                                                                      <  1:1122733/139‑1 (MQ=255)                                                               cccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACtt                                                                               <  2:1339904/125‑1 (MQ=255)                                                               cccGGCTACGTCACGACCCCCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACtt                                                                               >  1:1339904/1‑125 (MQ=255)                                                                         tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCt                                                       >  1:1370639/1‑139 (MQ=255)                                                                           aCGACCCGCCGGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCa                                                     <  2:702252/139‑1 (MQ=255)                                                                                 cGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCa                                                     <  1:437389/133‑1 (MQ=255)                                                                                 cGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCa                                                     >  2:437389/1‑133 (MQ=255)                                                                                 cGCCAGCGTCACCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtgt                                               <  1:1106264/139‑1 (MQ=255)                                                                                  gCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTa                                                                                                                 <  2:428590/72‑1 (MQ=255)                                                                                  gCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTa                                                                                                                 >  1:428590/1‑72 (MQ=255)                                                                                           aGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAggg                                                                                                                       <  1:1632129/57‑1 (MQ=255)                                                                                           aGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAggg                                                                                                                       >  2:1632129/1‑57 (MQ=255)                                                                                                 aGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGt                                                                                         >  2:138567/1‑81 (MQ=255)                                                                                                 aGGGACAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGt                                                                                         <  1:138567/81‑1 (MQ=255)                                                                                                  gggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTa                              >  1:889949/1‑139 (MQ=255)                                                                                                  gggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTa                              >  2:219818/1‑139 (MQ=255)                                                                                                          gCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTa                                                                                                                 <  2:1329468/48‑1 (MQ=255)                                                                                                          gCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTa                                                                                                                 >  1:1329468/1‑48 (MQ=255)                                                                                                               aaTGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACg                 <  2:160528/139‑1 (MQ=255)                                                                                                                      gTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATa                                                   <  2:839559/98‑1 (MQ=255)                                                                                                                      gTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATa                                                   >  1:839559/1‑98 (MQ=255)                                                                                                                        aCAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCtt                                                                                                          <  1:563673/41‑1 (MQ=255)                                                                                                                        aCAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCtt                                                                                                          >  2:563673/1‑41 (MQ=255)                                                                                                                         cAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATc       <  1:458377/139‑1 (MQ=255)                                                                                                                              aTGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATg                                                                                                    <  2:1022446/41‑1 (MQ=255)                                                                                                                              aTGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATg                                                                                                    >  1:1022446/1‑41 (MQ=255)                                                                                                                              aTGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACtt                                                                               <  2:1684960/62‑1 (MQ=255)                                                                                                                              aTGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACtt                                                                               >  1:1684960/1‑62 (MQ=255)                                                                                                                              aTGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTg  >  1:1094352/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                   AGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTG  >  NC_000913/3965416‑3965680

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