Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	670,055	A→G	intergenic (+88/‑94)	unknown → / → unknown	tRNA‑Pro‑TGG/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	670,055	0	A	G	100.0%	110.5 / NA	40	intergenic (+88/‑94)	unknown/unknown	tRNA‑Pro‑TGG/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/22);  total (18/22)

ATATGGTGGCTATAGCTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATAATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCTTAACAC  >  1/669808‑670205                                                                                                                                                                                                                                                        |                                                                                                                                                       atatGGTGGCTATAGCTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGAtt                                                                                                                                                     >  1:52309/1‑251 (MQ=255) atatGGTGGCTATAGCTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGAtt                                                                                                                                                     >  1:353134/1‑251 (MQ=255)     ggtggCTATAGCTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTaa                                                                                                                                                 >  1:275741/1‑251 (MQ=255)       tggCTATAGCTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGa                                                                                                                                               <  1:669715/251‑1 (MQ=255)        ggCTATAGCTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGAt                                                                                                                                              >  1:616093/1‑251 (MQ=255)          cTATAGCTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGAt                                                                                                                                              >  1:706993/1‑249 (MQ=255)          cTATAGCTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTg                                                                                                                                            >  1:504088/1‑251 (MQ=255)               gCTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAAc                                                                                                                                         >  1:412907/1‑249 (MQ=255)                cTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTa                                                                                                                                      >  1:87781/1‑251 (MQ=255)                     ttGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAAt                                                                                                                                    <  1:145799/248‑1 (MQ=255)                     ttGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAgg                                                                                                                                 <  1:285628/251‑1 (MQ=255)                     ttGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAgg                                                                                                                                 <  1:738451/251‑1 (MQ=255)                             gCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGa                                                                                                                          >  1:224232/1‑250 (MQ=255)                              ccccGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGa                                                                                                                          <  1:334893/249‑1 (MQ=255)                              ccccGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAg                                                                                                                        <  1:530161/251‑1 (MQ=255)                                   gATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCg                                                                                                                   >  1:226679/1‑251 (MQ=255)                                                 ttGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTttatttatatta                                                                                                      <  1:218491/250‑1 (MQ=255)                                                         ggTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                              <  1:342624/250‑1 (MQ=255)                                                                  cccATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAgg                                                                                                                                 >  1:419958/1‑206 (MQ=255)                                                                   ccATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAAt                                                                                                                                    >  1:56472/1‑202 (MQ=255)                                                                      ttAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTc                                                                                <  1:870935/251‑1 (MQ=255)                                                                          ccACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGaa                                                                                                                         <  1:859397/206‑1 (MQ=255)                                                                             ccccAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                              >  1:348975/1‑230 (MQ=255)                                                                              cccAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAAc                                                                                                                                         <  1:576738/186‑1 (MQ=255)                                                                              cccAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTg                                                                        <  1:744409/251‑1 (MQ=255)                                                                              cccAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTg                                                                        <  1:769996/251‑1 (MQ=255)                                                                               ccAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGaa                                                                                                                         <  1:255893/201‑1 (MQ=255)                                                                                cAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAg                                                                      <  1:605021/251‑1 (MQ=255)                                                                                                cAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTttga                                                      >  1:219767/1‑251 (MQ=255)                                                                                                     tGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTttgaattga                                                 >  1:773776/1‑251 (MQ=255)                                                                                                      ggTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCtt                                                           <  1:397822/240‑1 (MQ=255)                                                                                                             cAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTg                                         >  1:470007/1‑251 (MQ=255)                                                                                                                 tGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCa                                     <  1:813826/251‑1 (MQ=255)                                                                                                                              gTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGtttattt                                                                                                            >  1:203759/1‑167 (MQ=255)                                                                                                                                  gTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAAt                                                                                                                                    <  1:820150/139‑1 (MQ=255)                                                                                                                                  gTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACg                     >  1:161940/1‑250 (MQ=255)                                                                                                                                              aaTCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTttatttatatta                                                                                                      <  1:805384/157‑1 (MQ=255)                                                                                                                                               aTCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGattaatta                                                                                    <  1:318866/174‑1 (MQ=255)                                                                                                                                                    aTCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTCTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCTTAacac  <  1:584630/251‑1 (MQ=255)                                                                                                                                                                      aattctGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTttg                                                       <  1:603230/180‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                       ATATGGTGGCTATAGCTCAGTTGGTAGAGCCCCGGATTGTGATTCCGGTTGTCGTGGGTTCAAGTCCCATTAGCCACCCCAATCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGTAGGTTCAAATCCTATCTCGCCGACCAATTCTAATTCTGGGTTTTCTCAGAATTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATAATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCTTAACAC  >  1/669808‑670205

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