Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,494,539	T→C	F222S (TTC→TCC)	nanX →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,494,539	0	T	C	97.3%	235.2 / ‑3.1	75	F222S (TTC→TCC)	nanX	MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (33/40);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAG  >  NZ_CP009273/4494452‑4494625                                                                                        |                                                                                       aGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTccc                                                                                      >  2:704810/1‑90 (MQ=255) aGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTccc                                                                                      <  2:619588/90‑1 (MQ=255) aGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTccc                                                                                      <  1:35086/90‑1 (MQ=255)     gTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGt                                                                                  >  1:323762/1‑90 (MQ=255)           tGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGttt                                                                            <  2:1552748/90‑1 (MQ=255)            gAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGtttt                                                                           >  1:1469848/1‑90 (MQ=255)                aTAAATGTAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTc                                                                       <  1:570121/90‑1 (MQ=255)                 tAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCt                                                                      >  1:29467/1‑90 (MQ=255)                   aaTATAAGGCTAAATGAACATTTTTGTCTGTCTTCAGAAAACCAAATCTTGCAATCTCTATGATCGTTTCCCTGTTCTGTTTTTGTCTAt                                                                    >  1:732484/1‑90 (MQ=255)                    atatAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTAtt                                                                   <  2:559546/90‑1 (MQ=255)                     tataAGGATAAATCAACATTTTTGTCTGTCTTCAGAGAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTAttt                                                                  <  1:1310616/90‑1 (MQ=255)                         aGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTg                                                              >  2:416695/1‑90 (MQ=255)                           gATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCa                                                            >  2:17353/1‑90 (MQ=255)                            aTAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCaa                                                           >  1:1393251/1‑90 (MQ=255)                             tAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCaaa                                                          >  2:650620/1‑90 (MQ=255)                              aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGttttt                                                                          >  1:1452315/1‑73 (MQ=255)                              aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGttttt                                                                          <  2:1452315/73‑1 (MQ=255)                              aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                         <  1:1039169/90‑1 (MQ=255)                              aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                         >  2:1452075/1‑90 (MQ=255)                              aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                         >  2:1005286/1‑90 (MQ=255)                                    aCATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCg                                                   >  2:225165/1‑90 (MQ=255)                                    aCATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCg                                                   >  1:871670/1‑90 (MQ=255)                                    aCATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCg                                                   >  1:647973/1‑90 (MQ=255)                                         tttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATaaa                                              >  1:1535015/1‑90 (MQ=255)                                         tttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATaaa                                              >  2:513027/1‑90 (MQ=255)                                          ttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAAc                                             <  2:1560440/90‑1 (MQ=255)                                              ctgtctTCAGAAACCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         <  1:1291622/90‑1 (MQ=255)                                              ctgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         <  1:458730/90‑1 (MQ=255)                                              ctgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         <  2:1013665/90‑1 (MQ=255)                                               tgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGact                                        >  2:271609/1‑90 (MQ=255)                                                  ctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGactact                                     <  1:158230/90‑1 (MQ=255)                                                   ttCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActt                                    <  1:768161/90‑1 (MQ=255)                                                   ttCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActt                                    <  2:1022622/90‑1 (MQ=255)                                                   ttCAGAAAACCACATCTTTCAATCTCTAAGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActt                                    <  1:704810/90‑1 (MQ=255)                                                       gAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActtcctt                                <  2:1409075/90‑1 (MQ=255)                                                       gAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActtcctt                                <  2:389659/90‑1 (MQ=255)                                                              acaTCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCt                         >  1:1094028/1‑90 (MQ=255)                                                                aTCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTgg                       <  2:871670/90‑1 (MQ=255)                                                                  cTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCa                     <  1:246617/90‑1 (MQ=255)                                                                   tttCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAg                    >  2:128725/1‑90 (MQ=255)                                                                      cAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATa                 >  2:119184/1‑90 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGGCTACTTCCTTCCTACCTGGCAGATAATg              <  2:1099230/90‑1 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  2:1094028/90‑1 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  2:922563/90‑1 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:984864/90‑1 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:128725/90‑1 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:1439377/90‑1 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:1452075/90‑1 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:1551310/90‑1 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:162034/90‑1 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:78254/90‑1 (MQ=255)                                                                         tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:873382/90‑1 (MQ=255)                                                                           tctATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  2:678543/90‑1 (MQ=255)                                                                           tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         >  2:1258513/1‑61 (MQ=255)                                                                           tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         <  1:1258513/61‑1 (MQ=255)                                                                           tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  1:1145631/90‑1 (MQ=255)                                                                           tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  2:851440/90‑1 (MQ=255)                                                                           tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  2:1250067/90‑1 (MQ=255)                                                                           tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  2:296232/90‑1 (MQ=255)                                                                            ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  1:1375991/1‑90 (MQ=255)                                                                            ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  2:1408662/1‑90 (MQ=255)                                                                            ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  1:596260/1‑90 (MQ=255)                                                                            ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  2:1065878/1‑90 (MQ=255)                                                                            ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  1:833982/1‑90 (MQ=255)                                                                            ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCCGGCAGATAATGGAg           >  1:1030930/1‑90 (MQ=255)                                                                              aTGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGtt         <  1:1161139/90‑1 (MQ=255)                                                                              aTGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGtt         >  2:618734/1‑90 (MQ=255)                                                                                   cGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCCGATAATGGAGTTAATAc    <  1:1576399/90‑1 (MQ=255)                                                                                   cGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAc    <  2:1178305/90‑1 (MQ=255)                                                                                    gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTa                            <  2:28421/65‑1 (MQ=255)                                                                                    gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTa                            >  1:28421/1‑65 (MQ=255)                                                                                    gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              >  2:468043/1‑79 (MQ=255)                                                                                    gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:1290541/79‑1 (MQ=255)                                                                                    gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:468043/79‑1 (MQ=255)                                                                                    gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              >  2:1290541/1‑79 (MQ=255)                                                                                    gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa   >  2:429824/1‑90 (MQ=255)                                                                                    gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa   >  1:1485912/1‑90 (MQ=255)                                                                                    gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa   >  2:217319/1‑90 (MQ=255)                                                                                     tttcCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAg  <  1:8079/90‑1 (MQ=255)                                                                                        |                                                                                       AGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAG  >  NZ_CP009273/4494452‑4494625

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