Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,777,317	A→C	5.1%	K545T (AAG→ACG)	betA‑II →	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,777,317	0	A	C	5.1%	249.7 / 4.2	79	K545T (AAG→ACG)	betA‑II	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base C (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

TGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCA  >  NC_002947/5777170‑5777463                                                                                                                                                         |                                                                                                                                                    tGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAg                                                                                                                                                    <  2:219437/149‑1 (MQ=255)  gCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCCCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAg                                                                                                                                                    >  8:128670/1‑148 (MQ=255)      ggCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGtacta                                                                                                                                               <  7:128670/149‑1 (MQ=255)       gCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGtactac                                                                                                                                              <  8:254903/149‑1 (MQ=255)       gCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGtactac                                                                                                                                              <  7:73295/149‑1 (MQ=255)             aGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCg                                                                                                                                        >  3:126486/1‑149 (MQ=255)              tggCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCCCGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCCCTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCgc                                                                                                                                       >  7:196178/2‑148 (MQ=255)              gggCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCgg                                                                                                                                       >  3:60730/1‑149 (MQ=255)                      gcgTGGTCGATGCCTCGATCATGCCGCTGATCATGACCGGCAACCTCACCGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCATCCCCCGCCACGTAGTACTTGGCggtgctgcg                                                                                                                                <  6:60719/149‑8 (MQ=38)                       cgTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCAc                                                                                                                             <  2:346565/149‑1 (MQ=255)                        gTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCAc                                                                                                                             >  4:74527/1‑148 (MQ=255)                           gTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGt                                                                                                                         >  4:259464/1‑149 (MQ=255)                             cGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGa                                                                                                                       <  3:74527/149‑1 (MQ=255)                              gATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGAGGGGTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAg                                                                                                                          >  1:294157/1‑149 (MQ=255)                                tGCCTCGATCATGCCGATGATCATCACCGGAAAGCTCAACGCCACCACGATCATGATCGCCGATAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCTCAAGCCACTGCCGCGGA‑GCACCGCCAAGTACTACGTTGCGGGGG‑ATGCACCAGTGAAAg                                                                                                                    <  4:353720/149‑1 (MQ=255)                                   cTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGc                                                                                                                  >  7:61448/1‑148 (MQ=255)                                      gATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGc                                                                                                              >  1:49967/1‑149 (MQ=255)                                      gATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGc                                                                                                              <  5:12591/149‑1 (MQ=255)                                         gttGCCGGTGATGATCAGCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCgc                                                                                                           >  8:116354/3‑148 (MQ=255)                                          aTGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGa                                                                                                           <  7:171154/148‑1 (MQ=255)                                          aTGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGAt                                                                                                          <  1:125442/149‑1 (MQ=255)                                              cGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCgt                                                                                                      >  3:143255/1‑149 (MQ=255)                                              cGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCgt                                                                                                      <  8:46196/149‑1 (MQ=255)                                                   atcatcACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGC‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCgtgaa                                                                                                   <  4:312182/147‑1 (MQ=255)                                                   atcatcACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCgtgaa                                                                                                   <  5:357573/147‑1 (MQ=255)                                                     catcaCCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCgtgaagtg                                                                                                >  7:261743/1‑148 (MQ=255)                                                      atcaCCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCgtgaagtga                                                                                               <  8:61448/148‑1 (MQ=255)                                                       tcaCCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCgt                                                                                                      <  6:214164/140‑1 (MQ=255)                                                       tcaCCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCgt                                                                                                      >  5:214164/1‑140 (MQ=255)                                                          ccGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAg                                                                                          >  5:213878/1‑149 (MQ=255)                                                              cAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCt                                                                                      <  3:259464/149‑1 (MQ=255)                                                                  cTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAg                                                                                   >  3:217826/1‑148 (MQ=255)                                                                  cTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAg                                                                                   <  4:217826/148‑1 (MQ=255)                                                                   tCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGaaa                                                                                                                     >  7:4470/1‑113 (MQ=255)                                                                   tCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGaaa                                                                                                                     <  8:4470/113‑1 (MQ=255)                                                                   tCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGc                                                                                                                  <  1:125823/116‑1 (MQ=255)                                                                   tCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGc                                                                                                                  >  2:125823/1‑116 (MQ=255)                                                                     aaCGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGAACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCt                                                                               >  5:255519/1‑149 (MQ=255)                                                                           accacGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑At                                                                                                                                 >  3:185063/1‑93 (MQ=255)                                                                           accacGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑At                                                                                                                                 <  4:185063/93‑1 (MQ=255)                                                                            ccacGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAg                                                                        >  2:327382/1‑149 (MQ=255)                                                                             cacGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAgg                                                                       >  3:269189/1‑149 (MQ=255)                                                                                gATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGc                                                                    >  1:33207/1‑149 (MQ=255)                                                                                       aTCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCcgctgc                                                              <  4:143255/148‑4 (MQ=255)                                                                                             gagaAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGAGGCCGCGCAGCGGccc                                                       <  6:255519/149‑1 (MQ=255)                                                                                                  ga‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTc                                                  <  4:60730/149‑1 (MQ=255)                                                                                                  ga‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTc                                                  <  6:213878/149‑1 (MQ=255)                                                                                                   a‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCa                                                 >  5:166502/1‑149 (MQ=255)                                                                                                        tttCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCAc                                                >  1:360144/1‑149 (MQ=255)                                                                                                           cGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGt                                              <  4:269189/148‑1 (MQ=255)                                                                                                            ggACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGttt                                            >  7:342337/1‑149 (MQ=255)                                                                                                             gACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTg                                           >  6:80643/1‑149 (MQ=255)                                                                                                              aCAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCAAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGc                                          >  6:210474/1‑149 (MQ=255)                                                                                                               cAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCt                                         >  1:338603/1‑149 (MQ=255)                                                                                                                 aGATCCGTGTTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCa                                       <  8:342337/149‑1 (MQ=255)                                                                                                                  gATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAg                                      <  1:327382/149‑1 (MQ=255)                                                                                                                     ccGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCg                                   >  1:240041/1‑149 (MQ=255)                                                                                                                       gTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCggg                                 >  6:285893/1‑149 (MQ=255)                                                                                                                       gTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCggg                                 <  2:240041/149‑1 (MQ=255)                                                                                                                        tGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCggg                                 >  1:189941/1‑148 (MQ=255)                                                                                                                         ggTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCggg                                 >  2:34932/1‑147 (MQ=255)                                                                                                                          gTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAg                              <  2:49967/149‑1 (MQ=255)                                                                                                                              cAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAgcgc                           >  7:49585/1‑148 (MQ=255)                                                                                                                                   cACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGTCAGCGCaaaaaa                     <  1:34932/149‑1 (MQ=255)                                                                                                                                    aCTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAt                    >  5:345956/1‑149 (MQ=255)                                                                                                                                    aCTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAt                    <  4:126486/149‑1 (MQ=255)                                                                                                                                    aCTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAt                    >  6:229467/1‑149 (MQ=255)                                                                                                                                        ccGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtgt                  <  2:294157/147‑1 (MQ=255)                                                                                                                                        ccGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgc                <  2:189941/149‑1 (MQ=255)                                                                                                                                          gcgcA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgcg               >  7:57888/1‑148 (MQ=255)                                                                                                                                           cgcA‑GCACCGCCAAGTACTACGTGGCGGGGGAATCCACAAGTCAAAGGCAAGCCGAGGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGAAGCGCAAAAAAtgtgag               >  1:38158/1‑146 (MQ=255)                                                                                                                                              a‑gcaCCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCa          <  6:166502/149‑1 (MQ=255)                                                                                                                                              a‑gcaCCGCCAACACCTACGTGGCGGGGG‑ATGCACCAGTTAATGGCAAGGCGATGCGTCAAGGTAATCAGGGCTGAGGCCTTGATCAGGGGCCGGGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCa          <  8:196178/149‑1 (MQ=255)                                                                                                                                                gcaCCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAgg                                     >  4:115378/1‑121 (MQ=255)                                                                                                                                                gcaCCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAgg                                     <  3:115378/121‑1 (MQ=255)                                                                                                                                                  cccGCCACGTAGTACTTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAAc        >  8:9717/2‑148 (MQ=255)                                                                                                                                                   ccGCCACGTAGTACTTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACg       >  3:266566/1‑148 (MQ=255)                                                                                                                                                    cGCCACGTAGTACTTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGaa     <  7:9717/149‑1 (MQ=255)                                                                                                                                                       cAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAatca  >  7:5681/1‑149 (MQ=255)                                                                                                                                                         |                                                                                                                                                    TGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCA  >  NC_002947/5777170‑5777463

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