Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,624,386	T→C	E522G (GAA→GGA)	yhiJ ←	YhiJ/YhiL family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,624,386	0	T	C	98.8%	277.6 / ‑5.4	81	E522G (GAA→GGA)	yhiJ	YhiJ/YhiL family protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (50/30);  total (51/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTTCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCA  >  NZ_CP009273/3624299‑3624472                                                                                        |                                                                                        gTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCa                                                                                       >  2:465490/1‑90 (MQ=255) gTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCa                                                                                       >  1:604698/1‑90 (MQ=255)   aaCATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑Tg                                                                                    >  1:1002063/1‑90 (MQ=255)   aaCATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑Tg                                                                                    >  2:1053810/1‑90 (MQ=255)    aCATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGc                                                                                   <  1:1039498/90‑1 (MQ=255)         aaTATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTa                                                                              <  1:53846/90‑1 (MQ=255)           tatTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAg                                                                            <  1:1098539/90‑1 (MQ=255)             ttCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGAc                                                                          <  1:1056122/90‑1 (MQ=255)              tCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTTCATTGCTGGGAAGAc                                                                          >  2:629774/1‑90 (MQ=255)               cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAg                                                                            >  2:2791/1‑86 (MQ=255)               cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAg                                                                            <  1:2791/86‑1 (MQ=255)               cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACtt                                                                        >  2:1048651/1‑90 (MQ=255)               cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACtt                                                                        >  2:464703/1‑90 (MQ=255)                     gtCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAggg                                                                  >  2:1070010/1‑90 (MQ=255)                       cTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAg                                                                >  2:578921/1‑90 (MQ=255)                       cTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAg                                                                <  1:398397/90‑1 (MQ=255)                           ggTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAGGGCATAAGTCGATTTCTGATTTTTCCA‑TGCCGGTAAGGCTTCAAGGGAGTGGt                                                            >  1:77913/1‑90 (MQ=255)                           ggTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGt                                                            <  2:756116/90‑1 (MQ=255)                            gTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAGGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGtt                                                           <  2:771561/90‑1 (MQ=255)                            gTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGtt                                                           <  1:739347/90‑1 (MQ=255)                             ttCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTa                                                          <  1:983920/90‑1 (MQ=255)                              tCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTAt                                                         >  2:54087/1‑90 (MQ=255)                                 gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑Tg                                                                                    >  1:785446/1‑60 (MQ=255)                                 gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑Tg                                                                                    <  2:785446/60‑1 (MQ=255)                                 gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTAtctc                                                      >  2:546586/1‑90 (MQ=255)                                 gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTAtctc                                                      >  1:69244/1‑90 (MQ=255)                                 gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTAtctc                                                      >  2:102058/1‑90 (MQ=255)                                 gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTAtctc                                                      >  1:360811/1‑90 (MQ=255)                                 gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTAtctc                                                      >  1:454926/1‑90 (MQ=255)                                     tGTTACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTTTTTTTTCCA‑ATCTGGTAAGGCTTCCAGGGAGTGGTTATCTCAGGc                                                  <  1:464703/90‑1 (MQ=255)                                     tGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGc                                                  <  1:1238978/90‑1 (MQ=255)                                      gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAg                                                                >  1:670348/1‑75 (MQ=255)                                      gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAg                                                                <  2:670348/75‑1 (MQ=255)                                      gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGtt                                                           >  2:89310/1‑80 (MQ=255)                                      gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGt                                                            <  1:410453/79‑1 (MQ=255)                                      gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGt                                                            >  2:410453/1‑79 (MQ=255)                                      gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTa                                                          <  1:1074068/81‑1 (MQ=255)                                      gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTa                                                          >  2:1074068/1‑81 (MQ=255)                                      gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGcc                                                 >  2:853433/1‑90 (MQ=255)                                      gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGcc                                                 >  1:544094/1‑90 (MQ=255)                                      gTAACATTGTGCCCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCCGGTAAGACTTCCAGGGAGTGGtt                                                           <  1:89310/80‑1 (MQ=255)                                       tAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGccc                                                <  2:1209206/90‑1 (MQ=255)                                           aTTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCa                                                                                       <  2:465476/48‑1 (MQ=255)                                           aTTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCa                                            >  2:905314/1‑90 (MQ=255)                                           aTTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCa                                            >  1:1069744/1‑90 (MQ=255)                                            ttGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAc                                           >  2:518078/1‑90 (MQ=255)                                            ttGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAc                                           <  1:1037350/90‑1 (MQ=255)                                            ttGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAc                                           <  2:174217/90‑1 (MQ=255)                                            ttGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAc                                           >  1:635019/1‑90 (MQ=255)                                             tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAggg                                                                  <  2:1219907/66‑1 (MQ=255)                                             tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAggg                                                                  >  1:1219907/1‑66 (MQ=255)                                             tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                          >  2:149389/1‑90 (MQ=255)                                             tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                          <  1:927842/90‑1 (MQ=255)                                             tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                          >  2:897516/1‑90 (MQ=255)                                              gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                         >  2:453340/1‑90 (MQ=255)                                              gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                         >  2:912420/1‑90 (MQ=255)                                              gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                         >  2:78407/1‑90 (MQ=255)                                              gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                         >  2:507050/1‑90 (MQ=255)                                              gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                         >  1:561633/1‑90 (MQ=255)                                                gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAggg                                                                  >  1:619751/1‑63 (MQ=255)                                                gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAggg                                                                  <  2:619751/63‑1 (MQ=255)                                                gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAgg                                       >  2:204140/1‑90 (MQ=255)                                                gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAgg                                       >  1:347835/1‑90 (MQ=255)                                                gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAgg                                       >  2:624356/1‑90 (MQ=255)                                                gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAgg                                       >  2:525763/1‑90 (MQ=255)                                                     ttCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTaa                                  >  2:1134850/1‑90 (MQ=255)                                                              aaGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAg                                                                <  2:978218/51‑1 (MQ=255)                                                              aaGGCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAg                                                                >  1:978218/1‑51 (MQ=255)                                                                 gCATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTat                      >  2:620050/1‑90 (MQ=255)                                                                  cATAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatg                     <  1:373048/90‑1 (MQ=255)                                                                   aTAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGc                                             <  1:499416/65‑1 (MQ=255)                                                                   aTAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGc                                             >  2:499416/1‑65 (MQ=255)                                                                   aTAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatgc                    >  2:768111/1‑90 (MQ=255)                                                                   aTAAGTCCATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatgc                    >  1:212333/1‑90 (MQ=255)                                                                          cATTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGa             <  1:26525/90‑1 (MQ=255)                                                                           aTTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatg                     <  2:1219038/81‑1 (MQ=255)                                                                           aTTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatg                     >  1:1219038/1‑81 (MQ=255)                                                                           aTTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGaa            >  1:666472/1‑90 (MQ=255)                                                                           aTTTCTGTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGaa            >  2:222901/1‑90 (MQ=255)                                                                                 gTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                          <  1:802983/54‑1 (MQ=255)                                                                                 gTTTTTTCCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                          >  2:802983/1‑54 (MQ=255)                                                                                  ttttttcCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACAttctt     <  1:507050/90‑1 (MQ=255)                                                                                     tttcCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCa  <  2:778047/90‑1 (MQ=255)                                                                                     tttcCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCa  <  1:825264/90‑1 (MQ=255)                                                                                     tttcCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCa  <  2:604369/90‑1 (MQ=255)                                                                                        |                                                                                        GTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTTCA‑TGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCA  >  NZ_CP009273/3624299‑3624472

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