Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,497,915	T→C	L30L (TTG→CTG)	cysG →	siroheme synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,497,915	0	T	C	100.0%	224.7 / NA	60	L30L (TTG→CTG)	cysG	siroheme synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/34);  total (26/34)

AATTTAATCAAACTGAAACGGCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGTTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGA  >  NC_000913/3497778‑3498039                                                                                                                                          |                                                                                                                             aaTTTAATCAAACTGAAACGGCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTg                                                                                                                            <  2:112741/140‑1 (MQ=255) aaTTTAATCAAACTGAAACGGCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTg                                                                                                                            >  2:182044/1‑140 (MQ=255)            aCTGAAACGGCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGc                                                                                                                 >  2:924009/1‑140 (MQ=255)                 aaCGGCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGgcg                                                                                                            >  2:245092/1‑140 (MQ=255)                  aCGGCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGgcgc                                                                                                           <  1:94041/140‑1 (MQ=255)                  aCGGCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGgcgc                                                                                                           <  2:561343/140‑1 (MQ=255)                   cGGCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGcgct                                                                                                          <  1:619498/140‑1 (MQ=255)                    ggCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGcgctc                                                                                                         <  2:269006/140‑1 (MQ=255)                     gCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGcgctcg                                                                                                        >  1:593209/1‑140 (MQ=255)                         ccGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTa                                                                                                    <  1:49299/140‑1 (MQ=255)                          cGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTaa                                                                                                   <  1:182044/140‑1 (MQ=255)                            gTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACg                                                                                                 >  2:1025381/1‑140 (MQ=255)                            gTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACg                                                                                                 <  2:136271/140‑1 (MQ=255)                              tAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGt                                                                                               <  1:76351/140‑1 (MQ=255)                              tAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGt                                                                                               <  1:547230/140‑1 (MQ=255)                                aTTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGa                                                                                             >  2:33667/1‑140 (MQ=255)                                 ttACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGTAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGaa                                                                                            <  1:328551/140‑1 (MQ=255)                                       aGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCAtt                                                                                      <  2:186057/140‑1 (MQ=255)                                         gggTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGaa                                                                                            >  1:773270/1‑132 (MQ=255)                                         gggTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGaa                                                                                            <  2:773270/132‑1 (MQ=255)                                              tttACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGttt                                                                               >  2:788494/1‑140 (MQ=255)                                               ttACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTa                                                                              <  1:924009/140‑1 (MQ=255)                                               ttACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTa                                                                              <  1:953955/140‑1 (MQ=255)                                                    tGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCca                                                                         >  2:643730/1‑140 (MQ=255)                                                    tGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCca                                                                         <  2:883939/140‑1 (MQ=255)                                                    tGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCca                                                                         <  1:1004643/140‑1 (MQ=255)                                                      gATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCcaca                                                                       >  1:163880/1‑140 (MQ=255)                                                               ccTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGc                                                              <  1:786966/140‑1 (MQ=255)                                                                cTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCa                                                             <  1:353577/140‑1 (MQ=255)                                                                 tataTTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCAt                                                            >  2:36314/1‑140 (MQ=255)                                                                   tatTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATgg                                                          >  2:447938/1‑140 (MQ=255)                                                                        tGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGa                                                     >  2:140869/1‑140 (MQ=255)                                                                        tGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGa                                                     >  1:925339/1‑140 (MQ=255)                                                                        tGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTC‑CGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGAt                                                    <  2:987547/140‑1 (MQ=255)                                                                         gCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGAt                                                    >  2:200197/1‑140 (MQ=255)                                                                          ccAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATg                                                   >  1:929509/1‑140 (MQ=255)                                                                          ccAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATg                                                   <  2:853594/140‑1 (MQ=255)                                                                          ccAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATg                                                   <  1:200197/140‑1 (MQ=255)                                                                             aTTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAg                                                >  2:165874/1‑140 (MQ=255)                                                                              ttACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAgg                                               <  2:476091/140‑1 (MQ=255)                                                                               tACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGc                                              >  1:489418/1‑140 (MQ=255)                                                                                   cgATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGt                                          <  1:643730/140‑1 (MQ=255)                                                                                               tCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCtcgtcg                              <  1:33667/140‑1 (MQ=255)                                                                                               tCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGGGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCtcgtcg                              >  2:332208/1‑140 (MQ=255)                                                                                                cTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGa                             >  1:1011299/1‑140 (MQ=255)                                                                                                cTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGa                             <  1:420419/140‑1 (MQ=255)                                                                                                      gTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGcc                       >  2:196546/1‑140 (MQ=255)                                                                                                         ggcggTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCAtt                    >  2:425049/1‑140 (MQ=255)                                                                                                         ggcggTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCAtt                    >  1:626454/1‑140 (MQ=255)                                                                                                         ggcggTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCAtt                    <  1:1014784/140‑1 (MQ=255)                                                                                                          gcggTGGTGATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCAttt                   <  1:1025381/140‑1 (MQ=255)                                                                                                               ggtgATGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTtgatga              <  2:307315/140‑1 (MQ=255)                                                                                                                   aTGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGc          <  1:738461/140‑1 (MQ=255)                                                                                                                    tGTCGCGGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGcc         <  1:36314/140‑1 (MQ=255)                                                                                                                       cgcgGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTtc      <  1:513386/140‑1 (MQ=255)                                                                                                                        gcgGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTtct     <  1:165874/140‑1 (MQ=255)                                                                                                                        gcgGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTtct     <  2:882479/140‑1 (MQ=255)                                                                                                                        gcgGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTtct     >  1:953379/1‑140 (MQ=255)                                                                                                                        gcgGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAAAGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTtct     >  2:874557/1‑140 (MQ=255)                                                                                                                         cgGAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTtctc    >  1:480454/1‑140 (MQ=255)                                                                                                                           gAACGCAAAGCAAGGCTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGa  >  1:48393/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                             AATTTAATCAAACTGAAACGGCTGCCGGTTAATTACTAAGGGGTTTTTACGTGGATCATTTGCCTATATTTTGCCAATTACGCGATCGCGACTGTCTGATTGTCGGCGGTGGTGATGTCGCGGAACGCAAAGCAAGGTTGCTGTTAGACGCAGGCGCTCGCTTAACGGTGAATGCATTAGCGTTTATTCCACAGTTCACCGCATGGGCAGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGA  >  NC_000913/3497778‑3498039

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