Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,474,826	+A	intergenic (+1/+44)	yjgL → / ← argI	hypothetical protein/ornithine carbamoyltransferase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,474,823	1	.	A	97.5%	253.9 / 0.1	81	coding (1849/1851 nt)	yjgL	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (43/36);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAA‑TAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCAT  >  NC_000913/4474751‑4474893                                                                          |                                                                        aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaa                                                                      <  1:4851959/76‑1 (MQ=255) aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaa                                                                      >  1:5063918/1‑76 (MQ=255) aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaa                                                                      <  2:1312299/76‑1 (MQ=255) aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaa                                                                      >  2:4316496/1‑76 (MQ=255) aaaCTAAAGCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaa                                                                      <  2:1317976/76‑1 (MQ=255)   aCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaat                                                                    >  2:302565/1‑74 (MQ=255)       aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAAtg                                                                >  1:3449116/1‑76 (MQ=255)          ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGa                                                             >  1:3524755/1‑76 (MQ=255)           cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGAt                                                            >  2:4725805/1‑76 (MQ=255)           cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGAt                                                            >  2:5417071/1‑76 (MQ=255)           cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGAt                                                            <  1:311250/76‑1 (MQ=255)            aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATg                                                           <  2:1511079/76‑1 (MQ=255)             ccGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGc                                                          <  2:3524755/76‑1 (MQ=255)              cGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGcc                                                         >  1:5437992/1‑76 (MQ=255)               ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCg                                                        <  2:3627874/76‑1 (MQ=255)               ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCg                                                        >  1:2857855/1‑76 (MQ=255)               ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCg                                                        <  1:4104643/76‑1 (MQ=255)                  aaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGa                                                     >  1:426039/1‑76 (MQ=255)                  aaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGa                                                     >  1:1144610/1‑76 (MQ=255)                   aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGAt                                                    >  2:4396319/1‑76 (MQ=255)                   aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGAt                                                    <  2:795213/76‑1 (MQ=255)                    aTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATg                                                   <  1:4106289/76‑1 (MQ=255)                     tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGACGGGATgg                                                  >  1:803071/1‑76 (MQ=255)                      ttACGTATAGCAATCGACCACTTTGGGTGGACTAATATACCTATCACTGCAT‑AAACTCATGTGATGCCGGTATGGt                                                 <  1:839314/76‑1 (MQ=255)                      ttACGTATAGCAATAGAACACTTTGGGAGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGt                                                 >  2:5061648/1‑76 (MQ=255)                       tACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGtt                                                <  1:3729199/76‑1 (MQ=255)                       tACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGtt                                                >  1:90948/1‑76 (MQ=255)                        aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGttt                                               <  2:5074466/76‑1 (MQ=255)                         cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTg                                              <  2:2037453‑M1/76‑2 (MQ=255)                          gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgt                                             >  1:2500866‑M1/1‑74 (MQ=255)                          gTATAGCAATAGAACACTTTGGGCGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgt                                             >  2:1999166‑M1/1‑74 (MQ=255)                           tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgta                                            <  1:2697826‑M1/76‑4 (MQ=255)                            ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtat                                           <  2:3449116‑M1/76‑5 (MQ=255)                            ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtat                                           >  1:5286088‑M1/1‑72 (MQ=255)                             taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatt                                          <  1:1685334‑M1/76‑6 (MQ=255)                               gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttc                                        <  1:5047611‑M1/76‑8 (MQ=255)                               gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttc                                        <  2:5286087‑M1/76‑8 (MQ=255)                                cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcc                                       >  2:5170908‑M1/1‑68 (MQ=255)                                cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcc                                       >  2:1065605‑M1/1‑68 (MQ=255)                                 aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttccc                                      <  1:482835‑M1/76‑10 (MQ=255)                                 aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttccc                                      <  2:716706‑M1/76‑10 (MQ=255)                                 aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttccc                                      >  1:4083442‑M1/1‑67 (MQ=255)                                  aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccg                                     >  1:1273664‑M1/1‑66 (MQ=255)                                  aTAGAACACTTTGGGCGGAACAATAGAACTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccg                                     >  2:4153916‑M1/1‑66 (MQ=255)                                   tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccgg                                    <  1:290063‑M1/76‑12 (MQ=255)                                   tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccgg                                    >  1:2825938‑M1/1‑65 (MQ=255)                                    aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggc                                   <  1:4083800‑M1/76‑13 (MQ=255)                                    aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggc                                   >  2:5688683‑M1/1‑64 (MQ=255)                                    aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggc                                   <  2:2197713‑M1/76‑13 (MQ=255)                                     gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggca                                  >  2:4142820‑M1/1‑63 (MQ=255)                                      aaCACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcat                                 <  2:520571‑M1/76‑15 (MQ=255)                                         acTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                              >  1:5471607‑M1/1‑59 (MQ=255)                                           tttGTGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                            <  2:2236662‑M1/76‑20 (MQ=255)                                               ggTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                        <  1:4543673‑M1/76‑24 (MQ=255)                                               ggTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                        <  1:1307822‑M1/76‑24 (MQ=255)                                               ggTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                        >  1:415387‑M1/1‑53 (MQ=255)                                                 tGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                      >  1:1696397‑M1/1‑51 (MQ=255)                                                  gtaagaataGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                     <  1:2861330‑M1/74‑27 (MQ=255)                                                   gaagaaTAGACCTATCGCTGCAT‑AAAGTAATGTGATGCCGGGATGGTTTgtatttcccgtcatctttatagcgata                    <  2:1296643‑M1/76‑28 (MQ=255)                                                   gaagaaTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                    <  1:3289137‑M1/76‑28 (MQ=255)                                                   gaagaaTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                    <  1:510791‑M1/76‑28 (MQ=255)                                                     agaaTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                  >  2:3669994‑M1/1‑47 (MQ=255)                                                      gaaTAGACCTATCACTGCATAAAA‑TAATGTGGTGAAGGGGTGGTTTgtatttcacggcatctttatagcgatagca                 >  2:3446866‑M1/1‑46 (MQ=255)                                                      gaaTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                 <  1:2055278‑M1/76‑31 (MQ=255)                                                       aaTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatccttataccgacagcaa                >  2:738863‑M1/1‑45 (MQ=255)                                                        aTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat               >  1:1722902‑M1/1‑44 (MQ=255)                                                        aTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat               >  1:1685020‑M1/1‑44 (MQ=255)                                                          aGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatta             >  2:2719944‑M1/1‑42 (MQ=255)                                                            aCCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt           >  1:1971189‑M1/1‑40 (MQ=255)                                                             ccTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattt          >  2:2879384‑M1/1‑39 (MQ=255)                                                              cTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattta         >  2:417158‑M1/1‑38 (MQ=255)                                                               tATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac        >  2:5667976‑M1/1‑37 (MQ=255)                                                               tATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac        >  1:5258114‑M1/1‑37 (MQ=255)                                                                aTCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact       <  1:2693641‑M1/76‑41 (MQ=255)                                                                aTCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatt                                          >  1:663716‑M1/1‑36 (MQ=255)                                                                aTCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatt                                          <  2:663716‑M1/41‑6 (MQ=255)                                                                 tCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactg      >  2:3497901‑M1/1‑35 (MQ=255)                                                                 tCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagagatagcaattctttactg      >  2:3540028‑M1/1‑35 (MQ=255)                                                                  cACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga     >  1:1116388‑M1/1‑34 (MQ=255)                                                                  cACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga     >  1:4512257‑M1/1‑34 (MQ=255)                                                                  cACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga     <  1:382259‑M1/76‑43 (MQ=255)                                                                  cACTGCATAAAA‑TAATGGGATGCCGGGAGGGTTTgtattgcccggcatctttatagcgatagcaattatttactga     <  2:200589‑M1/76‑43 (MQ=255)                                                                   aCTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag    >  2:4259931‑M1/1‑33 (MQ=255)                                                                   aCTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag    <  2:415387‑M1/76‑44 (MQ=255)                                                                   aCTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag    <  1:4279716‑M1/76‑44 (MQ=255)                                                                     tGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcg  <  2:3442830‑M1/76‑46 (MQ=255)                                                                     tGCATAAAA‑TAATGGGATGCCGGGATGGTTTgtattgcccggcatctttatagcgatagcaattatttactgagcg  <  2:3483138‑M1/76‑46 (MQ=255)                                                                          |                                                                        AAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAA‑TAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCAT  >  NC_000913/4474751‑4474893

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