Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,777,321	C→G	5.3%	Y546* (TAC→TAG)	betA‑II →	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,777,321	0	C	G	5.3%	241.9 / 2.2	74	Y546* (TAC→TAG)	betA‑II	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (38/32);  new base G (2/2);  total (40/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 = 9.76e-01

GGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCG  >  NC_002947/5777175‑5777469                                                                                                                                                        |                                                                                                                                                      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                                                                                                            <  8:46196/149‑1 (MQ=255)                                         cGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCgt                                                                                                            >  3:143255/1‑149 (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                                                                                                                                       <  4:185063/93‑1 (MQ=255)                                                                      accacGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑At                                                                                                                                       >  3:185063/1‑93 (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                          <  4:126486/149‑1 (MQ=255)                                                                                                                               aCTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAt                          >  6:229467/1‑149 (MQ=255)                                                                                                                               aCTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAt                          >  5:345956/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                                           <  3:115378/121‑1 (MQ=255)                                                                                                                                           gcaCCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAgg                                           >  4:115378/1‑121 (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)                                                                                                                                                      tactacGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCt    >  2:229822/1‑149 (MQ=255)                                                                                                                                                      tactacGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCt    >  4:338429/1‑149 (MQ=255)                                                                                                                                                        ctacGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCg  >  7:63110/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                      GGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGA‑‑‑‑TTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCA‑GCACCGCCAAGTACTACGTGGCGGGGG‑ATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCG  >  NC_002947/5777175‑5777469

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