Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	196,868	A→G	N108S (AAC→AGC)	rseP →	intramembrane zinc metalloprotease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	196,868	0	A	G	100.0%	228.8 / NA	61	N108S (AAC→AGC)	rseP	intramembrane zinc metalloprotease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/31);  total (30/31)

GATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGCACCAGGTACGGAACTAAAAGCC  >  NC_000913/196734‑197004                                                                                                                                       |                                                                                                                                         gATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCa                                                                                                                                     >  2:1092345/1‑140 (MQ=255)      cGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCAttttt                                                                                                                                >  2:405136/1‑140 (MQ=255)       gTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTg                                                                                                                               <  2:691869/140‑1 (MQ=255)         tGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGATCTCCGCCGTCATGCTTTCTATAATAAATCTGTCGGCCAGCGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTgat                                                                                                                             <  2:857899/140‑3 (MQ=255)          gggCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTa                                                                                                                            <  2:113302/140‑1 (MQ=255)            gCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTc                                                                                                                                      >  1:706927/1‑128 (MQ=255)            gCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTc                                                                                                                                      <  2:706927/128‑1 (MQ=255)            gCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATc                                                                                                                          >  1:444107/1‑140 (MQ=255)             cGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCt                                                                                                                         <  2:790745/140‑1 (MQ=255)                         aaTGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTActggc                                                                                                             <  2:1027889/140‑1 (MQ=255)                          aTGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTActggct                                                                                                            <  2:516026/140‑1 (MQ=255)                           tGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTActggctg                                                                                                           <  2:444107/140‑1 (MQ=255)                             cTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGt                                                                                                         >  2:264849/1‑140 (MQ=255)                              tGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGtt                                                                                                        >  1:478347/1‑140 (MQ=255)                               ggATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGttt                                                                                                       >  1:8432/1‑140 (MQ=255)                                gATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCTTACTGGCTGGtttt                                                                                                      <  1:772073/140‑1 (MQ=255)                                        cgcAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTg                                                                                              >  1:839694/1‑140 (MQ=255)                                               ccGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGGCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGcc                                                                                       <  1:375265/140‑1 (MQ=255)                                               ccGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGcc                                                                                       <  2:8432/140‑1 (MQ=255)                                                cGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCt                                                                                      <  1:1165454/140‑1 (MQ=255)                                                  gTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTgg                                                                                    >  1:106690/1‑140 (MQ=255)                                                   tCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCACAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGc                                                                                   >  1:20583/1‑140 (MQ=255)                                                     gTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGt                                                                                 <  1:74328/140‑1 (MQ=255)                                                     gTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGt                                                                                 <  2:106690/140‑1 (MQ=255)                                                      ttCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTa                                                                                <  2:844661/140‑1 (MQ=255)                                                       tCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTAc                                                                               >  2:1184512/1‑140 (MQ=255)                                                                 cGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCggtggt                                                                     >  1:116037/1‑140 (MQ=255)                                                                 cGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCggtggt                                                                     >  1:1075794/1‑140 (MQ=255)                                                                 cGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCACAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCggtggt                                                                     <  2:20583/140‑1 (MQ=255)                                                                      ccaTGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCg                                                                >  2:427469/1‑140 (MQ=255)                                                                      ccaTGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCg                                                                >  2:656577/1‑140 (MQ=255)                                                                         tGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGaaa                                                             >  2:973166/1‑140 (MQ=255)                                                                         tGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGaaa                                                             >  2:531049/1‑140 (MQ=255)                                                                           ccTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATa                                                           >  2:497179/1‑140 (MQ=255)                                                                           ccTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATa                                                           >  2:489078/1‑140 (MQ=255)                                                                              tCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATagca                                                        <  1:264849/140‑1 (MQ=255)                                                                                aataatAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATagcagc                                                      >  2:1003380/1‑140 (MQ=255)                                                                                aataatAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATagca                                                        <  2:559034/138‑1 (MQ=255)                                                                                aataatAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATagca                                                        >  1:559034/1‑138 (MQ=255)                                                                                 ataataAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGcc                                                     <  1:1003380/140‑1 (MQ=255)                                                                                       aaTCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCg                                               <  2:428849/140‑1 (MQ=255)                                                                                            gTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGc                                          >  2:571320/1‑140 (MQ=255)                                                                                             tCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCt                                         >  2:60537/1‑140 (MQ=255)                                                                                             tCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCt                                         >  2:67216/1‑140 (MQ=255)                                                                                             tCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCt                                         >  2:463661/1‑140 (MQ=255)                                                                                             tCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCt                                         >  1:569378/1‑140 (MQ=255)                                                                                              cGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTg                                        <  1:343737/140‑1 (MQ=255)                                                                                                gCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCg                                      >  1:1147557/1‑140 (MQ=255)                                                                                                gCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCg                                      <  2:142889/140‑1 (MQ=255)                                                                                                     cGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGc                                 <  1:67216/140‑1 (MQ=255)                                                                                                      gAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGca                                <  1:497179/140‑1 (MQ=255)                                                                                                       aGCGGCGATTATTGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGcac                               >  1:698912/1‑140 (MQ=255)                                                                                                               ttattGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGCa                       <  1:489078/140‑1 (MQ=255)                                                                                                               ttattGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGCa                       <  1:427469/140‑1 (MQ=255)                                                                                                               ttattGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGCa                       <  2:1147557/140‑1 (MQ=255)                                                                                                                 attGCCGCAGGTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCGCAAATTGCAcc                     <  2:478347/140‑1 (MQ=255)                                                                                                                          ggTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGc                        <  1:155319/128‑1 (MQ=255)                                                                                                                          ggTCCGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGc                        >  2:155319/1‑128 (MQ=255)                                                                                                                             ccGGTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGCACCAGGTACGGAACt         >  1:27038/1‑140 (MQ=255)                                                                                                                               ggTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGCACCAGGTACGGAACTaa       <  1:405136/140‑1 (MQ=255)                                                                                                                               ggTTGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTACTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGCACCAGGTACGGAACTaa       <  2:1075794/140‑1 (MQ=255)                                                                                                                                 ttGCAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGCACCAGGTACGGAACTaaaa     <  1:463661/140‑1 (MQ=255)                                                                                                                                    cAAGCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGCTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGCACCAGGTACGGAACTAAAAGcc  <  2:839694/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         GATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAATTCGATAGCTGCGGAAGCACAAATTGCACCAGGTACGGAACTAAAAGCC  >  NC_000913/196734‑197004

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