Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	654,408	T→C	intergenic (‑227/‑245)	unknown ← / → unknown	cell filamentation‑like protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	654,408	0	T	C	100.0%	114.5 / NA	38	intergenic (‑227/‑245)	unknown/unknown	cell filamentation‑like protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/24);  total (14/24)

TATAGCGTTAAGTGCATCGTCCTCTCGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGTTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAGTGCTCAAGAAACAAACAGTGAAAAGCAAATTTTCAACCAAAGCCCATAGGGGATTTAAT  >  1/654287‑654651                                                                                                                          |                                                                                                                                                                                                                                                    tgtaacgacaaGTGCATCGTCCTCTCGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGt                                                                                                                                                                                                                                   >  1:484379/10‑140 (MQ=255) tgtaacgacaaGTGCATCGTCCTCTCGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAg                                                                                                                    <  1:93709/242‑1 (MQ=255)       gacaaGTGCATCGTCCTCTCGAATATAATTAAGTAGCTCATTGAGCTGCTCTCCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATgg                                                                                                                                                                                                                         <  1:710457/141‑1 (MQ=255)            gTGCATCGTCCTCTCGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGt                                                                                                                                                                                                                                   >  1:283631/1‑129 (MQ=255)            gTGCATCGTCCTCTCGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATgg                                                                                                                                                                                                                         <  1:480309/139‑1 (MQ=255)            gTGCATCGTCCTCTCGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTcagc                                                                                                                                                                                                  <  1:664973/162‑1 (MQ=255)             tGCATCGTCCTCTCGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTa                                                                                                                                                                            >  1:564451/1‑183 (MQ=255)                        ctcGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTcagcag                                                                                                                                                                                                <  1:55586/152‑1 (MQ=255)                         tcGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAgg                                                                                            <  1:715071/251‑1 (MQ=255)                          cGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGc                                                                                                           >  1:47876/1‑235 (MQ=255)                                taATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAg                                                                                                                                                                                                                                      >  1:180351/1‑106 (MQ=255)                                taATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGt                                                                                                                                                                              <  1:540067/162‑1 (MQ=255)                                taATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTAc                                                                                                                                               <  1:294717/193‑1 (MQ=255)                                taATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTcc                                                                                                                                         <  1:87433/199‑1 (MQ=255)                                taATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCtt                                                                                                         >  1:452819/1‑231 (MQ=255)                                taATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAg                                                                                     <  1:720731/251‑1 (MQ=255)                                taATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGAACTCTTTAGCGTTATTTAATGGGATATAtgt                                                                                                                                                                                                                >  1:519194/1‑128 (MQ=255)                                taATTAAGTAGCTCATTGAGCTGCTCTTCG‑TTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATAtgt                                                                                                                                                                                                                <  1:72391/127‑1 (MQ=255)                                            tCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTc                                                                         <  1:181026/251‑1 (MQ=255)                                             cATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGc                                                                                                           <  1:809885/216‑1 (MQ=255)                                                 gAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATAtg                                                                                                                                                                                                                 <  1:15548/110‑1 (MQ=255)                                                 gAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCCGGAGGGGAAATTGGCATGTTATCTGCTAGaa                                                                                                                                                                  >  1:729029/1‑157 (MQ=255)                                                 gAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCAcc                                                                      >  1:872700/1‑249 (MQ=255)                                                        tctTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAGGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAg                                                             <  1:165483/251‑1 (MQ=255)                                                        tctTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAg                                                             <  1:290231/251‑1 (MQ=255)                                                        tctTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAg                                                             <  1:824369/251‑1 (MQ=255)                                                                          tttGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAGTGCTCAAGAAACAAACAg                                           <  1:897685/251‑1 (MQ=255)                                                                                      aaGTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAg                                                             <  1:561639/221‑1 (MQ=255)                                                                                        gTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAAtt                                                                                                                                                                                     <  1:845044/99‑1 (MQ=255)                                                                                        gTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCcaac                                                                                                                                      <  1:893223/146‑1 (MQ=255)                                                                                        gTTTTATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAGTGCTCAAGAAACAAACAGTGAAAAGCAAAttt                             >  1:295806/1‑251 (MQ=255)                                                                                           ttATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTa                                                                                                                                                                            >  1:780500/1‑105 (MQ=255)                                                                                           ttATATGACCGAATTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGcc                                                                                   >  1:517745/1‑194 (MQ=255)                                                                                                      aaTTTTTGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAAtt                                                                                                                                                                                     >  1:269484/1‑85 (MQ=255)                                                                                                          tttGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAGTGCTCAAGAAACAAACAGTGAAAAGCAAATTTTCAACCAAAGCCCATAgg           <  1:538432/251‑1 (MQ=255)                                                                                                           ttGCATTGCCGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAGTGCTCAAGAAACAAACAGTGAAAAGCAAATTTTCAACCAAAGCCCATAggg          <  1:41439/251‑1 (MQ=255)                                                                                                                   ccGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAGTGCTCAAGAAACAAACAGTGAAAAGCAAATTTTCAACCAAAGCCCATAGGGGATTTAAt  <  1:470850/251‑1 (MQ=255)                                                                                                                    cGCTTGCTGTGGACCTCTTTAGCGTTATTTAATGGGATATAtgt                                                                                                                                                                                                                >  1:862668/1‑44 (MQ=255)                                                                                                                          |                                                                                                                                                                                                                                                    TATAGCGTTAAGTGCATCGTCCTCTCGAATATAATTAAGTAGCTCATTGAGCTGCTCTTCGTTTTTTTTCGGCTTTGCCTGAAAAAAGTTTTATATGACCGAATTTTTGCATTGCCGCTTGTTGTGGACCTCTTTAGCGTTATTTAATGGGATATATGTTGTAAGATTTCAGCAGGAGGGGAAATTGGCATGTTATCTGCTAGAAAGGTCGGCACGCCCAGTACGTTTCCAACAACAAGCAAGATCCGCAGGTTTACCGCTTCGCTGCATTAAGGCTTATAGCCTTAATAATTCACCCTTTTAAAGTGCTCAAGAAACAAACAGTGAAAAGCAAATTTTCAACCAAAGCCCATAGGGGATTTAAT  >  1/654287‑654651

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