Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,396,962	T→G	14.4%	R278R (CGA→CGC)	rng ←	RNase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,396,962	0	T	G	14.4%	88.7 / 9.3	52	R278R (CGA→CGC)	rng	RNase G
Reads supporting (aligned to +/- strand):  ref base T (16/26);  new base G (4/4);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01

AATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCG  >  NC_000913/3396821‑3397102                                                                                                                                              |                                                                                                                                             aaTATTGGTATTGAAAATGTTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctc                                                                                                                                             <  2:488514/143‑1 (MQ=255)   tatTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctcgc                                                                                                                                           >  1:293735/1‑143 (MQ=255)        gTATTGAAAATGGTGTCTTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGAGCGCTGGaa                                                                                                                                      <  1:106276/143‑2 (MQ=255)        gTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTTCACGGTGGCCATCGCTTCGGTATGGTCGTTAATGAGATCAGCACCGGATTTCAGTTCTACTGTGCGTTCCAGcgctcgtcggtt                                                                                                                                      >  2:369083/1‑137 (MQ=255)         tATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACCAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGAtt                                                                                                                                     <  2:266241/143‑1 (MQ=255)         tATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTGCAGCGCTCGCTGGAtt                                                                                                                                     >  1:282371/1‑143 (MQ=255)         tATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGAtt                                                                                                                                     >  1:568678/1‑143 (MQ=255)            tGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAAGGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCg                                                                                                                                  <  1:576964/143‑1 (MQ=255)             gAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGACATCGCTTAGGTCTGATCGATAATGAGATAACCACCTGATCTCAGCTGTACTCTGCGTTCCAGCGCTCGCTggccttggg                                                                                                                                 >  2:1026258/1‑136 (MQ=255)                agtGGTGTCGTCCAGGTTGCGATGACCCACCAACGCCCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGtttt                                                                                                                              <  1:470339/141‑1 (MQ=255)                 ttGGGGTCGTCGAGACTGCCAAGACCGACAAACGCCCCGGGACCGACGCCCACGGCGGTCATCGCCTCGGTCTGGTCGTTAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTc                                                                                                                             <  1:791605/142‑1 (MQ=255)                    gtgtCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGAc                                                                                                                          <  1:200341/143‑1 (MQ=255)                       tcgtcCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCATCGGTCTGGTCGATAATGAGATAACCACCGGAATCCCGTTCAGGCAATCGTTGCAGAGCCCGCGCGATTGCGTTTTCGACATc                                                                                                                       >  2:464052/1‑143 (MQ=255)                       tcgtcCAGATCGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATc                                                                                                                       <  1:251518/143‑1 (MQ=255)                          ggcAGTTTGCGAGGACTGACGTGCGGTCGGGTGTTGATGTCCACGGTGGTCGTCGCTTCGGTCTGGTAGGTAATGGGATAACTACGGGATTTCAGTTCTACTTTGCGTTCCAGCGCTGGCTGGATTTCGTTTTCGACATCaaa                                                                                                                    <  2:691607/141‑1 (MQ=255)                           ccAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGACTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGACCCCGCCATCAAag                                                                                                                   >  1:229770/1‑143 (MQ=255)                                 tGCGATGACCGCCACACGCTGCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCg                                                                                                             <  2:326945/143‑1 (MQ=255)                                    gATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGAGCGCTGGaa                                                                                                                                      >  1:831378/1‑114 (MQ=255)                                    gATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGAGCGCTGGaa                                                                                                                                      <  2:831378/115‑2 (MQ=255)                                              aaaCGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGc                                                                                                <  1:176600/143‑1 (MQ=255)                                              aaaCGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGACCGACACTCGGCTgcc                                                                                                >  2:965658/1‑141 (MQ=255)                                                    tCCGGTATTGCTGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCtg                                                                                          <  2:63020/143‑1 (MQ=255)                                                     ccGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCtgt                                                                                         <  1:316680/143‑1 (MQ=255)                                                           tctaTGGCCACGGTGGTCATCGCGTCGGTCTGGTCGATACTGATATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATg                                                                                   <  2:885987/140‑1 (MQ=255)                                                             gaggTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAACGAGATAGCCACCGGATTTCAGTGCTACTTGGCGTCCGAGCGGGCGCTGGATTTCGTTTTCGAGATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCt                                                                                 <  2:238876/140‑1 (MQ=255)                                                                  ccACGGTTGTGATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgc                                                                            <  1:462754/143‑1 (MQ=255)                                                                   cACTGTGGTCATCGCCTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgct                                                                           <  1:881409/143‑1 (MQ=255)                                                                           tCATCGCTTCGTTCTGGGCGATAATGAAATAACCACCGGATTACAGTTCTCCTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTAGACATCAAAGCGGTCGAGACGCGGCTGGCGGCGTGTGTAATGCTGACGCTTGATTGTc                                                                   >  2:819408/1‑143 (MQ=255)                                                                           tCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACATTGCGTTCCAGCGATCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAAACGGCTGGCGGCCTGTGTAGTGGTGCAgctggcaatgc                                                                   >  2:819282/1‑138 (MQ=255)                                                                             aTCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAt                                                                 >  1:560638/1‑143 (MQ=255)                                                                              tCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtc                                                                >  1:84212/1‑143 (MQ=255)                                                                                   tCGGTCTGGTCGATAATGAGGTAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCggg                                                           <  2:828496/143‑1 (MQ=255)                                                                                    cGGTCTGGGCGATAATGAGATAACGGCCGGATTTGAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTCCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGa                                                          <  2:572662/143‑1 (MQ=255)                                                                                      gTCTTGTCCATAATGAGATAACCACCTGATTTCAGTTTTACTTTGCGTTTCTGCGGGCGCTGGATCTCCTCTTCGCCATCCAAGAGCTCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAAt                                                        <  1:767955/143‑1 (MQ=255)                                                                                       tcggggcGATAATTAGATCACCACCGGATTTCAGTCCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGACCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATg                                                       <  1:369083/137‑1 (MQ=255)                                                                                           gTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACt                                                   >  1:920437/1‑143 (MQ=255)                                                                                                 aaTGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAgg                                             <  1:493875/143‑1 (MQ=255)                                                                                                  aTGAGATAACCACCGGATTTCACTTCTACGTTGCGTTCCAGCGCTTGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGt                                            <  1:143774/143‑1 (MQ=255)                                                                                                      gaTCACCACCGGATTTCAGTTCGACCTGGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGCCATCAAATAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTTCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAAc                                        <  1:128296/143‑1 (MQ=255)                                                                                                          agcactGGATTGGAGTTCTACGTTGCGTTCCAGCGCTCGCGGGATTTCGTTTTCCGCATCAAAGGGATCGCAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCaa                                    <  2:68275/141‑1 (MQ=255)                                                                                                            caccGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGt                                  >  2:589831/1‑143 (MQ=255)                                                                                                             accGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGACGTGACCTCAAGTa                                 >  1:78364/1‑143 (MQ=255)                                                                                                             accGGATTTCAGTTCTACTTTGCGTTCCAGCGCGCGCTGGAGTTCGTTTTCGGCATCAAAGAGATCGACAATCGCCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCCGGAATGTCCCCCGCGGTGCACTCAAGTa                                 >  1:477091/1‑143 (MQ=255)                                                                                                                         ttCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAAATCCTCGTAAg                     >  1:247896/1‑143 (MQ=255)                                                                                                                           cTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTCATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGCAGAGCACGAAAGTAAGCCATCGTAAGTc                   >  2:684167/1‑143 (MQ=255)                                                                                                                            tACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCa                  <  1:882956/143‑1 (MQ=255)                                                                                                                             aCTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAATTCAAGTAACGCTTCGTAAGTCAg                 >  1:413319/1‑143 (MQ=255)                                                                                                                             aCTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAg                 <  1:63962/143‑1 (MQ=255)                                                                                                                               tttGCGTTCCAGTGCTCGCTGGGTTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTACTGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGc               <  2:406072/143‑1 (MQ=255)                                                                                                                                   cGTTCCAGCACTCGCTGGATTTCGGGTTCGACATCAAAGAGATCGAAAATCGTCTGTCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGa           <  2:58997/143‑1 (MQ=255)                                                                                                                                   cGTGCCAGCGCCTGCTGGATTTCGTTTTCGACATCAAAGAGAGCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGa           <  2:722750/143‑1 (MQ=255)                                                                                                                                   aaatCCAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCCTATAGGTTGTGTACTCCGAGGCGGAGCCAAGTAACGCTTCGTAAGTCTGGCGTGa           >  2:417676/4‑143 (MQ=255)                                                                                                                                   aaatCCAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGCGAATGTACTCCGCGGTGAACTCAAGTAACGCTTCCTAAGTCAGGCATGa           >  2:93023/4‑143 (MQ=255)                                                                                                                                      tCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAGCGCTTCGTAAGTCAGGCGTGAGTc        >  1:643253/1‑143 (MQ=255)                                                                                                                                         aGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGACTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAc     <  1:371882/143‑1 (MQ=255)                                                                                                                                            gctcgctGGATTTCGTTTTTGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGCATGTCCTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAcgcg  >  2:427064/1‑143 (MQ=255)                                                                                                                                              |                                                                                                                                             AATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCG  >  NC_000913/3396821‑3397102

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