Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,226,865	G→T	8.1%	L11L (CTC→CTA)	PP_5623 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,226,865	0	G	T	8.1%	142.1 / 2.4	49	L11L (CTC→CTA)	PP_5623	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (22/22);  new base T (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

GGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGATTTATTTCTAGGAACCACAGTTTCC  >  NC_002947/4226721‑4227008                                                                                                                                                 |                                                                                                                                                ggCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCag                                                                                                                                                 <  5:180437/145‑1 (MQ=255)       aaaGAGGCGCCGCGACTCTTCAAGCGCATCCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGaa                                                                                                                                       <  1:286511/149‑1 (MQ=255)       aaaGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTAttt                                                                                                                                          <  7:319690/146‑1 (MQ=255)       aaaGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATAATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTAttt                                                                                                                                          >  8:319690/1‑146 (MQ=255)        aaGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGaac                                                                                                                                      >  3:76950/1‑148 (MQ=255)                 cgcgACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCttt                                                                                                                             <  3:100537/149‑1 (MQ=255)                    gACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTaaa                                                                                                                          <  5:336923/149‑1 (MQ=255)                          tCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTgg                                                                                                                     >  7:78003/1‑148 (MQ=255)                             aGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTc                                                                                                                                >  5:51692/1‑134 (MQ=255)                             aGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTc                                                                                                                                <  6:51692/134‑1 (MQ=255)                               cgcAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGacac                                                                                                                 >  3:361291/1‑147 (MQ=255)                                    ccTCGCTGTCCTTAGACAAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCAACCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCtt                                                                                                           >  2:112804/1‑148 (MQ=255)                                      tCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCa                                                                                                        >  3:236003/1‑149 (MQ=255)                                              cTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTGTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAAc                                                                                                <  4:85005/149‑1 (MQ=255)                                              cTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATATTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAAc                                                                                                <  6:27519/149‑1 (MQ=255)                                                 aGACCAATAACTCGACTGGTCAGACAATCCCTTATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGATTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACAGTTAACTg                                                                                              <  8:78003/148‑1 (MQ=255)                                                            tCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGc                                                                                  >  8:238024/1‑149 (MQ=255)                                                             cGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATATTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGcc                                                                                 <  4:6141/149‑1 (MQ=255)                                                                 tGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTGTCACCGTTAAGTGATCCACGATTGcgg                                                                                >  8:305333/1‑144 (MQ=255)                                                                       gacaATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTg                                                                       <  2:135293/149‑1 (MQ=255)                                                                        acaATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCATAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTg                                                                       >  8:126550/1‑148 (MQ=255)                                                                        acaATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCATAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTg                                                                       >  8:132047/1‑148 (MQ=255)                                                                        acaATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTg                                                                       >  5:85953/1‑148 (MQ=255)                                                                            tCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATaa                                                                  >  2:384358/1‑149 (MQ=255)                                                                                     gcgTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGa                                                         <  1:384358/149‑1 (MQ=255)                                                                                          tCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTc                                                    <  6:85953/149‑1 (MQ=255)                                                                                             aCCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCaa                                                  >  2:170831/1‑148 (MQ=255)                                                                                             aCCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAAt                                                 >  8:303969/1‑149 (MQ=255)                                                                                               ccATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCGGTTAGATTTTCAATc                                                >  8:315849/1‑148 (MQ=255)                                                                                                 atcatATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACCGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCAc                                             <  3:139573/149‑1 (MQ=255)                                                                                                    atatTTTTTCTGATTTTCATGTGTGAGTGGAATCAGCTGGGGCGAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTTGATTTTCAATCCACGc                                           <  4:76950/148‑1 (MQ=255)                                                                                                     tatTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGcc                                          >  4:245315/1‑148 (MQ=255)                                                                                                            tCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGAcc                                  >  7:14855/1‑149 (MQ=255)                                                                                                             cTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAATCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGAcc                                  <  3:185165/148‑1 (MQ=255)                                                                                                              tGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGt                                >  1:261718/1‑149 (MQ=255)                                                                                                              tGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGt                                >  1:189393/1‑149 (MQ=255)                                                                                                               gATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGtt                               >  3:84656/1‑149 (MQ=255)                                                                                                                 ttttCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTg                             <  1:112804/149‑1 (MQ=255)                                                                                                                    tCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGa                           <  2:261718/148‑1 (MQ=255)                                                                                                                     cATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGat                          <  1:170831/148‑1 (MQ=255)                                                                                                                      aTCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGattt                        <  3:245315/149‑1 (MQ=255)                                                                                                                             aGTGGAATCAGCTGGGGCCATAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGATTTATTTCTa                 <  7:132047/149‑1 (MQ=255)                                                                                                                             aGTGGAATCAGCTGGGGCCACAGTTTTTGAAGGCAGTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACACGTTAGATTTTCATTCGACGCCCACGGCCCGTTTGGATTTAGTTCTa                 >  6:345793/1‑149 (MQ=255)                                                                                                                              gTGGAATCAGCTGGGGCCATAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGATTTATTTCTa                 <  7:126550/148‑1 (MQ=255)                                                                                                                                  aaTCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGATTTATTTCTAGGAAc            >  6:170228/1‑149 (MQ=255)                                                                                                                                  aaTCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGATTTATTTCTAGGAAc            <  4:361291/149‑1 (MQ=255)                                                                                                                                       gCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGATTTATTTCTAGGAACCACAg       <  1:201540/149‑1 (MQ=255)                                                                                                                                            ggCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGATTTATTTCTAGGAACCACAGTTTcc  <  7:238024/149‑1 (MQ=255)                                                                                                                                              ccAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGATTTATTTCTAGGAACCACAGTTTcc  >  6:6665/1‑147 (MQ=255)                                                                                                                                              ccAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGATTTATTTCTAGGAACCACAGTTTcc  >  5:140810/1‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                                GGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTATCCACCCATCATATTTTTTCTGATTTTCATCTGTGAGTGGAATCAGCTGGGGCCAGAGTATTTGAAGGCATTCTTTAAATCTGGACACCAGCTTTCACCGTTAACTGATCCACGATTGCCGAAGCAATTGGATAACCCGTTAGATTTTCAATCCACGCCCACTGACCGTTTGGATTTATTTCTAGGAACCACAGTTTCC  >  NC_002947/4226721‑4227008

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