Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,435,815	A→T	8.2%	L50Q (CTG→CAG)	PP_1256 ←	alpha‑ketoglutarate semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,435,815	0	A	T	8.2%	170.4 / 5.7	61	L50Q (CTG→CAG)	PP_1256	alpha‑ketoglutarate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base T (2/3);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

TGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCA  >  NC_002947/1435669‑1435962                                                                                                                                                    |                                                                                                                                                    tGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAgc                                                                                                                                                   <  6:36846/149‑1 (MQ=255) tGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAgc                                                                                                                                                   >  8:56138/1‑149 (MQ=255) tGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAgc                                                                                                                                                   >  8:51755/1‑149 (MQ=255)        cTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAgg                                                                                                                                            >  6:208325/1‑149 (MQ=255)             gCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTgg                                                                                                                                       <  7:165604/149‑1 (MQ=255)             gCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTgg                                                                                                                                       <  6:259878/149‑1 (MQ=255)              cAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCTGCGCGCTGGCCTggg                                                                                                                                      <  1:27442/149‑1 (MQ=255)                ggCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGcg                                                                                                                                    <  4:283107/149‑2 (MQ=255)                     gTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACgtg                                                                                                                               >  7:194060/1‑149 (MQ=255)                     gTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACgtg                                                                                                                               >  5:126369/1‑149 (MQ=255)                       ttCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGtt                                                                                                                             >  1:22365/1‑149 (MQ=255)                            tCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCt                                                                                                                        >  2:148047/1‑149 (MQ=255)                                          aTCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGcg                                                                                                                                    <  4:311431/123‑2 (MQ=255)                                          aTCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGcg                                                                                                                                    >  3:311431/1‑122 (MQ=255)                                             aGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGc                                                                                                       <  2:22365/149‑1 (MQ=255)                                               ggCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCAc                                                                                                                                  >  3:201182/1‑120 (MQ=255)                                               ggCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCAc                                                                                                                                  <  4:201182/120‑1 (MQ=255)                                                         aGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGt                                                                                           <  1:148047/149‑1 (MQ=255)                                                             ccTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAAc                                                                                             >  3:74144/1‑143 (MQ=255)                                                             ccTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGACGCCGAGGTAGGCCGGCTCCAAc                                                                                             <  4:74144/143‑1 (MQ=255)                                                              cTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGAc                                                                                      >  2:12380/1‑149 (MQ=255)                                                              cTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGAc                                                                                      >  2:12382/1‑149 (MQ=255)                                                                      gtgtgGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGtt                                                                              >  6:203430/1‑149 (MQ=255)                                                                       tgtgGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCTCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGtt                                                                              <  8:194060/148‑1 (MQ=255)                                                                       tgtgGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGGGGTAGGCCGGCTCCAACGTCTGACCGGTGGtt                                                                              <  5:208325/148‑1 (MQ=255)                                                                       tgtgGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCgg                                                                                   <  6:266795/143‑1 (MQ=255)                                                                       tgtgGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGtt                                                                              <  1:122632/148‑1 (MQ=255)                                                                           gCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGt                                                                          <  6:107877/148‑1 (MQ=255)                                                                           gCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGt                                                                          <  5:203430/148‑1 (MQ=255)                                                                                  cTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCggtg                                                                                 >  4:118256/1‑134 (MQ=255)                                                                                  cTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCggtg                                                                                 <  3:118256/134‑1 (MQ=255)                                                                                   ttCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATGAGAATGCCGCACAGCCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCa                                                                  >  1:72343/1‑149 (MQ=255)                                                                                      aaGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGa                                                              >  5:84546/1‑149 (MQ=255)                                                                                      aaGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGa                                                              >  5:27183/1‑149 (MQ=255)                                                                                       aGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGAt                                                             <  3:200773/149‑1 (MQ=255)                                                                                        ggAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATg                                                            <  7:51755/149‑1 (MQ=255)                                                                                        ggAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATg                                                            <  7:56138/149‑1 (MQ=255)                                                                                                 cTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCgg                                                   >  3:198070/1‑149 (MQ=255)                                                                                                   cGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCt                                                 >  7:261350/1‑149 (MQ=255)                                                                                                     ttgttCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGc                                               <  1:12380/149‑1 (MQ=255)                                                                                                     ttgttCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGc                                               <  1:12382/149‑1 (MQ=255)                                                                                                       gttCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCg                                             >  1:193470/1‑149 (MQ=255)                                                                                                        ttCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCg                                             <  6:27183/148‑1 (MQ=255)                                                                                                           gAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTa                                         >  6:144430/1‑149 (MQ=255)                                                                                                            aGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTAc                                        >  6:183084/1‑149 (MQ=255)                                                                                                               gATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGa                                                                                       <  5:75200/99‑1 (MQ=255)                                                                                                               gATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGa                                                                                       >  6:75200/1‑99 (MQ=255)                                                                                                               gATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCgg                                     >  3:181900/1‑149 (MQ=255)                                                                                                                aTGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTATGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGc                                    >  7:181736/1‑149 (MQ=255)                                                                                                                       gcgGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTgg                              <  6:126369/148‑1 (MQ=255)                                                                                                                       gcgGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGc                             <  5:144430/149‑1 (MQ=255)                                                                                                                       gcgGGAGGCAT‑CGAATGCCGCCCCGGCCAGCTCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACTGCTGACCGGTGGTTGGGTAGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCctcttg                              <  6:271638/148‑3 (MQ=255)                                                                                                                                 t‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagcag                   <  8:165974/149‑1 (MQ=255)                                                                                                                                 t‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagcag                   >  1:32597/1‑149 (MQ=255)                                                                                                                                     aaTGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGtt                <  4:198070/149‑1 (MQ=255)                                                                                                                                      aTGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTg               >  1:246893/1‑149 (MQ=255)                                                                                                                                       tGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGc              <  5:183084/149‑1 (MQ=255)                                                                                                                                          cgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCtg           >  8:61395/1‑149 (MQ=255)                                                                                                                                              cAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAg       >  7:47872/1‑149 (MQ=255)                                                                                                                                                   cAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCa  <  2:209373/149‑1 (MQ=255)                                                                                                                                                   cAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTCCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCa  >  1:209373/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCAT‑CGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCA  >  NC_002947/1435669‑1435962

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