Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,684	+A	coding (400/501 nt)	ytfI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,681	1	.	A	98.8%	262.7 / ‑3.2	82	coding (397/501 nt)	ytfI	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (39/42);  total (39/43)
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.85e-01

TACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑A‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGT  >  NC_000913/4435606‑4435751                                                                             |                                                                        tACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac                                                                          >  1:2293702/1‑76 (MQ=255) tACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac                                                                          <  1:2084368/76‑1 (MQ=255)  aCGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a                                                                        >  1:5058080/1‑76 (MQ=255)  aCGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a                                                                        >  2:2648186/1‑76 (MQ=255)  aCGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a                                                                        >  1:2267685/1‑76 (MQ=255)    gaagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATCTCCGTTAGaac‑a‑aa                                                                     >  1:3375308/1‑76 (MQ=255)     aagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aaa                                                                    <  2:4447842/76‑2 (MQ=255)      agaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aaac                                                                   <  2:2811962/76‑3 (MQ=255)      agaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aaac                                                                   <  2:315392/76‑3 (MQ=255)       gaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aaaca                                                                  <  2:4760026/76‑4 (MQ=255)       gaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aaaca                                                                  <  2:3356476/76‑4 (MQ=255)          aTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGAc                                                                >  1:4589162/1‑76 (MQ=255)           tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACg                                                               >  2:720202/1‑76 (MQ=255)           tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACg                                                               <  2:3674838/76‑1 (MQ=255)            aCCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGa                                                              <  1:4792171/76‑1 (MQ=255)            aCCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGa                                                              >  1:1012038/1‑76 (MQ=255)            aCCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGa                                                              >  2:898365/1‑76 (MQ=255)            aCCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGa                                                              >  1:1062157/1‑76 (MQ=255)              ccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGAtt                                                            <  1:795408/76‑1 (MQ=255)               cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTg                                                           <  2:2875092/76‑1 (MQ=255)                aTAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGt                                                          <  1:1620103/76‑1 (MQ=255)                 tAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTa                                                         <  1:4969962/76‑1 (MQ=255)                  aaCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTAt                                                        <  1:2485898/76‑1 (MQ=255)                  aaCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTAt                                                        <  1:2444779/76‑1 (MQ=255)                     aaTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGcc                                                     >  1:4712084/1‑76 (MQ=255)                       tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTg                                                   <  1:1052811/76‑1 (MQ=255)                        aaaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTgg                                                  <  1:4808181/76‑1 (MQ=255)                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGa                                                <  2:2950197/76‑1 (MQ=255)                           ggAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGat                                               <  2:2494247/75‑1 (MQ=255)                            gAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGata                                              >  1:1971318/1‑76 (MQ=255)                            gAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGata                                              <  1:4731134/76‑1 (MQ=255)                            gAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGata                                              >  1:1243631/1‑76 (MQ=255)                            gAATATAACTCATGCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGAATGTATGCCTGGGata                                              >  1:4121613/1‑76 (MQ=255)                              atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATAtt                                            <  2:4616447/76‑1 (MQ=255)                                ataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTTTATGCCTGGGATATTCt                                          >  2:2662010/1‑76 (MQ=255)                                 taACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTa                                         <  2:5128913/76‑1 (MQ=255)                                  aaCTCATTCATGGCTTTCAAAGACCGTGAATTTCCGTAAGATCAA‑AACCGGCTATTGTCTGCCTAGGATATTCTAc                                        <  1:1556011/76‑1 (MQ=255)                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTAc                                        >  2:3560307/1‑76 (MQ=255)                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTAc                                        <  2:1062157/76‑1 (MQ=255)                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTAc                                        >  2:3598829/1‑76 (MQ=255)                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTAc                                        <  2:1012038/76‑1 (MQ=255)                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTAc                                        >  1:1800349/1‑76 (MQ=255)                                     tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGa                                     <  1:41251/76‑1 (MQ=255)                                     tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGACTGGGATATTCTCCAGa                                     >  2:209578/1‑76 (MQ=255)                                      cattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaa                                    <  2:5058079/76‑1 (MQ=255)                                      cattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaa                                    >  1:3357038/1‑76 (MQ=255)                                      cattcatGGATTTCAAAAACCGTGAATCTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaa                                    <  2:3375308/76‑1 (MQ=255)                                       attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTa                                                         >  1:1262110/1‑54 (MQ=255)                                       attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTa                                                         <  2:1262110/54‑1 (MQ=255)                                       attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaaa                                   >  1:2746279/1‑76 (MQ=255)                                       attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaaa                                   <  1:4092893/76‑1 (MQ=255)                                         tcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaaaaa                                 <  2:755380/76‑1 (MQ=255)                                         tcatGGATTTCAAAAACCGTGAATTTCCGCTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaaaaa                                 <  1:3232948/76‑1 (MQ=255)                                          catGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaaaaaa                                >  1:5155951/1‑76 (MQ=255)                                          catGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGaaaaaa                                <  1:2642957/76‑1 (MQ=255)                                          catGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGCATATTCTACAGaaaaaa                                >  2:5537559/1‑76 (MQ=255)                                           atGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAAc                               <  1:475130/76‑1 (MQ=255)                                            tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACa                              >  2:1097188/1‑76 (MQ=255)                                             ggATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTTTATGCCTGGGATATTCTACAGAAAAAACAg                             <  1:2662010/76‑1 (MQ=255)                                              gATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGt                            >  2:640526/1‑76 (MQ=255)                                               aTTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                           >  2:1326218/1‑76 (MQ=255)                                                tttCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                          >  1:928759/1‑76 (MQ=255)                                                 ttCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTg                         >  1:3630314/1‑76 (MQ=255)                                                 ttCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTg                         >  1:3513524/1‑76 (MQ=255)                                                  tCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGa                        >  2:4126304/1‑76 (MQ=255)                                                  tCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGa                        >  1:1815967/1‑76 (MQ=255)                                                   cAAAAACCGTGAATTTCCGTTAGATCAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGAc                       <  2:1619322/76‑1 (MQ=255)                                                   cAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGAc                       <  1:1969589/76‑1 (MQ=255)                                                   cAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGAc                       >  1:4048162/1‑76 (MQ=255)                                                     aaaaCCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                     >  2:3457086/1‑76 (MQ=255)                                                     aaaaCCGTGAAGTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                     >  2:3287642/1‑76 (MQ=255)                                                        aCCGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAcc                  <  1:3545121/76‑1 (MQ=255)                                                         ccGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                 <  2:1500561/76‑1 (MQ=255)                                                         ccGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                 <  2:1545834/76‑1 (MQ=255)                                                         ccGTGAATTGCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                 <  2:2423038/76‑1 (MQ=255)                                                          cGTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGa                >  2:549455/1‑76 (MQ=255)                                                           gTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa               <  2:3272433/76‑1 (MQ=255)                                                           gTGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa               >  1:4068451/1‑76 (MQ=255)                                                            tGAATTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa              <  2:928759/76‑1 (MQ=255)                                                               aTTTCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg           >  2:2909088/1‑76 (MQ=255)                                                                 ttCCGTTAGAACAA‑AACAGCCGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt         <  1:1866419/76‑1 (MQ=255)                                                                  tCCGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt        <  2:4589161/76‑1 (MQ=255)                                                                   ccGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa       <  2:1478652/76‑1 (MQ=255)                                                                   ccGTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa       >  1:820182/1‑76 (MQ=255)                                                                   ccGTTAGAAC‑ACAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAACCAGTCTGCCATACCGAAAAGGATTa       <  1:504948/76‑1 (MQ=255)                                                                     gTTAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa     <  1:549455/76‑1 (MQ=255)                                                                      ttAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAg    <  2:70616/76‑1 (MQ=255)                                                                       tAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAgg   >  2:4921635/1‑76 (MQ=255)                                                                       tAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAgg   >  2:241664/1‑76 (MQ=255)                                                                       tAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAgg   >  2:2060945/1‑76 (MQ=255)                                                                       tAGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAgg   >  2:2041746/1‑76 (MQ=255)                                                                        aGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt  <  1:3332347/76‑1 (MQ=255)                                                                        aGAACAA‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt  >  2:5784020/1‑76 (MQ=255)                                                                             |                                                                        TACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑A‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGT  >  NC_000913/4435606‑4435751

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