Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	1,252,163	T→G	28.1%	intergenic (+28/‑169)	U712_06305 → / → U712_06310	Uncharacterized protein yjgD/putative lipoprotein yjhA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	1,252,163	0	T	G	28.1%	46.5 / 4.3	41	intergenic (+28/‑169)	U712_06305/U712_06310	Uncharacterized protein yjgD/putative lipoprotein yjhA
Reads supporting (aligned to +/- strand):  ref base T (8/21);  new base G (8/3);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 1.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATACAAAGGAGCTGGAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATGAGATAGAAATAATCGTTACTTTCTTACCTGTTAAGAAAAGGGAAAATCTGATTATAATAG  >  CP006881/1252040‑1252301                                                                                                                            |                                                                                                                                           gATACAAAGGAACTGGAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGACCGGGTTCAGCCAGCGCCGGGTTTTtgatg                                                                                                                              >  2:210061/1‑138 (MQ=255)  aTACAAAGGAGCTGGAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGaa                                                                                                                            <  3:61314/139‑1 (MQ=255)   tACAACGGAGCTGGACCGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGaa                                                                                                                            <  3:389610/138‑1 (MQ=255)   tACAAAGGAGCTGGAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGACCCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGaaa                                                                                                                           <  7:303520/139‑1 (MQ=255)   tACAAAGGAGCTGGAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGaaa                                                                                                                           <  3:181105/139‑1 (MQ=255)   acaaaaCGAGCTGGAACCCACGTCGCAGCCTCCGGAAAACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGTATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGGTCGGGTTTTTTGATGaaa                                                                                                                           <  4:38420/136‑1 (MQ=255)      aaaGGAGCTGGAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCCGGTTTTTTGATGAAAAgg                                                                                                                        >  5:375682/1‑139 (MQ=255)           aGCTGGAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGCGCCGCTTCCGCCCCCGCGGGGTTTTTGATGTAAAGGtggaa                                                                                                                   >  4:259589/1‑135 (MQ=255)             cTGGAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGACCGGCTTCAGCCCGCGCGGGTTTTTTGTTGAAAATGTCTTATg                                                                                                                 >  6:208816/1‑139 (MQ=255)             cTGGAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCCCGCCGCATCAGCCCGCTCGGGGTTTTTTATGAAAAAGTGTAAtt                                                                                                                 >  6:322058/1‑138 (MQ=255)             cTGGAACGCACGACGCAGCCTCCGGAACACC‑AAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGACCGGCTTCAGCCAGCGCGGGGTTTTTGATGAAAAGGTCGAATgg                                                                                                                >  7:79242/1‑139 (MQ=255)             cTGGAACGCACGACGCAGCCGCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATg                                                                                                                 <  5:322058/139‑1 (MQ=255)                gAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGa                                                                                                               <  4:91420/138‑1 (MQ=255)                    gCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGCGCGGCTTAAGCCAGCCGGGGTTTTTTGTTGAAAAGGTGTAATgtgcgtg                                                                                                           >  1:241152/1‑132 (MQ=255)                                    ggAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGCGCGGCTTCAGCCAGCGCGGGGGTTTTGATGAAAAAGTCTAATGTCATATTTAAAAATATTAAAAt                                                                                           >  8:347270/1‑138 (MQ=255)                                         acCAAAAACATCATCATGAACCCCGAGATAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAgg                                                                                     <  8:341744/139‑1 (MQ=255)                                         acCAAAAACATCATCAGTAACCCCGATAGAAAAGAGTAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAgg                                                                                     <  5:347017/139‑1 (MQ=255)                                          cctaaaaCATCATCAGTAACCCCGAGACAAAAGACGAATGAAACAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATCGACTAGTAAAAAATATTCAAATACTAGga                                                                                    <  8:182692/136‑1 (MQ=255)                                          ccAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGga                                                                                    <  7:13301/139‑1 (MQ=255)                                          ccAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCGCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGga                                                                                    <  1:53472/139‑1 (MQ=255)                                               aaCATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGaagaag                                                                               <  8:192562/139‑1 (MQ=255)                                                 catcatCACGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACACGAGCTGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTa                                                                             <  5:208816/139‑1 (MQ=255)                                                    catcaGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGCGCGGCGTCAGCCCGCGCGGGGTTTTTTTTGAAAAGGTCGTATGGACTATTAAAAAATATTCAAATAATATGAAAAAGt                                                                              >  6:341586/1‑135 (MQ=255)                                                    catcaGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCGTCAGCCAGCGCGGGGTTTTTGATGAAAAGGTGTAATGGACTAGTAAAAAATATTCAAATACTATGAAGAAGt                                                                              >  3:192501/1‑135 (MQ=255)                                                     atcaGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCGTCAGCCAGCTCGGGTTTTTTGTTGAAAAGGTCTTATGGTCTAGTAAAAAATATTAAAATACTAGGAAGAAGTAAGAt                                                                         >  8:384961/1‑139 (MQ=255)                                                         ggAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATtgtg                                                                     <  6:375682/139‑1 (MQ=255)                                                              cccGAGATAAAAGCGGAATGATTAAGCGAGAATGACCTCCCCCCCGAGCGGCTTCAGCCAGCGCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGatgaa                                                                <  2:157280/139‑1 (MQ=255)                                                              cccGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGatgaa                                                                <  1:301506/139‑1 (MQ=255)                                                               ccGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGatgaat                                                               <  3:310291/139‑1 (MQ=255)                                                                cGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGCGCGGCTTCCGCCCGCTCGGGGTTTTTGTTGAAAATGTCTAATGGGCTATTATAAAATATTGAAATAATAGGAATAAGTAATATTTTGatggat                                                               >  4:305913/1‑138 (MQ=255)                                                                  agagaAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGGTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATgag                                                            <  8:337769/139‑1 (MQ=255)                                                                     gaAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCGTCCGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATgaga                                                           >  1:41797/1‑137 (MQ=255)                                                                           agGAATGAATAAGCGAGATTGACCTCCCACCCGCCCCGCGTCACCCAGCGCGGGGTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAACTAAGATTGTGATGAATGAGATAGAaat                                                    >  8:388089/1‑138 (MQ=255)                                                                                  aataatCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATGAGATAGAAATAATCGtt                                             <  1:225157/138‑1 (MQ=255)                                                                                  aatAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATGAGATAGAAATAATCGtt                                             <  3:259589/138‑1 (MQ=255)                                                                                  aatAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCGCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATGAGATAGAAATAATCGtt                                             <  2:241152/138‑1 (MQ=255)                                                                                           gagtGACCTCCCACCCGAGCGGCTGCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATGAGATAGAAATAATCGTTACTTTCTTAc                                   <  5:29295/136‑1 (MQ=255)                                                                                                aCCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATGAGATAGAAATAATCGTTACTTTCTTACCTGtt                              >  8:388063/1‑139 (MQ=255)                                                                                                                 cTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATGAGATAGAAATAATCGTTACTTTCTTACCTGTTAAGAAAAGGGAAAATCt             <  1:210061/139‑1 (MQ=255)                                                                                                                   tCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGCCTAGTAAAAAATATTCAAATACTAGGAAGATGTAAGATTGTGATGAATGAGATAGAAATAATCGTTACTTTCTTAACTGTTAAGAAAAGGGAAAATCTGa           >  2:261046/1‑139 (MQ=255)                                                                                                                           cTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATGAGATAGAAATAATCGTTACTTTCTTACCTGTTAAGAAAAGGGAAAATCTGATTataata   >  2:187493/1‑139 (MQ=255)                                                                                                                            tCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATGAGATAGAAATAATCGTTACTTTCTTACCTGTTAAGAAAAGGGAAAATCTGATTATAATAg  <  5:31349/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                           GATACAAAGGAGCTGGAACGCACGACGCAGCCTCCGGAACACCAAAAACATCATCAGGAACCCCGAGAGAAAAGAGGAATGAATAAGCGAGATTGACCTCCCACCCGAGCGGCTTCAGCCAGCTCGGGTTTTTTGATGAAAAGGTCTAATGGACTAGTAAAAAATATTCAAATACTAGGAAGAAGTAAGATTGTGATGAATGAGATAGAAATAATCGTTACTTTCTTACCTGTTAAGAAAAGGGAAAATCTGATTATAATAG  >  CP006881/1252040‑1252301

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