Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,013,706	T→C	L150L (TTG→CTG)	pqiB →	paraquat‑inducible, SoxRS‑regulated MCE domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,013,706	0	T	C	100.0%	165.5 / NA	48	L150L (TTG→CTG)	pqiB	paraquat‑inducible, SoxRS‑regulated MCE domain protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/25);  total (23/25)

CGTCTTTTGGGTGGTGAAACCGCAGATTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGTTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTAGGTTCGGTTGAAACCAGCACCTTCGATACACAAAAA  >  NC_000913/1013564‑1013843                                                                                                                                               |                                                                                                                                          cGTCTTTTGGGTGGTGAAACCGCAGATTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCg                                                                                                                                   >  2:204843/1‑151 (MQ=255) cGTCTTTTGGGTGGTGAAACCGCAGATTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCg                                                                                                                                   <  1:457257/151‑1 (MQ=255)    cTTTTGGGTGGTGAAACCGCAGATTGGTCGCGGAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCt                                                                                                                                >  1:237661/1‑151 (MQ=255)        tGGGTGGTGAAACCGCAGATTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATg                                                                                                                            >  1:599253/1‑151 (MQ=255)         gggTGGTGAAACCGCAGATTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATg                                                                                                                            <  2:282121/150‑1 (MQ=255)            tggtgAAACCGCAGGGTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGATGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATGAATACGATTTGCTGGACTAGCCACCGCTGGCCCCGCCTGATGCg                                                                                                                          >  2:253602/1‑149 (MQ=255)            tggtgAAACCGCAGATTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGa                                                                                                                         >  2:87655/1‑150 (MQ=255)             ggtgAAACCGCAGATTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGaa                                                                                                                        >  1:57092/1‑150 (MQ=255)                 aaaCCGCAGATTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGt                                                                                                                    >  1:32149/1‑150 (MQ=255)                     cGCAGATTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCg                                                                                                               <  1:329122/151‑1 (MQ=255)                           ttGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGACAGCAAAATGGATAACTACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGAt                                                                                                         >  2:391788/1‑151 (MQ=255)                            tGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGAt                                                                                                         <  1:47492/150‑1 (MQ=255)                                  cgAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGAt                                                                                                  <  2:496904/151‑1 (MQ=255)                                   gAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATa                                                                                                 <  1:84558/151‑1 (MQ=255)                                     tgggATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGAATCCGTGTGATTCTCGATAg                                                                                                >  2:291510/2‑150 (MQ=255)                                      gggATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGc                                                                                               >  1:390444/1‑150 (MQ=255)                                         aTTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCaaaa                                                                                           <  2:178994/151‑1 (MQ=255)                                         aTTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCaaaa                                                                                           <  2:5316/151‑1 (MQ=255)                                          ttAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCaaaa                                                                                           >  2:145001/1‑150 (MQ=255)                                              cGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGcc                                                                                      <  2:14765/151‑1 (MQ=255)                                                    gggAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGc                                                                                  <  2:122732/149‑1 (MQ=255)                                                      gAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGct                                                                              >  2:597599/1‑151 (MQ=255)                                                                   tGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAgg                                                                     >  1:471253/1‑147 (MQ=255)                                                                      aGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCg                                                              >  2:397559/1‑151 (MQ=255)                                                                             atCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGt                                                       >  2:59614/1‑151 (MQ=255)                                                                               cGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTc                                                     <  2:59612/151‑1 (MQ=255)                                                                                     gCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTg                                                 <  1:407018/149‑1 (MQ=255)                                                                                     gCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGc                                               >  1:276858/1‑151 (MQ=255)                                                                                      cAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCt                                              <  2:27499/151‑1 (MQ=255)                                                                                          cAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCg                                          <  2:660494/151‑1 (MQ=255)                                                                                           aGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCgg                                         >  2:427860/1‑151 (MQ=255)                                                                                                   aaGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTc                                 <  2:631179/151‑1 (MQ=255)                                                                                                    aGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCg                                <  2:426502/151‑1 (MQ=255)                                                                                                          caAAATGGATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGaa                          <  2:230231/151‑1 (MQ=255)                                                                                                                ggATAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAg                     <  1:673972/150‑1 (MQ=255)                                                                                                                ggATAAATACGATTTGCTGGACTCGCAACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGc                    <  1:284017/151‑1 (MQ=255)                                                                                                                  aTAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCAc                  >  1:158748/1‑151 (MQ=255)                                                                                                                   tAAATACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCAcc                 >  2:83815/1‑151 (MQ=255)                                                                                                                    aaaTACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCACCt                <  1:2656/151‑1 (MQ=255)                                                                                                                    aaaTACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCACCt                <  1:99007/151‑1 (MQ=255)                                                                                                                      aTACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCACCTTc              <  1:511428/151‑1 (MQ=255)                                                                                                                      aTACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCACCTTc              <  1:135967/151‑1 (MQ=255)                                                                                                                      aTACGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCACCTTc              >  1:66729/1‑151 (MQ=255)                                                                                                                         cGATTTGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCACCTTCGAt           >  1:671209/1‑151 (MQ=255)                                                                                                                            tttGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCACCTTCGATac         >  2:313479/1‑150 (MQ=255)                                                                                                                             ttGCTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCACCTTCGATaca        <  1:19486/150‑1 (MQ=255)                                                                                                                                cTGGACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCACCTTCGATACACaa     >  2:59336/1‑150 (MQ=255)                                                                                                                                  ggACTCGCCACCGCTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTGGGTTCGGTTGAAACCAGCACCTTCGATACACaaaaa  <  1:548803/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                          CGTCTTTTGGGTGGTGAAACCGCAGATTGGTCGCGAAGGGATTAGCGGCCTGGGAACGCTGCTGTCTGGAGTTTATATCGAACTGCAGCCAGGCGCGAAAGGCAGCAAAATGGATAAATACGATTTGCTGGACTCGCCACCGTTGGCCCCGCCTGATGCGAAAGGTATCCGTGTGATTCTCGATAGCAAAAAAGCCGGGCAGCTCTCGCCAGGAGATCCGGTGCTGTTCCGTGGCTATCGGGTAGGTTCGGTTGAAACCAGCACCTTCGATACACAAAAA  >  NC_000913/1013564‑1013843

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