Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,116,243	T→A	9.9%	L151Q (CTG→CAG)	lpd →	dihydrolipoyl dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,116,243	0	T	A	9.9%	161.7 / 5.0	58	L151Q (CTG→CAG)	lpd	dihydrolipoyl dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (23/29);  new base A (3/3);  total (26/32)
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

CTGTTTCGCAAGAACAAGGTCGAGTGGATCAAAGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCTGGCGCCGCCTGGGTAGCCAG  >  NC_002947/6116095‑6116388                                                                                                                                                     |                                                                                                                                                   cTGTTTCGCAAGAACAAGGTCGAGTGGATCAAAGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCctct                                                                                                                                                    >  1:206900/1‑149 (MQ=255) cTGTTTCGCAAGAACAAGGTCGAGTGGATCAAAGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCACTGAGCCCACCCctct                                                                                                                                                    >  5:71192/1‑149 (MQ=255) cTGTTTCGCAAGAACAAGGTCGAGTGGATCAAAGGATGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGAGTTGCAGGCCAAGGACAACGTCATTGCCACCGCCTCTGAGCCCACCCctct                                                                                                                                                    >  5:40414/1‑149 (MQ=255) cTGATTCGAAAGACCAAGCGCGAGTGCATCAAAGGCTGCGGCCGCCTGGCTGGCATCGGCCAGGTCGTCGTCAAGGCCGAAGACGGAAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCGGAGCCCACCCctct                                                                                                                                                    <  6:257326/149‑1 (MQ=255)  tGTTTCGCAAGAACAAGGTCGAGTGGATCAAAGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTg                                                                                                                                                   >  7:162378/1‑149 (MQ=255)   gTTTCGCAAGAACAAGGTCGAGTGGATCAAAGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGc                                                                                                                                                  >  7:173974/1‑149 (MQ=255)        gCAAGAACAAGGTCGAGTGGATCAAAGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCAACCCTCTGCC‑AGgt                                                                                                                                            >  5:159317/1‑149 (MQ=255)               cAAGGTCGAGTGGATCAAAGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCg                                                                                                                                     <  5:306050/149‑1 (MQ=255)                      gAGTGGATCAAAGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCTAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACaa                                                                                                                              >  6:89637/1‑149 (MQ=255)                           gATCAAAGGCTGGGGCCGACTGGATGGCATATGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGAGAAACCCCGTTGCAGCCCAAGGACATCCTCATTACCACCGGCACTGAGCCCACCACTCTGCA‑AGGTGTGACCGACGACagaacg                                                                                                                          >  6:269847/1‑143 (MQ=255)                             tCAAAGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTTCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAgcg                                                                                                                        >  6:173555/1‑148 (MQ=255)                               aaaGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGcat                                                                                                                     >  6:217082/1‑149 (MQ=255)                               aaaGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGcat                                                                                                                     >  6:221852/1‑149 (MQ=255)                                aaGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAAc                                                                                                                             <  7:87137/140‑1 (MQ=255)                                aaGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAAc                                                                                                                             >  8:87137/1‑140 (MQ=255)                                aaGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGcatc                                                                                                                    <  1:2453/149‑1 (MQ=255)                                aaGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGcat                                                                                                                     >  2:27962/1‑148 (MQ=255)                                   gCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCatcatc                                                                                                                 <  2:206900/149‑1 (MQ=255)                                   gCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCatcatc                                                                                                                 <  8:134190/149‑1 (MQ=255)                                                  aTGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGgcgcg                                                                                                  <  1:27962/149‑1 (MQ=255)                                                     gCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCt                                                                                                <  6:113367/148‑1 (MQ=255)                                                      cATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGt                                                                                              >  2:27554/1‑149 (MQ=255)                                                          ggCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCt                                                                                          <  2:156630/149‑1 (MQ=255)                                                                 tCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCaa                                                                                   <  8:173974/149‑1 (MQ=255)                                                                      gTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACACGACCGGCGCGCTGTCGCTGCCGCAAGTGc                                                                               <  1:27554/148‑1 (MQ=255)                                                                         aaGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCaa                                                                           >  7:199312/1‑149 (MQ=255)                                                                         aaGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCaa                                                                           >  5:22030/1‑149 (MQ=255)                                                                           ggCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAg                                                                          >  6:114164/1‑148 (MQ=255)                                                                           ggCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAg                                                                          >  8:98951/1‑148 (MQ=255)                                                                               gAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCt                                                                     >  8:33567/1‑149 (MQ=255)                                                                               gAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCAGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAATTGCCCAAGCACCt                                                                     >  2:298440/1‑149 (MQ=255)                                                                               gAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCAGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCt                                                                     >  2:298438/1‑149 (MQ=255)                                                                               gAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCAGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCt                                                                     >  2:303056/1‑149 (MQ=255)                                                                                   aCGGTAGCGAAACCGCGTTGCAGGCTAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGtt                                                                 <  5:89637/149‑1 (MQ=255)                                                                                      gTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTc                                                              <  6:168925/149‑1 (MQ=255)                                                                                       tAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGtt                                                                 <  1:88309/145‑1 (MQ=255)                                                                                       tAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGtt                                                                 >  2:88309/1‑145 (MQ=255)                                                                                          cGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCg                                                          >  8:135145/1‑149 (MQ=255)                                                                                                 gcgTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTgg                                                   >  6:250610/1‑149 (MQ=255)                                                                                                 gcgTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCCACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTgg                                                   <  5:114164/149‑1 (MQ=255)                                                                                                       cAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGAt                                             >  4:150620/1‑149 (MQ=255)                                                                                                       cAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCAGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGAt                                             <  1:303056/149‑1 (MQ=255)                                                                                                       cAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCAGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGAt                                             <  1:298438/149‑1 (MQ=255)                                                                                                       cAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCAGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCCGTGCTGGCGTGAt                                             <  1:298440/149‑1 (MQ=255)                                                                                                        aGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATc                                            <  8:199312/149‑1 (MQ=255)                                                                                                            cAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTc                                        <  5:221852/149‑1 (MQ=255)                                                                                                            cAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTc                                        <  5:217082/149‑1 (MQ=255)                                                                                                              aGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCtt                                      <  3:280549/149‑1 (MQ=255)                                                                                                              aGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCtt                                      <  7:33567/149‑1 (MQ=255)                                                                                                                 aCATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAg                                   <  3:150620/149‑1 (MQ=255)                                                                                                                   aTCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCt                                 >  2:109089/1‑149 (MQ=255)                                                                                                                      gTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTggg                              <  5:173555/149‑1 (MQ=255)                                                                                                                      gTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTgg                               <  6:159317/148‑1 (MQ=255)                                                                                                                            gCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCgg                         <  8:162378/148‑1 (MQ=255)                                                                                                                             ccACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTc                           >  4:11632/1‑145 (MQ=255)                                                                                                                             ccACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTc                           <  3:11632/145‑1 (MQ=255)                                                                                                                             ccACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTc                       <  2:250451/149‑1 (MQ=255)                                                                                                                              cACGGGCTCTGAGCCCTACCCTCTGCC‑TGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCATGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCt                      <  6:40414/149‑1 (MQ=255)                                                                                                                               aCCGGCTCTGAGCCCACCCCTCTGCCGAGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCt                      >  1:177685/1‑149 (MQ=255)                                                                                                                                         aGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCTGGCGCCGCCt            <  5:250610/148‑1 (MQ=255)                                                                                                                                           cccACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCTGGCGCCGCCTggg         <  2:177685/149‑1 (MQ=255)                                                                                                                                             cACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCTGGCGCCGCCTGGGTa       <  2:153215/149‑1 (MQ=255)                                                                                                                                                  ctctGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCTGGCGCCGCCTGGGTAGCCAg  <  5:217844/149‑1 (MQ=255)                                                                                                                                                  ctctGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCTGGCGCCGCCTGGGTAGCCAg  >  1:99466/1‑149 (MQ=255)                                                                                                                                                  ctctGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCTGGCGCCGCCTGGGTAGCCAg  >  3:125682/1‑149 (MQ=255)                                                                                                                                                  ctctGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCTGGCGCCGCCTGGGTAGCCAg  <  2:99466/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CTGTTTCGCAAGAACAAGGTCGAGTGGATCAAAGGCTGGGGCCGCCTGGATGGCATCGGCAAGGTCGTGGTCAAGGCCGAAGACGGTAGCGAAACCGCGTTGCAGGCCAAGGACATCGTCATTGCCACCGGCTCTGAGCCCACCCCTCTGCC‑AGGTGTGACCGTCGACAACCAGCGCATCATCGACTCGACCGGCGCGCTGTCGCTGCCGCAAGTGCCCAAGCACCTGGTTGTCATCGGTGCTGGCGTGATCGGTCTTGAGCTGGGTTCGGTCTGGCGCCGCCTGGGTAGCCAG  >  NC_002947/6116095‑6116388

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: ‑