Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	886,189	A→G	G7G (GGT→GGC)	grxA ←	glutaredoxin 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	886,189	0	A	G	100.0%	189.2 / NA	57	G7G (GGT→GGC)	grxA	glutaredoxin 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/21);  total (36/21)

ATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGT  >  NZ_CP009273/886105‑886274                                                                                     |                                                                                      atactgatactgAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCaaa                                                                                  <  1:447670/90‑1 (MQ=255)  tactgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCaaaa                                                                                 >  2:1205568/1‑90 (MQ=255)    ctgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAAt                                                                               >  1:349671/1‑90 (MQ=255)     tgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATa                                                                              >  2:1271143/1‑90 (MQ=255)      gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCGGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  2:456031/1‑90 (MQ=255)      gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  1:519014/1‑90 (MQ=255)      gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  1:1241321/1‑90 (MQ=255)      gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  1:1387185/1‑90 (MQ=255)      gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  1:1539583/1‑90 (MQ=255)      gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  1:921411/1‑90 (MQ=255)      gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  1:701555/1‑90 (MQ=255)      gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             <  1:242922/90‑1 (MQ=255)       atactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAAc                                                                            >  2:848120/1‑90 (MQ=255)       atactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAAc                                                                            >  1:467057/1‑90 (MQ=255)         actgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACgg                                                                          >  2:1107493/1‑90 (MQ=255)           tgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGtt                                                                        <  2:216872/90‑1 (MQ=255)            gaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGttt                                                                       <  2:519014/90‑1 (MQ=255)                atcatcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCat                                                                   <  1:990334/90‑1 (MQ=255)                atcatcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCat                                                                   <  1:1271143/90‑1 (MQ=255)                atcatcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTATGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCat                                                                   <  2:790107/90‑1 (MQ=255)                   atcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCattat                                                                >  2:1513786/1‑90 (MQ=255)                      gcgTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTtc                                                             >  2:1061264/1‑90 (MQ=255)                        gTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                           <  1:812783/90‑1 (MQ=255)                        gTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                           >  1:766437/1‑90 (MQ=255)                               gCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTGTTTCTCTCCTCAt                                                    <  1:427613/90‑1 (MQ=255)                                cTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATa                                                   <  2:596226/90‑1 (MQ=255)                                     tttCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGAttt                                              <  1:80925/90‑1 (MQ=255)                                            gCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTg                                       >  2:600401/1‑90 (MQ=255)                                                     tttGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTACATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCa                              <  1:631007/90‑1 (MQ=255)                                                      ttGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCAc                             >  2:664295/1‑90 (MQ=255)                                                        gcacgcacACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGc                           <  2:467057/90‑1 (MQ=255)                                                        gcacgcacACAGTAAGGGCAACCCGAACGGCCAAAAATAACGG‑TTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCt                          <  2:1125796/90‑1 (MQ=255)                                                          acgcacACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTa                         >  1:1241166/1‑90 (MQ=255)                                                                   gTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTAtt                <  1:600401/90‑1 (MQ=255)                                                                       gggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCt            >  1:1103603/1‑90 (MQ=255)                                                                       gggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCt            >  2:1153953/1‑90 (MQ=255)                                                                       gggCAACCCGAACGGCCAAAAATAACGGCTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCt            >  2:724047/1‑90 (MQ=255)                                                                         gCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                         >  1:456353/1‑59 (MQ=255)                                                                         gCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                         <  2:456353/59‑1 (MQ=255)                                                                         gCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGa          >  1:121395/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:1410539/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:698936/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:771891/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:1560691/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:105364/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:1472644/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:1355031/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:358573/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:341255/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:225170/1‑90 (MQ=255)                                                                           aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:1147025/1‑90 (MQ=255)                                                                            aCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAg       <  1:1153953/90‑1 (MQ=255)                                                                             cccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGt      <  2:1241321/90‑1 (MQ=255)                                                                             cccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGt      <  2:1472644/90‑1 (MQ=255)                                                                             cccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGt      <  1:1255489/90‑1 (MQ=255)                                                                              ccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTa     <  2:1341448/90‑1 (MQ=255)                                                                                 aaCGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGt  >  1:191951/1‑90 (MQ=255)                                                                                     |                                                                                      ATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGT  >  NZ_CP009273/886105‑886274

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