Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,621,173	A→T	10.7%	intergenic (‑37/+82)	U712_18670 ← / ← U712_18675	50S ribosomal protein L31/Transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,621,173	0	A	T	10.7%	123.4 / 3.7	45	intergenic (‑37/+82)	U712_18670/U712_18675	50S ribosomal protein L31/Transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base T (3/2);  total (23/22)
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

AACGCGTACCTCTTCTTTTACTGAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACCTTCTTGCAGATGATAGATTTGCCTGTTTCCATTCTTGATTGAG  >  CP006881/3621038‑3621302                                                                                                                                        |                                                                                                                                  aaCGCGTACCTCTTCTTTTACTGAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGaaa                                                                                                                                  <  3:231757/137‑1 (MQ=255)   cgcgTACCTCTTCTTTTACTGAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACaga                                                                                                                              <  1:100492/139‑1 (MQ=255)                  ttACTGAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGTATGAATTCCTGCTTTCATTGTATCCATCTCCTTCAGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAg                                                                                                               >  3:247404/1‑139 (MQ=255)                    aCTGAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGc                                                                                                              <  4:112539/138‑1 (MQ=255)                    aCTGAGCCTGTTTCAAATTCATTTCCGAAAGCGCATTTAACTGTTGCTTTTTTGAAATTATGATGAATTCCTCCTTACATTGAATCCATCTCCTTCCGCCCTGAATCCAGTGTGGAAACAGAGTTATTATAACGGAGCa                                                                                                             >  5:78260/1‑139 (MQ=255)                     cTGAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCa                                                                                                             <  5:125004/138‑1 (MQ=255)                       gAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGATACAGAGTTATTATTTCGGAGCATTg                                                                                                          >  7:271865/1‑139 (MQ=255)                       gAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGATACAGAGTTATTATTTCGGAGCATTg                                                                                                          >  4:24086/1‑139 (MQ=255)                       gAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTg                                                                                                          >  1:101526/1‑139 (MQ=255)                       gAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTg                                                                                                          >  6:148181/1‑139 (MQ=255)                              tttCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTCTCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGTAGCATTGTCCGaaa                                                                                                   >  1:267805/1‑139 (MQ=255)                              tttCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGaaa                                                                                                   <  4:75525/139‑1 (MQ=255)                              tttCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATATCGGAGCATTGTCCGaa                                                                                                    <  8:62861/138‑1 (MQ=255)                              tttCAAATTCATTTCCGCAAGCGCAGTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGaaa                                                                                                   >  6:17781/1‑139 (MQ=255)                              attCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGaa                                                                                                    <  4:23819/137‑1 (MQ=255)                                        aTTTCCGCAAGCTCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATgag                                                                                         <  3:163047/139‑1 (MQ=255)                                        aTTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATgag                                                                                         <  8:264272/139‑1 (MQ=255)                                        aTTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATga                                                                                          <  6:280680/138‑1 (MQ=255)                                               cAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATta                                                                                  <  2:101526/139‑1 (MQ=255)                                                  gcgcATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTatat                                                                               <  3:24600/139‑1 (MQ=255)                                                            cTGTTGCTTTTTTGAAATTAGGATGAATACCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTa                                                                     <  4:161932/139‑1 (MQ=255)                                                            cTGTTGCTGTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTa                                                                     <  6:152992/139‑1 (MQ=255)                                                                tGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGaa                                                                                       <  2:81580/117‑1 (MQ=255)                                                                tGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGaa                                                                                       >  1:81580/1‑117 (MQ=255)                                                                tGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTActtc                                                                 >  1:214349/1‑139 (MQ=255)                                                                tGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTActtc                                                                 >  1:28968/1‑139 (MQ=255)                                                                                ggATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAAt                                                  <  7:114998/138‑1 (MQ=255)                                                                                             cTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGATACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACCTTCTTGCa                                     <  3:24086/138‑1 (MQ=255)                                                                                               ttCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACCTTCTTGCAgat                                  >  4:133061/1‑139 (MQ=255)                                                                                                  aTTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCg                                                                                                      >  4:281626/1‑68 (MQ=255)                                                                                                  aTTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCg                                                                                                      <  3:281626/68‑1 (MQ=255)                                                                                                        tCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACATTCTTGCAGATGATAGAttt                         >  6:212511/1‑139 (MQ=255)                                                                                                        tCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCAGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACCTTCTTGCAGATGATAGAttt                         >  7:163065/1‑139 (MQ=255)                                                                                                            tctcCTTCCGCCCTGAATCCCGTGTGGATACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCaaa                                                        <  5:257346/104‑1 (MQ=255)                                                                                                            tctcCTTCCGCCCTGAATCCCGTGTGGATACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCaaa                                                        >  6:257346/1‑104 (MQ=255)                                                                                                            tctcCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGaa                                                                                                    >  8:208458/1‑60 (MQ=255)                                                                                                            tctcCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGaa                                                                                                    <  7:208458/60‑1 (MQ=255)                                                                                                            tctcCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCa                                                                             <  7:78027/83‑1 (MQ=255)                                                                                                            tctcCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCa                                                                             >  8:78027/1‑83 (MQ=255)                                                                                                             ctcCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACCTTCTTGCAGATGATAGATTTGCCTg                    >  6:10339/1‑139 (MQ=255)                                                                                                                cTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACCTTCTTGCAGATGATAGATTTGCGTGttt                 >  1:214366/1‑139 (MQ=255)                                                                                                                  tCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACCTTCTTGCAGATGATAGATTTGCCTGTTTc                >  4:188976/1‑138 (MQ=255)                                                                                                                       ccTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACCTTCTTGCAGATGATAGATTTGCCTGTTTCCATTCt          <  5:10339/139‑1 (MQ=255)                                                                                                                                cGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACCTTCTTGCAGATGATAGATTTGCCTGTTTCCATTCTTGATTgag  >  3:30205/1‑138 (MQ=255)                                                                                                                                 gtgtGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAAtt                                                 >  6:275006/1‑90 (MQ=255)                                                                                                                                 gtgtGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAAtt                                                 <  5:275006/90‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  AACGCGTACCTCTTCTTTTACTGAGCCTGTTTCAAATTCATTTCCGCAAGCGCATTTAACTGTTGCTTTTTTGAAATTAGGATGAATTCCTGCTTTCATTGTATCCATCTCCTTCCGCCCTGAATCCCGTGTGGAAACAGAGTTATTATTTCGGAGCATTGTCCGAAACGCATATGAGAAAATTATATCAGTGATTTACTTCTTTTGCAAATACAATTTACCTTCTTGCAGATGATAGATTTGCCTGTTTCCATTCTTGATTGAG  >  CP006881/3621038‑3621302

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