Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	4,546,084	G→T	H18Q (CAC→CAA)	yjiC ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,546,084	0	G	T	98.3%	190.8 / ‑5.8	59	H18Q (CAC→CAA)	yjiC	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (22/36);  minor base C (0/1);  total (22/37)
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.79e-01

GTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGC  >  CP009273/4545939‑4546226                                                                                                                                                  |                                                                                                                                               gTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGa                                                                                                                                             <  2:175956/149‑1 (MQ=255)       tAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCt                                                                                                                                       <  2:110783/149‑1 (MQ=255)       tAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCACCt                                                                                                                                       >  2:518236/1‑149 (MQ=255)              gACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGCGTAAACTCACGGTCGGCTATTGGGTGATGAAGGAATTGAGCATCGGtaataa                                                                                                                                >  2:499857/1‑149 (MQ=255)                 aaCGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGtaataa                                                                                                                                >  1:270087/1‑146 (MQ=255)                     cTTTTTCCATCGTAATTATGGTTATTGGCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGCct                                                                                                                         <  2:119754/149‑1 (MQ=255)                     cTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGAct                                                                                                                         <  2:910990/149‑1 (MQ=255)                        tttCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGAct                                                                                                                         <  1:631406/146‑1 (MQ=255)                            cATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGa                                                                                                                     >  1:839280/1‑146 (MQ=255)                            cATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACAATTGGGTGGTGAAGGAATTGAGCACATGTACTAATTTGACTCTGa                                                                                                                     >  1:542401/1‑146 (MQ=255)                             aTCGTATTTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGTGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGaa                                                                                                                    <  1:250594/146‑1 (MQ=255)                             aTCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGaa                                                                                                                    >  1:985145/1‑146 (MQ=255)                               cGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCa                                                                                                               <  2:850392/149‑1 (MQ=255)                               cGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCa                                                                                                               <  2:850697/149‑1 (MQ=255)                                   aTTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTACACTCACTGTCGGCTAGTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCAt                                                                                                              <  1:904026/146‑1 (MQ=255)                                       tGGTTATTGGCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCt                                                                                                          <  1:5071/146‑1 (MQ=255)                                        ggTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAAt                                                                                                      <  2:365778/149‑1 (MQ=255)                                        ggTTATTGCCACCACTAAGAGCGGGAATATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAAt                                                                                                      <  2:22902/149‑1 (MQ=255)                                         gTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATg                                                                                                     <  2:1109533/149‑1 (MQ=255)                                            aTTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCACGCGTGCATATTTTAATTTACGTGTAAACCCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATg                                                                                                     <  1:302445/146‑1 (MQ=255)                                             ttGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATgg                                                                                                    <  1:77397/146‑1 (MQ=255)                                             ttGCCACGAGTAAGAGGCGGAACATAGATGTTTTCCGAATTCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTTAAGGGCATTGCTTAATgg                                                                                                    <  1:481213/146‑1 (MQ=255)                                                     aGTAAGAGGCGGAAGATAGATGTTTTCCGCATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTGGAGTATTGGGTGATGAAGGAATTTGGCATCTTTAATAATTTGACTCTGACGGGCATTGCTTAGTGTCATAGACa                                                                                            >  1:284288/1‑146 (MQ=255)                                                     aGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAATCTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACataa                                                                                         <  2:568102/149‑1 (MQ=255)                                                     aGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGGAGGGCATTGCTTAATGGCATAGACa                                                                                            >  1:421204/1‑146 (MQ=255)                                                          gagGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACataata                                                                                       <  1:146342/146‑1 (MQ=255)                                                           agGCGGAAGATGGATGTTTTCCGAATCCATGCGTGGATATTTTAATTTACGGGTAAACTCACTGTCGTCTATTGGGTGATGAGGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATTGACATAATAg                                                                                      <  1:719631/146‑1 (MQ=255)                                                             gCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGtt                                                                                    <  1:1000706/146‑1 (MQ=255)                                                               tgAAGAGAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTAc                                                                                  <  1:467012/145‑1 (MQ=255)                                                               ggAAGAGAGATGTTTTCCCAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTAc                                                                                  <  1:810613/146‑1 (MQ=255)                                                                   gatagatGTTTTCCGAATCCATGCGGGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTAt                                                                           >  2:938469/1‑149 (MQ=255)                                                                        atGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGa                                                                         <  1:425367/146‑1 (MQ=255)                                                                           ttttCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACAGAATAGTTACCTTTTATGAGTa                                                                      >  1:835634/1‑146 (MQ=255)                                                                                 gAATCCAGGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATGGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCAc                                                                <  1:326908/146‑1 (MQ=255)                                                                                  aaTCCATTCGTTCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGTGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAATAttcccct                                                               <  1:913226/146‑6 (MQ=255)                                                                                       aTGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTa                                                       >  2:18905/1‑149 (MQ=255)                                                                                               atatTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATGTGTAATCATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGGTGTTTAGaaaa                                                  >  1:772996/1‑146 (MQ=255)                                                                                                 atTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAACTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAAt                                                <  1:570973/146‑1 (MQ=255)                                                                                                         ttACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAtttt                                     >  2:181937/1‑149 (MQ=255)                                                                                                         ttACGTGTAAACTCACTGTCGACGATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAt                                        <  1:995387/146‑1 (MQ=255)                                                                                                            cGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTTCCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCt                                  <  2:670525/149‑1 (MQ=255)                                                                                                            cGTGTAAACTCACGGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCt                                  <  2:506430/149‑1 (MQ=255)                                                                                                             gtgtAACCTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAttttt                                    <  1:449078/146‑1 (MQ=255)                                                                                                              tgtAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTc                                   >  1:894247/1‑146 (MQ=255)                                                                                                                tAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgtag                              >  2:1077003/1‑149 (MQ=255)                                                                                                                    cTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGCTAAATCTTGCTgtagt                             <  1:910368/146‑1 (MQ=255)                                                                                                                     tCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTaaaa                         >  2:1003318/1‑149 (MQ=255)                                                                                                                      cACTGCCAACTATTTGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGAATAAAAATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTaaaaa                        <  2:245828/149‑1 (MQ=255)                                                                                                                       aCTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTCGTaaa                          >  1:25765/1‑146 (MQ=255)                                                                                                                        cTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAATAATAGATAAATTTTTCTGTAGTAAAAaga                      >  2:30796/1‑149 (MQ=255)                                                                                                                         ttgcgatgaTTGGTGGATGATGGAATTAAGCATCTGTAATAATTTGACTCTTAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTaaaaa                        <  1:100656/138‑1 (MQ=255)                                                                                                                         gggcGACTATTGGGTCATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAagag                     <  2:393997/146‑1 (MQ=255)                                                                                                                           tCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAaga                      <  1:706835/146‑1 (MQ=255)                                                                                                                           tCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGaa                   <  2:856801/149‑1 (MQ=255)                                                                                                                                     ggTGATGAAGGAATTGAGCATCTGTCATAATTTGACTCTGAAGGGCATTACTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTCCTGTAGTAACAAGAGAAGTAAACa            >  1:895784/1‑146 (MQ=255)                                                                                                                                     ggTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAATATAGATACATTTTTCGGGAGTAAAAAGAGAAGTAAACAAAt         >  2:264318/1‑149 (MQ=255)                                                                                                                                     ggTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACa            <  1:859655/146‑1 (MQ=255)                                                                                                                                       tgatgaAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGa       >  2:1112874/1‑149 (MQ=255)                                                                                                                                         atgaAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAca     >  2:877666/1‑149 (MQ=255)                                                                                                                                            aaggaatTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatgc  >  2:104713/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                               GTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGC  >  CP009273/4545939‑4546226

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