Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,375,713	T→A	7.7%	I280I (ATT→ATA)	sbcC →	exodeoxyribonuclease SbcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,375,713	0	T	A	7.7%	157.0 / 2.8	52	I280I (ATT→ATA)	sbcC	exodeoxyribonuclease SbcC
Reads supporting (aligned to +/- strand):  ref base T (21/27);  new base A (2/2);  total (23/29)
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.87e-01

GATTGAAAATATACCACTTTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATATTGAAAATAGA  >  CP000730/1375575‑1375842                                                                                                                                           |                                                                                                                                  gattgaAAATATACCACTTTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATAAAAAGATAATATTaa                                                                                                                                 >  1:89590/1‑141 (MQ=255)    tgaAAATATACCACTTTTACAAGCTAGGTCGAAAGAAAGACCAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAAtca                                                                                                                              <  1:513293/141‑1 (MQ=255)      aaaaTATACCACTTTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAAtcaac                                                                                                                            >  1:549728/1‑141 (MQ=255)           ataCCACTTTTACAAGCTTGTGCTAAATAAATACTAGCATCTGGAAATGAAAGTAAAGAAACGGCTATTAAAGAATTTGAAATAATAGAAGAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAAtcaactcaac                                                                                                                       <  2:298917/141‑1 (MQ=255)           ataCCACATTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAAtcaactcaac                                                                                                                       <  1:282790/141‑1 (MQ=255)              ccACTTTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACaaa                                                                                                                    <  1:597379/141‑1 (MQ=255)                      aCAAACGAGGTAGAAAGAAATAATAACATGTGTAAATGCAAGTAAAGAAACTGATATGAACGAACTTGAAATATTAGAAAAAACAACAGTAGAAAATAATATATTAAAAGATAATATTAATCAACACAAAAaaaataaaat                                                                                                            >  2:232688/1‑141 (MQ=255)                       cAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATGATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAAtt                                                                                                           >  1:74577/1‑141 (MQ=255)                        aaGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATCTTAAAAGATAATATAAATCAACTCAACAAAAATAAAATTg                                                                                                          >  1:512963/1‑141 (MQ=255)                        aaGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATAAATCAACTCAACAAAAATAAAATTg                                                                                                          >  1:513026/1‑141 (MQ=255)                          gCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGAt                                                                                                        >  1:502329/1‑141 (MQ=255)                              ggTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCg                                                                                                    >  2:582079/1‑141 (MQ=255)                              ggTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCg                                                                                                    >  1:591689/1‑141 (MQ=255)                                 cGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTc                                                                                                 >  1:201631/1‑141 (MQ=255)                                   aaagaaaTACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCaa                                                                                               <  2:326911/141‑1 (MQ=255)                                      gaacTATTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAAGATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTg                                                                                            <  1:41661/141‑1 (MQ=255)                                       aaaTACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGa                                                                                           >  2:168674/1‑141 (MQ=255)                                          tACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAg                                                                                        <  1:613694/141‑1 (MQ=255)                                          tACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAg                                                                                        >  1:280181/1‑141 (MQ=255)                                           aCTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGa                                                                                       >  1:309414/1‑141 (MQ=255)                                                   cTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGaacaacaac                                                                               <  1:555651/141‑1 (MQ=255)                                                    tGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAATATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGaacaacaaca                                                                              >  1:10120/1‑140 (MQ=255)                                                    tGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAAcc                                                                              <  2:318576/141‑1 (MQ=255)                                                       aaaTGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGa                                                                           >  2:287778/1‑141 (MQ=255)                                                         aTGAAAGTAAAGAAACGTGTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGaaa                                                                         <  1:143367/141‑1 (MQ=255)                                                            aaaGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATag                                                                      <  1:449022/141‑1 (MQ=255)                                                                tAAAGAAACTGCTATTAAAAAATTTGAAATAATAGAAAATAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCttacatcgggat                                                                  <  2:335956/141‑12 (MQ=255)                                                                 aaaGAAACTGCTATTAAAGGATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGa                                                                 >  1:398004/1‑141 (MQ=255)                                                                  aaGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAACGAACAACAACCTGAAATAGAGGGaa                                                                >  2:16733/1‑141 (MQ=255)                                                                                aaaGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTa                                                  <  1:410100/141‑1 (MQ=255)                                                                                   aaaattGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAAATGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAg                                               <  2:123988/137‑1 (MQ=255)                                                                                    aaTTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAgt                                              <  2:74577/141‑1 (MQ=255)                                                                                     aTTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAgtt                                             <  2:385398/141‑1 (MQ=255)                                                                                       ttGAAATAATATAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAgttgt                                           <  2:264487/141‑1 (MQ=255)                                                                                          aaaTAATAGAAAACAAAACATTAGAAAATAATATATTAAAAGATAAGATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTAc                                        <  1:606725/141‑1 (MQ=255)                                                                                                taGAAAAGAAATCAATAGAAAATAATAGATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGata                                  <  2:304231/141‑1 (MQ=255)                                                                                                    aaagaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTAc                              <  1:171803/141‑1 (MQ=255)                                                                                                     aagaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACa                             <  1:615806/141‑1 (MQ=255)                                                                                                     aagaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACa                             >  2:296829/1‑141 (MQ=255)                                                                                                     aagaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACa                             >  2:189718/1‑141 (MQ=255)                                                                                                       gaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACaaa                           <  1:287778/141‑1 (MQ=255)                                                                                                       gaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACaaa                           >  1:429676/1‑141 (MQ=255)                                                                                                        aaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAAc                          <  1:163744/141‑1 (MQ=255)                                                                                                        aaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAAc                          <  1:234305/141‑1 (MQ=255)                                                                                                                   aaaTAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATta               >  2:583175/1‑141 (MQ=255)                                                                                                                    aataatATATTAAAAGATAATATAAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATtat              <  2:512963/141‑1 (MQ=255)                                                                                                                    aataatATATTAAAAGATAATATAAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATtat              <  2:513026/141‑1 (MQ=255)                                                                                                                    aataatATATAAAAAGATAATCTTATGCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATtat              >  1:661074/1‑141 (MQ=255)                                                                                                                       aataTATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATAtt           <  2:251367/141‑1 (MQ=255)                                                                                                                       aataTATTAAAAGATAATATTAATCAACTCAACAAAAATAAAAATGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATAtt           <  2:658243/141‑1 (MQ=255)                                                                                                                             ttAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATATTGAAAAt     <  1:168674/141‑1 (MQ=255)                                                                                                                                aaaGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATATTGAAAATaga  >  2:588704/1‑141 (MQ=255)                                                                                                                                           |                                                                                                                                  GATTGAAAATATACCACTTTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATATTGAAAATAGA  >  CP000730/1375575‑1375842

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: ‑