Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	445,967	Δ1 bp	intergenic (‑29145/‑51496)	unknown ← / → unknown	SSU rRNA ## 16S rRNA, small subunit ribosomal RNA/tRNA‑SeC(p)‑TCA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	445,967	0	A	.	100.0%	151.2 / NA	36	intergenic (‑29145/‑51496)	unknown/unknown	SSU rRNA ## 16S rRNA, small subunit ribosomal RNA/tRNA‑SeC(p)‑TCA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (25/11);  total (25/11)

TTATTAAACTCTTCTATTATTGGAGTATCTTCCTGATAAGAAAATAATGCACTACATAATTTCAAAAATAATCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTCATATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGAATTATCACTAACAAACACTCTCAATGCATCCATATTCTTTAAAATTATTCTAACTCCAGCTAAATCTCCCATTCGAGATAATTTCATATC  >  1/445738‑446151                                                                                                                                                                                                                                      |                                                                                                                                                                                         ttattaAACTCTTCTATTATTGGAGTATCTTCCTGATAAGAAAATAATGCACTACATAATTTCAAAAATAATCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAccc                                                                                                                                                                                   >  1:158594/1‑236 (MQ=255)         ctctTCTATTATTGGAGTATCTTCCTGATAAGAAAATAATGCACTACATAATTTCAAAAATAATCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCtatat                                                                                                                                                            >  1:849639/1‑251 (MQ=255)            ttctATTATTGGAGTATCTTCCTGATAAGAAAATAATGCACTACATAATTTCAAAAATAATCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAAt                                                                                                                                                         <  1:670784/251‑1 (MQ=255)                attattGGAGTATCTTCCTGATAAGAAAATAATGCACTACATAATTTCAAAAATAATCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAg                                                                                                                                                                     >  1:809215/1‑235 (MQ=255)                attattGGAGTATCTTCCTGATAAGAAAATAATGCACTACATAATTTCAAAAATAATCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATg                                                                                                                                                     >  1:691269/1‑251 (MQ=255)                           aTCTTCCTGATAAGAAAATAATGCACTACATAATTTCAAAAATAATCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTtttattta                                                                                                                                          >  1:449292/1‑251 (MQ=255)                                                 gCACTACATAATTTCAAAAATAATCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAc                                                                                                                                                                                     >  1:804842/1‑186 (MQ=255)                                                                 aaaaTAATCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAc                                                                                                                                                                                     <  1:761651/170‑1 (MQ=255)                                                                        tCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAAt                                                                                                                                                         >  1:259790/1‑191 (MQ=255)                                                                        tCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTg                                                                                             >  1:369311/1‑251 (MQ=255)                                                                         cTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAccc                                                                                                                                                                                   >  1:528749/1‑164 (MQ=255)                                                                         cTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAcc                                                                                                                                                                                    <  1:883818/163‑1 (MQ=255)                                                                          ttATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGaa                                                                                           <  1:550500/251‑1 (MQ=255)                                                                                 cTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAc                                                                                                                                                                                     <  1:474039/154‑1 (MQ=255)                                                                                 cTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAc                                                                                                                                                                                     >  1:574511/1‑154 (MQ=255)                                                                                 cTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAc                                                                                                                                                                                     <  1:638106/154‑1 (MQ=255)                                                                                  tcatcaCCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAAt                                                                                                                                                         >  1:627471/1‑181 (MQ=255)                                                                                       aCCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAAt                                                                                                                                                         >  1:278026/1‑176 (MQ=255)                                                                                            gTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAc                                                                                                                                                                                     >  1:318208/1‑143 (MQ=255)                                                                                                 aCTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGAATTATCACTaacaaac                                                                            <  1:553941/243‑1 (MQ=255)                                                                                                                         gagCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAAt                                                                                                                                                         >  1:611182/1‑142 (MQ=255)                                                                                                                           gCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAgg                                                                                                                                                                       >  1:538993/1‑126 (MQ=255)                                                                                                                                       aaCTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGAATTATCACTAACAAacac                                                                          <  1:771604/207‑1 (MQ=255)                                                                                                                                        aCTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGAATTATCACTAACAAACACTCTCAATGCATCCATATTCTTTAAAATTATTCTAACTCCAGCTaa                             >  1:474368/1‑251 (MQ=255)                                                                                                                                           aTGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCa                                                                                                                                                       >  1:653292/1‑126 (MQ=255)                                                                                                                                                 gCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGAATTATCACTaacaaac                                                                            <  1:231597/195‑1 (MQ=255)                                                                                                                                                   ccACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAAt                                                                                                                                                         >  1:350718/1‑116 (MQ=255)                                                                                                                                                      cTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAc                                                                                                                                                                                     >  1:492216/1‑85 (MQ=255)                                                                                                                                                           gttgtAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATAc                                                                                                                                                                                     >  1:306805/1‑80 (MQ=255)                                                                                                                                                                                                 ttATAATTTTTATCTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGAATTATCACTaacaaac                                                                            >  1:90044/1‑147 (MQ=255)                                                                                                                                                                                                              cTTCAAATACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAAt                                                                                                                                                         >  1:424787/1‑57 (MQ=255)                                                                                                                                                                                                                     tACTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAAt                                                                                                                                                         >  1:686427/1‑50 (MQ=255)                                                                                                                                                                                                                      aCTTGATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGAATTATCACTAACAAACACTCTCAATg                                                                  <  1:851198/136‑1 (MQ=255)                                                                                                                                                                                                                          gATGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGaa                                                                                           >  1:616444/1‑107 (MQ=255)                                                                                                                                                                                                                           aTGTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGAATTATCACTAACAAACACTCTCAATGCATCCATATTCTTTAAAATTATTCTAACTCCAGCTAAATCTCCCATTCGAGATAATTTCATATc  >  1:536718/1‑195 (MQ=255)                                                                                                                                                                                                                             gTACACTTC‑TATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGAATTATCACTAACAAACACTCTCAATGCATCCATATTCTTTAAAATTATTCTAACTCCAGCTAAATCTCCCATTCgaga              <  1:322025/181‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                         TTATTAAACTCTTCTATTATTGGAGTATCTTCCTGATAAGAAAATAATGCACTACATAATTTCAAAAATAATCTTATTCCCTCATCACCTTGTCCAACTTTAAATGATGTTTTCATCAAAGAGCCCAATATTTCAACTATGGTTGCCCACTCATGTTGTAACTGAGTTCTAATTTGCAACTCTAATTTTAGTTTATAATTTTTATCTTCAAATACTTGATGTACACTTCATATACCCATCTTTTTTAGGAGTTTCTATATAATCATGCGTTTTATTTACATAACAGAAGAAATTTTCCTTCCTAAGTTCACAGATTGTATTTGAATTATCACTAACAAACACTCTCAATGCATCCATATTCTTTAAAATTATTCTAACTCCAGCTAAATCTCCCATTCGAGATAATTTCATATC  >  1/445738‑446151

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