Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP028101	761,062	A→G	K479K (AAA→AAG)	C7Y58_03755 →	hemolysin activator protein precursor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP028101	761,062	0	A	G	98.0%	132.5 / NA	50	K479K (AAA→AAG)	C7Y58_03755	hemolysin activator protein precursor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (22/27);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATAGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAATTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTT  >  CP028101/761013‑761207                                                  |                                                                                                                                                  caanggaaatagaaaATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTatat                                               >  1:3364608/6‑150 (MQ=0) caanggaaatagaaaATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTatat                                               >  1:3364609/6‑150 (MQ=255)                              agagTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATattta                  <  2:3314015/150‑1 (MQ=11)                              agagTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATattta                  <  2:3294064/150‑1 (MQ=11)                              agagTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATattta                  <  1:3502690/150‑1 (MQ=11)                              agagTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATattta                  <  1:2658745/150‑1 (MQ=11)                              agagTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATattta                  <  2:2078350/150‑1 (MQ=11)                              agagTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATattta                  <  1:2855324/150‑1 (MQ=11)                               gagTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttag                 <  1:1422050/150‑1 (MQ=14)                               gagTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttag                 <  2:3082940/150‑1 (MQ=14)                                agTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttaga                <  1:1633381/150‑1 (MQ=17)                                agTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttaga                <  1:2160952/150‑1 (MQ=17)                                 gtgtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagat               <  1:3260318/150‑1 (MQ=17)                                 gtgtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagat               >  2:2724771/1‑150 (MQ=17)                                 gtgtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagat               >  2:3263273/1‑150 (MQ=17)                                  tgtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagatt              >  2:2416890/1‑150 (MQ=18)                                   gtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCGTATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattt             >  1:1887807/1‑150 (MQ=19)                                   gtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattt             >  2:250478/1‑150 (MQ=21)                                   gtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattt             >  1:3468665/1‑150 (MQ=21)                                   gtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattt             <  1:2902529/150‑1 (MQ=21)                                   gtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattt             >  1:1887808/1‑150 (MQ=21)                                   gtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagactt             >  1:3468666/1‑147 (MQ=17)                                   gtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTAGTCAGAGGGAAAAGTAAGAAATAACAAAAAGAATTCAAAATAGAGAAAAGGAGATATGTCAGGTGCTATATTAGGGTTAAGATATAATATGAAATatttagattt             <  2:3393534/150‑1 (MQ=11)                                    ttAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattta            <  2:1663789/150‑1 (MQ=25)                                    ttAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattta            <  2:881539/150‑1 (MQ=25)                                    ttAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattta            <  1:2873720/150‑1 (MQ=25)                                    ttAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattta            <  1:2873725/150‑1 (MQ=25)                                    ttAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattta            <  1:2873726/150‑1 (MQ=25)                                    ttAGACTCAGAAAAGTTTGGAAAAACAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCATGTGCTATATTAGGTTTAAGATATAATATGAAATatttagattta            <  1:2873719/150‑1 (MQ=255)                                     tAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagatttag           <  1:2576968/150‑1 (MQ=25)                                     tAGACTCAGAAAAATTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATatttagatttag           <  1:3354577/150‑1 (MQ=35)                                      aGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATAtttagatttaga          <  2:1159609/150‑1 (MQ=37)                                      aGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATAtttagatttaga          <  2:1583849/150‑1 (MQ=37)                                      aGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATAtttagatttaga          <  1:1470424/150‑1 (MQ=37)                                       gACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATctagatttagat         >  2:3419903/1‑150 (MQ=255)                                        aCTCAGAAGAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATg        <  1:660001/150‑1 (MQ=255)                                        aCTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGGTTTAGATg        >  1:3254427/1‑150 (MQ=255)                                        aCTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATg        >  1:3254428/1‑150 (MQ=255)                                        aCTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATg        <  2:794642/150‑1 (MQ=255)                                        aCTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATg        <  2:31417/150‑1 (MQ=255)                                         cTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGt       >  2:1643951/1‑150 (MQ=255)                                         cTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGt       >  2:1643950/1‑150 (MQ=255)                                         cTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGt       >  1:942711/1‑150 (MQ=255)                                         cTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGt       >  2:368829/1‑150 (MQ=255)                                          tCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTa      >  1:2307735/1‑150 (MQ=255)                                          tCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTa      <  1:859746/150‑1 (MQ=255)                                             gAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAACTAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCt   >  1:2684226/1‑150 (MQ=255)                                             gAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCt   >  2:2911394/1‑150 (MQ=255)                                             gAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCt   >  1:2684231/1‑150 (MQ=255)                                             gAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCt   >  1:1158219/1‑150 (MQ=255)                                              aaaagTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAAGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCtt  >  1:2066761/1‑150 (MQ=255)                                                  |                                                                                                                                                  CAATAGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAATTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTT  >  CP028101/761013‑761207

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