Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,961,849	C→T	R391R (CGC→CGT)	rep →	ATP‑dependent DNA helicase Rep

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,961,849	0	C	T	100.0%	172.1 / NA	49	R391R (CGC→CGT)	rep	ATP‑dependent DNA helicase Rep
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/23);  total (26/23)

TCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCA  >  NC_000913/3961717‑3961970                                                                                                                                     |                                                                                                                          tCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGa                                                                                                                    >  2:306107/1‑140 (MQ=255)   tttATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACt                                                                                                                  <  2:453692/140‑1 (MQ=255)     tATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTaa                                                                                                                >  1:354178/1‑140 (MQ=255)       tCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAAcc                                                                                                              <  2:45412/140‑1 (MQ=255)          cgGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCgg                                                                                                           >  1:140948/1‑140 (MQ=255)          cgGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCgg                                                                                                           >  1:139622/1‑140 (MQ=255)          cgGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCgg                                                                                                           >  1:255147/1‑140 (MQ=255)          cgGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCgg                                                                                                           >  1:367370/1‑140 (MQ=255)            gTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAAcc                                                                                                              <  2:320915/135‑1 (MQ=255)            gTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAAcc                                                                                                              >  1:320915/1‑135 (MQ=255)            gTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGAc                                                                                                         >  1:197467/1‑140 (MQ=255)             tAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGGCTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACg                                                                                                        <  2:533224/140‑1 (MQ=255)             tAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACg                                                                                                        >  2:22590/1‑140 (MQ=255)               aCCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGAt                                                                                                      >  1:449336/1‑140 (MQ=255)               aCCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGAt                                                                                                      <  1:904648/140‑1 (MQ=255)                  aTCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGCCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGAc                                                                                                   >  1:549274/1‑140 (MQ=255)                     aGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAgc                                                                                                >  1:168518/1‑140 (MQ=255)                         gcgGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCAt                                                                                            >  1:779408/1‑140 (MQ=255)                          cgGGTGTTTGAAAAATTTCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAgcg                                                                                               <  2:882817/136‑1 (MQ=255)                          cgGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAgcg                                                                                               >  1:882817/1‑136 (MQ=255)                             gtgtTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCt                                                                                        >  2:139979/1‑140 (MQ=255)                               gtTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGc                                                                                      >  1:340356/1‑140 (MQ=255)                                tttGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCg                                                                                     >  1:872409/1‑140 (MQ=255)                                tttGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCg                                                                                     <  2:139622/140‑1 (MQ=255)                                 ttGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGt                                                                                    <  2:140948/140‑1 (MQ=255)                                 ttGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGt                                                                                    <  2:197467/140‑1 (MQ=255)                                 ttGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGt                                                                                    <  2:449336/140‑1 (MQ=255)                                  tGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTa                                                                                   <  2:872409/140‑1 (MQ=255)                                                     aaCCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCg                                                                <  1:139979/140‑1 (MQ=255)                                                      aCCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCga                                                               <  1:806775/140‑1 (MQ=255)                                                                 tACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCg                                                                <  1:191423/128‑1 (MQ=255)                                                                 tACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCg                                                                >  2:191423/1‑128 (MQ=255)                                                                 tACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGccc                                                    >  2:934110/1‑140 (MQ=255)                                                                    aaaaTATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGc                                                 <  2:558454/140‑1 (MQ=255)                                                                        tatCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACg                                             >  2:409288/1‑140 (MQ=255)                                                                        tatCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACg                                             >  1:359617/1‑140 (MQ=255)                                                                             ggtggtACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGaa                                        <  1:145214/140‑1 (MQ=255)                                                                                    cgtcgtTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAGCACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTg                                 >  2:901374/1‑140 (MQ=255)                                                                                       cgtTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGt                              <  2:103573/140‑1 (MQ=255)                                                                                       cgtTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGt                              <  2:354178/140‑1 (MQ=255)                                                                                         tttttCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGa                            <  2:168518/140‑1 (MQ=255)                                                                                         tttttCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGa                            >  1:582862/1‑140 (MQ=255)                                                                                                  cGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGAt                   >  2:485278/1‑140 (MQ=255)                                                                                                    tCCTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGa                 >  1:890940/1‑140 (MQ=255)                                                                                                      cTGAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg               <  1:485278/140‑1 (MQ=255)                                                                                                        gAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgcg             <  2:150906/140‑1 (MQ=255)                                                                                                        gAAATCAAGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgcg             <  2:582862/140‑1 (MQ=255)                                                                                                               aGGACTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATaaa      <  1:934110/140‑1 (MQ=255)                                                                                                                   cTTGCTGGCTTATCTGCGTGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCa  >  1:789534/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                          TCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCA  >  NC_000913/3961717‑3961970

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