Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	913,007	T→C	T270A (ACC→GCC)	hcp ←	protein S‑nitrosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	913,007	0	T	C	100.0%	159.6 / NA	44	T270A (ACC→GCC)	hcp	protein S‑nitrosylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/25);  total (19/25)

CGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGACCGT  >  NC_000913/912870‑913138                                                                                                                                          |                                                                                                                                    cGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGgcgt                                                                                                                                   >  1:1017911/1‑140 (MQ=255)   aGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAg                                                                                                                                 >  1:95124/1‑140 (MQ=255)   aGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAg                                                                                                                                 <  2:589578/140‑1 (MQ=255)         tcaCGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTa                                                                                                                           >  1:827770/1‑140 (MQ=255)         tcaCGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTa                                                                                                                           >  2:38525/1‑140 (MQ=255)                ggggCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCg                                                                                                                    >  2:306743/1‑140 (MQ=255)                ggggCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCg                                                                                                                    <  1:306743/140‑1 (MQ=255)                     cAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCt                                                                                                               <  1:619480/140‑1 (MQ=255)                       gggAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTc                                                                                                             >  2:631773/1‑140 (MQ=255)                        ggAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCg                                                                                                            <  1:574009/140‑1 (MQ=255)                             cGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTg                                                                                                       >  1:1089818/1‑140 (MQ=255)                             cGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTg                                                                                                       >  1:755740/1‑140 (MQ=255)                                  gAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGTCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCca                                                                                                  >  2:1102314/1‑140 (MQ=255)                                  gAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCca                                                                                                  <  2:500959/140‑1 (MQ=255)                                  gAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCca                                                                                                  <  2:95124/140‑1 (MQ=255)                                  gAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCca                                                                                                  >  2:1073457/1‑140 (MQ=255)                                  gAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCca                                                                                                  <  2:903675/140‑1 (MQ=255)                                    aCTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCcagc                                                                                                <  1:396277/140‑1 (MQ=255)                                    aCTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCcagc                                                                                                <  2:1221473/140‑1 (MQ=255)                                    aCTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTAGGTCTGTTCcagc                                                                                                <  1:719450/140‑1 (MQ=255)                                     cTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCcagca                                                                                               <  1:1102314/140‑1 (MQ=255)                                        cACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGt                                                                                            >  2:1195236/1‑140 (MQ=255)                                              cTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTaga                                                                                      >  2:1110167/1‑140 (MQ=255)                                               tGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTagag                                                                                     <  1:631773/140‑1 (MQ=255)                                               tGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTagag                                                                                     <  2:755740/140‑1 (MQ=255)                                                gATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAgaga                                                                                    >  2:70267/1‑140 (MQ=255)                                                                      gTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGaaa                                                              <  1:1000232/140‑1 (MQ=255)                                                                        tACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATc                                                            >  2:299756/1‑140 (MQ=255)                                                                         aCCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCa                                                           >  1:843191/1‑140 (MQ=255)                                                                           cGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGa                                                         <  1:70267/140‑1 (MQ=255)                                                                           cGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGa                                                         <  2:1089818/140‑1 (MQ=255)                                                                            gACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGTAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGaa                                                        <  1:1110167/140‑1 (MQ=255)                                                                              ccAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGaa                                                               >  2:512318/1‑131 (MQ=255)                                                                              ccAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGaa                                                               <  1:512318/131‑1 (MQ=255)                                                                                 gATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGaa                                                        >  2:1020420/1‑135 (MQ=255)                                                                                 gATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGaa                                                        <  1:1020420/135‑1 (MQ=255)                                                                                    gCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTa                                                <  1:1195236/140‑1 (MQ=255)                                                                                            ttACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGt                                        <  1:299756/140‑1 (MQ=255)                                                                                                  aGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCttt                                  >  2:640620/1‑140 (MQ=255)                                                                                                      ccGGGGAGCCATGCGCAGGCACCATTTCGCCGTGGGCGTAAACGTTAACGCCCGTGCTTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGAAAATCAGAATGCTTTTACCCGCCGTCGCTTTGACg                              <  2:657110/140‑1 (MQ=255)                                                                                                                   cgcAGGCAGCATTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTg                 >  2:1147128/1‑140 (MQ=255)                                                                                                                         cagcaTTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCg           <  1:1045278/140‑1 (MQ=255)                                                                                                                           gcaTTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCggg         <  1:334074/140‑1 (MQ=255)                                                                                                                             aTTTCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTg       <  1:1073457/140‑1 (MQ=255)                                                                                                                              tttCGCCGTGGGCGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGa      <  2:827770/140‑1 (MQ=255)                                                                                                                                  gCCGTGGGCGTAGACATTAACGCCCCTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGACCGt  >  2:226287/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                    CGAGGTCATCACGATGGGGCCAGGGAAACGAGCGAACTCCACTTGCTGATTCTGCCAGCCGCTGCCGTAGTTACCGACCAGATGCTTGAATTTACGCAGCTCCGGGTAGCCATGCGCAGGCAGCATTTCGCCGTGGGTGTAGACATTAACGCCCGTGCCTTCGGTCTGTTCCAGCAGGTTGTAGAGATCTTTGAGATCGTGACCGGAAATCAGAATGCATTTACCCGCCGTCGCTTTGACGTTGACCTGGGTTGGCGTCGGGTGACCGT  >  NC_000913/912870‑913138

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