Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,073,603	T→C	T137A (ACT→GCT)	rutA ←	pyrimidine oxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,073,603	0	T	C	100.0%	154.3 / NA	45	T137A (ACT→GCT)	rutA	pyrimidine oxygenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/20);  total (25/20)

TAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCC  >  NC_000913/1073468‑1073738                                                                                                                                        |                                                                                                                                        tAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACg                                                                                                                                     <  2:183928/140‑1 (MQ=255)         cTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGAc                                                                                                                             >  2:787845/1‑140 (MQ=255)              ccccGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCaa                                                                                                                        >  1:947059/1‑140 (MQ=255)                  gTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAAcg                                                                                                                    <  2:260576/140‑1 (MQ=255)                             gTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAAcgcc                                                                                                                  >  1:983200/1‑131 (MQ=255)                             gTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAAcgcc                                                                                                                  <  2:983200/131‑1 (MQ=255)                                gcgcAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATgg                                                                                                      >  1:596707/1‑140 (MQ=255)                                   cagcaCCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATGCCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGt                                                                                                   >  2:56277/1‑140 (MQ=255)                                     gcaCCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCg                                                                                                 <  2:1051095/140‑1 (MQ=255)                                        ccTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATg                                                                                              >  1:323200/1‑140 (MQ=255)                                            aaCATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTg                                                                                          >  2:903023/1‑140 (MQ=255)                                                tatTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGc                                                                                      <  2:480618/140‑1 (MQ=255)                                                 atTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGcc                                                                                     >  2:1180393/1‑140 (MQ=255)                                                   tCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGa                                                                                                <  1:422411/127‑1 (MQ=255)                                                   tCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGa                                                                                                >  2:422411/1‑127 (MQ=255)                                                   tCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCAGCCAt                                                                                   >  2:463123/1‑140 (MQ=255)                                                      gTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACg                                                                                >  1:156801/1‑140 (MQ=255)                                                           aTAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGcga                                                                           >  2:1243027/1‑140 (MQ=255)                                                                   aaCGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATa                                                                                  <  2:658255/125‑1 (MQ=255)                                                                   aaCGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATa                                                                                  >  1:658255/1‑125 (MQ=255)                                                                   aaCGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGcc                                                                   <  1:1180393/140‑1 (MQ=255)                                                                    acgacgGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCg                                                                  <  1:1243027/140‑1 (MQ=255)                                                                    acgacgGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCg                                                                  <  1:63629/140‑1 (MQ=255)                                                                         gggAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGt                                                             <  2:596707/140‑1 (MQ=255)                                                                         gggAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGt                                                             >  1:527044/1‑140 (MQ=255)                                                                         gggAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGt                                                             >  2:39494/1‑140 (MQ=255)                                                                          ggAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTa                                                            <  1:39494/140‑1 (MQ=255)                                                                             gaAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACg                                                         <  2:850973/140‑1 (MQ=255)                                                                               aaTAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGtt                                                       <  2:1017357/140‑1 (MQ=255)                                                                                  aGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAg                                                    <  2:1009428/140‑1 (MQ=255)                                                                                   gTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAgg                                                   >  1:1164978/1‑140 (MQ=255)                                                                                                cagATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGc                                      <  2:156801/140‑1 (MQ=255)                                                                                                 agATGCCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCg                                     <  1:56277/140‑1 (MQ=255)                                                                                                   aTACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGc                                      <  2:970388/137‑1 (MQ=255)                                                                                                   aTACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGc                                      >  1:970388/1‑137 (MQ=255)                                                                                                             gCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATg                         >  1:304828/1‑140 (MQ=255)                                                                                                             gCTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATg                         >  2:537926/1‑140 (MQ=255)                                                                                                              cTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATgc                        >  2:1142254/1‑140 (MQ=255)                                                                                                              cTCATACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATgc                        >  1:298591/1‑140 (MQ=255)                                                                                                                 aTACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGcgcg                     >  2:526005/1‑140 (MQ=255)                                                                                                                 aTACTCGGGCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGcgcg                     >  1:769346/1‑140 (MQ=255)                                                                                                                       gggCTTTTGCCAGCCATCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCa               <  2:1130094/140‑1 (MQ=255)                                                                                                                        ggCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCAc              >  2:1016692/1‑140 (MQ=255)                                                                                                                        ggCTTTTGCCAGCCAGCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCAc              >  1:1166290/1‑140 (MQ=255)                                                                                                                                  agccagccACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGc    <  1:537926/140‑1 (MQ=255)                                                                                                                                   gccagccaCGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGcc   >  1:105196/1‑140 (MQ=255)                                                                                                                                    ccagCCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGccc  <  1:1016692/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        TAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCC  >  NC_000913/1073468‑1073738

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: ‑