Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,337,024	G→A	S3L (TCG→TTG)	adiC ←	arginine:agmatine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,337,024	0	G	A	98.6%	235.9 / NA	71	S3L (TCG→TTG)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (39/31);  total (40/31)
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

GGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGG  >  NC_000913/4336929‑4337120                                                                                                |                                                                                                 ggTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGa                                                                                              <  2:921167/100‑1 (MQ=255) ggTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGa                                                                                              >  2:115081/1‑100 (MQ=255)  ttttGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGAc                                                                                             <  1:1459832/99‑1 (MQ=255)  gTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGAc                                                                                             <  1:1306274/100‑1 (MQ=255)   tttGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACa                                                                                            >  2:1459843/1‑100 (MQ=255)       cAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAAt                                                                                        <  1:311130/100‑1 (MQ=255)         ggTAACAGAAAAACACCTGACCGCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                      >  1:1053174/1‑100 (MQ=255)         ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                      >  2:1502211/1‑100 (MQ=255)          gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaa                                                                                     <  1:803688/100‑1 (MQ=255)            aaCGGAACAACACCTGACCCCAGAATATGCCCCGACACCATCATGGTGACGGGGATTAATCCCACTTTGTGATCATCAGCATCCAAAGACATAATTAAAc                                                                                   <  1:1805271/100‑1 (MQ=255)            aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAAc                                                                                   <  1:1594147/100‑1 (MQ=255)            aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAAc                                                                                   <  2:512314/100‑1 (MQ=255)             aCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAAc                                                                                   >  2:1091391/1‑99 (MQ=255)             aCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACt                                                                                  <  2:889644/100‑1 (MQ=255)              cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                                 >  2:1705436/1‑100 (MQ=255)              cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                                 >  2:1366444/1‑100 (MQ=255)               aGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTcc                                                                                >  2:185964/1‑100 (MQ=255)               aGAAAAACACCTGAACCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTcc                                                                                <  1:1536715/100‑1 (MQ=255)                     acacCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGaa                                                                          >  1:1643332/1‑100 (MQ=255)                       acCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAgg                                                                        <  2:1157106/100‑1 (MQ=255)                                 ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                              >  1:60420/1‑100 (MQ=255)                                  taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             <  2:1795495/100‑1 (MQ=255)                                  taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             <  2:1163333/100‑1 (MQ=255)                                   cttaatCCCCGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGa                                                            >  1:1168241/6‑100 (MQ=25)                                   cttaatCCCCGTCACCCTGATGGTGTCGGGGAATAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTACACTCCTGCGAAGGCGAGCTTCGTGa                                                            >  2:1362296/6‑100 (MQ=21)                                   cttaatCCCCGTCACCCTGATGGTGTCAGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGa                                                            >  1:563123/6‑100 (MQ=25)                                    atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                           <  2:1280995/100‑1 (MQ=255)                                      atTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACaa                                                         >  1:497447/1‑100 (MQ=255)                                      atTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCATGCGAAGGCGAGCTTCGTGACaa                                                         >  2:733206/1‑100 (MQ=255)                                        tCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATa                                                       >  2:116212/1‑100 (MQ=255)                                        tCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATa                                                       <  1:1009196/100‑1 (MQ=255)                                           ccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAt                                                    <  2:935929/100‑1 (MQ=255)                                                acCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGt                                               <  2:1819100/100‑1 (MQ=255)                                                 ccATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                              >  1:379367/1‑100 (MQ=255)                                                       gggTGACGGGGATTAAGCCCAATTTGAGAGAGTTAGCATAAGAGGAGATAATTGAGCTCTTGAGAAGGGGATGTTGGTGACAATAAATTACGtact                                            >  1:1542559/1‑93 (MQ=255)                                                         gTGACGGGGATTAAGCCCACTTTGTGCGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAGTTACGtcatatcata                                      >  2:1788473/1‑100 (MQ=255)                                                            aCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                   >  2:847106/1‑100 (MQ=255)                                                             cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAGAATAAATTACGTTATATCATACGcc                                  >  2:1295557/1‑100 (MQ=255)                                                             cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  >  2:1265814/1‑100 (MQ=255)                                                              ggggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  <  2:626221/99‑1 (MQ=255)                                                              ggggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCt                                 <  1:25931/100‑1 (MQ=255)                                                               gggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTAATATCATACACCt                                 >  2:871022/1‑99 (MQ=255)                                                                ggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGc                               >  2:11357/1‑100 (MQ=255)                                                                ggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGACGAGCTACGTGACAATAAATTACGTCATATCATACGCCTGc                               >  2:460371/1‑100 (MQ=255)                                                                 gATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACTATAAATTACGTCATATCATACGCCTg                                <  1:87106/98‑1 (MQ=255)                                                                  aTTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAt                             <  2:499276/100‑1 (MQ=255)                                                                   ttAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAtt                            >  2:1490828/1‑100 (MQ=255)                                                                    tAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                           >  2:925807/1‑100 (MQ=255)                                                                     aaGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCg                          >  2:499930/1‑100 (MQ=255)                                                                      aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                         >  2:1454722/1‑100 (MQ=255)                                                                       ncccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                        >  2:1041909/2‑100 (MQ=255)                                                                       gCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                        >  1:180244/1‑100 (MQ=255)                                                                            cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTACGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  1:115081/100‑1 (MQ=255)                                                                             tttGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  2:257302/99‑1 (MQ=255)                                                                                gtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt               <  2:1780730/100‑1 (MQ=255)                                                                                gtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt               >  1:1765599/1‑100 (MQ=255)                                                                                 tgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatta              <  1:1389609/100‑1 (MQ=255)                                                                                  gAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat             >  2:362565/1‑100 (MQ=255)                                                                                  gAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat             >  1:490446/1‑100 (MQ=255)                                                                                  gAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat             >  1:1471336/1‑100 (MQ=255)                                                                                   agcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattatt            <  1:1362286/100‑1 (MQ=255)                                                                                   agcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattatt            >  2:995527/1‑100 (MQ=255)                                                                                    gcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAAttatta           >  1:410691/1‑100 (MQ=255)                                                                                     catcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGATCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATtattat          <  2:511948/100‑1 (MQ=255)                                                                                      atcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATc         <  2:180244/100‑1 (MQ=255)                                                                                      atcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATc         <  1:1139669/100‑1 (MQ=255)                                                                                        cagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGa       <  1:1654217/100‑1 (MQ=255)                                                                                          gcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGa       >  2:30771/1‑98 (MQ=255)                                                                                            atcCAAAGACATAATTAATCTCCTGCGAAGGCGAGCTTCGGGACAATAAATTACGTCATATCTTACGCCTGCATTCGATAGGTGAATTATTATCGACCgg   <  1:1765737/100‑1 (MQ=255)                                                                                             tcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCggg  >  2:837744/1‑100 (MQ=255)                                                                                             tcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGACTAAATTAGATAGTTGAATTATTATCGACCggg  >  2:151586/1‑100 (MQ=255)                                                                                                |                                                                                                 GGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGG  >  NC_000913/4336929‑4337120

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