Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	583,398	G→A	8.3%	A422T (GCC→ACC)	selB →	selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	583,398	0	G	A	8.3%	176.3 / 3.3	60	A422T (GCC→ACC)	selB	selenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base A (2/3);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCTGCGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGCCTTG  >  NC_002947/583254‑583544                                                                                                                                                  |                                                                                                                                                   gTGATGCGCGATGCAGCGGATCTGGAGCAGGCCCTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTC‑CGGGTGGTGGCCAc                                                                                                                                                >  6:68997/1‑149 (MQ=255)   gctgcGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCAccc                                                                                                                                              >  1:182189/1‑149 (MQ=255)    ctgcGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCAccc                                                                                                                                              <  1:158662/148‑1 (MQ=255)      gcgcGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGCCACCACCCgca                                                                                                                                           <  7:277500/149‑3 (MQ=255)        gcgATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGCCACCACCCgcaca                                                                                                                                         <  7:4617/149‑5 (MQ=255)         cgATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAgg                                                                                                                                        <  1:270760/149‑1 (MQ=255)           aTGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGc                                                                                                                                      <  5:287440/149‑1 (MQ=255)               aGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAAc                                                                                                                                  <  6:267662/149‑1 (MQ=255)                   gATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAActgct                                                                                                                              >  1:85401/1‑149 (MQ=255)                        ggAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCg                                                                                                                         >  8:33337/1‑149 (MQ=255)                        ggAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCg                                                                                                                         >  2:97709/1‑149 (MQ=255)                         gAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACTTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTc                                                                                                                          <  8:60288/147‑1 (MQ=255)                         gAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTc                                                                                                                          <  2:182189/147‑1 (MQ=255)                         gAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTc                                                                                                                          <  7:28677/147‑1 (MQ=255)                         gAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTc                                                                                                                          <  2:85401/147‑1 (MQ=255)                              ggCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCCCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACg                                                                                                                   <  1:97709/149‑1 (MQ=255)                                     ccAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATg                                                                                                            >  6:205714/1‑149 (MQ=255)                                            tGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCCGGGTGGCCACCACCCgcac                                                                                                                                          >  2:252476/1‑110 (MQ=255)                                            tGCTCGCCAGTGCACCGCGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTGCCGGGTGGCCACCACCCgcac                                                                                                                                          <  1:252476/113‑4 (MQ=255)                                               tCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCt                                                                                                   >  3:136223/1‑148 (MQ=255)                                               tCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCt                                                                                                   >  4:118087/1‑148 (MQ=255)                                                   cAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGACACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAag                                                                                               >  7:180769/1‑148 (MQ=255)                                                        cACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGcagcag                                                                                         <  1:86489/149‑1 (MQ=255)                                                         aCCGGGCGGCATTAGCCCGCAGCGACTCGANCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAgg                                                                                        >  1:23694/1‑149 (MQ=255)                                                         aCCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAgg                                                                                        >  1:18873/1‑149 (MQ=255)                                                           cGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTc                                                                                                  >  2:53072/1‑137 (MQ=255)                                                           cGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTc                                                                                                  <  1:53072/137‑1 (MQ=255)                                                              ctttgcatAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTg                                                                                   <  7:71732/142‑1 (MQ=255)                                                                ggcATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGa                                                                                 <  4:136223/149‑1 (MQ=255)                                                                ggcATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGGGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGa                                                                                 <  7:133813/149‑1 (MQ=255)                                                                    ttAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAAt                                                                                                             <  4:83833/117‑1 (MQ=255)                                                                    ttAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAAt                                                                                                             >  3:83833/1‑117 (MQ=255)                                                                       gCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGtt                                                                           <  8:109414/148‑1 (MQ=255)                                                                         ccGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGgc                                                                        <  3:118087/149‑1 (MQ=255)                                                                           gCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGgcgc                                                                      <  5:37127/149‑1 (MQ=255)                                                                                     gAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            <  2:18873/149‑1 (MQ=255)                                                                                     gAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            <  5:205714/149‑1 (MQ=255)                                                                                     gAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            <  5:15017/149‑1 (MQ=255)                                                                                     gAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            <  2:23694/149‑1 (MQ=255)                                                                                       gCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGGTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGc                                                          <  4:162510/149‑1 (MQ=255)                                                                                            aGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            <  5:184613/142‑1 (MQ=255)                                                                                            aGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGa                                                            >  6:184613/1‑142 (MQ=255)                                                                                            aGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAagag                                                     >  7:6181/1‑149 (MQ=255)                                                                                                   ccGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTc                                                                <  7:3557/131‑1 (MQ=255)                                                                                                   ccGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTc                                                                >  8:3557/1‑131 (MQ=255)                                                                                                    cGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGa                                                                                 >  5:142519/1‑113 (MQ=255)                                                                                                    cGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGa                                                                                 <  6:142519/113‑1 (MQ=255)                                                                                                             gAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGa                                                <  5:94514/137‑1 (MQ=255)                                                                                                             gAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGa                                                >  6:94514/1‑137 (MQ=255)                                                                                                             gAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGccc                                    <  7:147377/149‑1 (MQ=255)                                                                                                             gAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGcc                                     >  8:147377/1‑148 (MQ=255)                                                                                                                  cTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGAcc                               >  4:14078/1‑149 (MQ=255)                                                                                                                        gTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATc                         <  7:33337/149‑1 (MQ=255)                                                                                                                         ttACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCAc                                                              >  2:289100/1‑111 (MQ=255)                                                                                                                         ttACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCAc                                                              <  1:289100/111‑1 (MQ=255)                                                                                                                           aCCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCt                      >  5:125951/1‑149 (MQ=255)                                                                                                                                    tgtg‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCt             >  8:175125/1‑149 (MQ=255)                                                                                                                                    cctgccGGGTGGCCACCACCC‑GCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCt             >  2:241951/6‑149 (MQ=255)                                                                                                                                      tg‑CGGGTGGTGGCTACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAATAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCttt           <  5:68997/149‑1 (MQ=255)                                                                                                                                            gtggtgGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGc      >  5:97284/1‑149 (MQ=255)                                                                                                                                             tggtggCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGcc     >  6:292183/1‑149 (MQ=255)                                                                                                                                                tggCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGCCTTg  <  8:6181/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GTGCTGCGCGATGCAGCGGATCTGGAGCAGGCACTGCCAGCACTGCTCGCCAGTGCACCGGGCGGCATTAGCCCGCAGCGACTCGAGCGGCAGTTCAACCGCCTGCGGGAAACCTGGCAGTTACCGGATGATGTG‑CGGGTGGTGGCCACCCGGCAGGGCCAACTGCTGTTCGCCAACGCTCAATGGCAGGCGCTCAAGCAGCAGGTCCTGGAGCAGTTGGCGCGCTTTCACGAGCAAGAGCCCGACCAGCTTGGCCCTGACCGTGATCGCTTGCGGCGCTTTGTAGCCTTG  >  NC_002947/583254‑583544

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