Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	1,530,173	A→C	100%	*399S (TAA→TCA)  noncoding (62/244 nt)	ftsZ → PP_mr19 →	Cell division protein FtsZ Xoo8 ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,530,173	0	A	C	91.7%	110.1 / 3.8	38	*399S (TAA→TCA)  noncoding (62/244 nt)	ftsZ PP_mr19	Cell division protein FtsZ Xoo8 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base C (11/23);  total (12/26)
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.94e-02
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGTTGACACCG  >  AE015451/1529959‑1530392                                                                                                                                                                                                                       |                                                                                                                                                                                                                            gACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCATCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGtata                                                                                                                                                                                                        >  1:286854/1‑236 (MQ=255)   cAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATaaa                                                                                                                                                                                                      <  1:206704/236‑1 (MQ=255)      tACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGt                                                                                                                                                                                                   <  1:277370/236‑1 (MQ=255)                  tGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATtca                                                                                                                                                                                       <  1:333162/236‑1 (MQ=255)                     ccAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGa                                                                                                                                                                                                                        >  2:234545/1‑200 (MQ=255)                      cagcagGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAAtt                                                                                                                                                                                                                     >  2:452936/1‑202 (MQ=255)                      cagcagGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAAtt                                                                                                                                                                                                                     <  1:452940/202‑1 (MQ=255)                                                  cgcCGGTGCGCCAGGAGCAGCCAGCGGTCACCTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGGCAGGCTGATTCATGGAATTTATCAGGGGt                                                                                                                                                                                                           <  2:93372/184‑1 (MQ=255)                                                    ccGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                  <  1:88340/175‑1 (MQ=255)                                                     cGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAg                                                                                                                                                          >  2:243205/1‑230 (MQ=255)                                                     cGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTttg                                                                                                                                                    >  1:294448/1‑236 (MQ=255)                                                       gTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTAtt                                                                                                                                                                                         <  1:7658/197‑1 (MQ=255)                                                          cgcCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAgggg                                                                                                                                                                                                            >  2:488970/1‑175 (MQ=255)                                                             cAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGt                                                                                                                                                                                                   <  2:198853/181‑1 (MQ=255)                                                              aGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                  <  2:178521/165‑1 (MQ=255)                                                                                       tGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGt                                                                                                                                                                                                   <  2:18889/155‑1 (MQ=255)                                                                                         aCCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCa                                                                                                                                                           <  1:4257/193‑1 (MQ=255)                                                                                          ccTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCaa                                                                                                                <  1:257212/235‑1 (MQ=255)                                                                                                  gcCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                  <  2:456144/129‑1 (MQ=255)                                                                                                  gcCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                  >  1:456148/1‑129 (MQ=255)                                                                                                      gACCGTCATGCGCAACCAGGCGCATGCAGATGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCTATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGTGATAAAg                                                                                                                                                                                                     >  2:181800/1‑138 (MQ=255)                                                                                                                      cAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTatcatc                                                                                                                                                               <  1:133805/160‑1 (MQ=255)                                                                                                                         gcgcATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTgaagac                                                                                  >  1:439781/1‑233 (MQ=255)                                                                                                                          cgcATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATCTCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTgaa                                                                                     <  2:227693/230‑1 (MQ=255)                                                                                                                                               ccGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATAGCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCaa                                                                                                                <  2:7786/182‑1 (MQ=255)                                                                                                                                                            cccTCAGGATGACCTGTACTGCCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATGTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATAcc                                                                              <  2:409622/203‑1 (MQ=255)                                                                                                                                                              cTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCggg                                             <  2:286850/234‑1 (MQ=255)                                                                                                                                                                   gATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGgaga                                         <  1:512290/233‑1 (MQ=255)                                                                                                                                                                           ggACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAAtt                                                                                                                                                                                                                     <  2:15141/53‑1 (MQ=255)                                                                                                                                                                             aCTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGAc                            <  1:317384/236‑1 (MQ=255)                                                                                                                                                                                          cccGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGAtt                                                                                                      >  1:366319/1‑149 (MQ=255)                                                                                                                                                                                          cccGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGc                                                                                                                      >  2:25070/1‑133 (MQ=255)                                                                                                                                                                                              gCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                  <  1:276508/37‑1 (MQ=255)                                                                                                                                                                                                 ttCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGtt         >  1:472592/1‑235 (MQ=255)                                                                                                                                                                                                          cttcAGGCTGATTCATGGAATTGATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTGCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGgtg                                                              <  2:7361/171‑1 (MQ=255)                                                                                                                                                                                                          cgtcAGTCTGTTTTATGGAATTTATCAGGGGTATAAAGGTTATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCaa                                                                                                                <  2:8242/123‑1 (MQ=255)                                                                                                                                                                                                                   gATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGAc                            <  2:294444/198‑1 (MQ=255)                                                                                                                                                                                                                   gATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGTTGACACCg  <  1:462014/224‑1 (MQ=255)                                                                                                                                                                                                                       |                                                                                                                                                                                                                            GACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGTTGACACCG  >  AE015451/1529959‑1530392

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