Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,458,751	A→G	6.9%	intergenic (+43/‑225)	pcaJ → / → trkH‑I	3‑oxoadipate CoA‑transferase subunit A/potassium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,458,751	0	A	G	6.9%	167.5 / 2.2	58	intergenic (+43/‑225)	pcaJ/trkH‑I	3‑oxoadipate CoA‑transferase subunit A/potassium transporter
Reads supporting (aligned to +/- strand):  ref base A (25/29);  new base G (2/2);  total (27/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 = 9.97e-01

GCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGCCGGC  >  NC_002947/4458604‑4458897                                                                                                                                                    |                                                                                                                                                     gCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCAcaatcaa                                                                                                                                                     >  4:108009/1‑149 (MQ=255) gCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCAcaatca                                                                                                                                                      <  3:315916/148‑1 (MQ=255)   cGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCAcaatcaa                                                                                                                                                     >  4:77231/1‑147 (MQ=255)   cGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATa                                                                                                                                                   <  8:57772/149‑1 (MQ=255)          gAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCg                                                                                                                                             >  5:73618/1‑148 (MQ=255)            aGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGcac                                                                                                                                          <  1:93018/149‑1 (MQ=255)             gTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGcac                                                                                                                                          >  2:132477/1‑148 (MQ=255)             gTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGcac                                                                                                                                          >  1:58006/1‑148 (MQ=255)                    cTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTATGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCAGCACAAGcg‑ccg                                                                                                                                 >  3:139920/1‑149 (MQ=255)                     tGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACTATTGATTGTGCCGcaga                                                                                                                                         >  6:222696/1‑139 (MQ=255)                     tGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACTATTGATTGTGCCGcaga                                                                                                                                         <  5:222696/141‑3 (MQ=255)                        aGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGa                                                                                                                              <  4:9670/148‑1 (MQ=255)                             tGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAgg                                                                                                                         >  8:225672/1‑148 (MQ=255)                             tGAAAGTGGTCGAAATCTGCGCGGACATCGACTTAGACGAGCTGCAGAAAATCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAgg                                                                                                                         >  1:50344/1‑148 (MQ=255)                                   tGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCAcg                                                                                                                 <  3:77231/149‑1 (MQ=255)                                     gTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCAcgcg                                                                                                               <  6:73618/149‑1 (MQ=255)                                        gAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTg                                                                                                            <  3:262379/149‑1 (MQ=255)                                         aaaTCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCATAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGc                                                                                                           <  1:274289/149‑1 (MQ=255)                                           aTCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCtt                                                                                                         <  1:329548/149‑1 (MQ=255)                                                cgcgGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTgaga                                                                                                     <  3:117307/148‑1 (MQ=255)                                                cgcgGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTga                                                                                                       <  1:132477/146‑1 (MQ=255)                                                cgcgGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTga                                                                                                       <  6:185961/146‑1 (MQ=255)                                                cgcgGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTga                                                                                                       <  2:58006/146‑1 (MQ=255)                                                         cGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTa                                                                                            >  2:187669/1‑148 (MQ=255)                                                          gACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACa                                                                                          <  1:187669/149‑1 (MQ=255)                                                           aCTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAg                                                                                         <  3:78904/149‑1 (MQ=255)                                                            cTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGc                                                                                        >  2:129677/1‑149 (MQ=255)                                                                  cGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTc                                                                                  <  5:357059/149‑1 (MQ=255)                                                                  cGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTc                                                                                  <  4:139920/149‑1 (MQ=255)                                                                  cGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTc                                                                                  <  4:369212/149‑1 (MQ=255)                                                                  cGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTc                                                                                  <  7:114376/149‑1 (MQ=255)                                                                          acttgatCAGCGGCACGCCACTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGa                                                                          >  2:375621/7‑149 (MQ=255)                                                                            aaaCTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATg                                                                        <  7:225672/149‑1 (MQ=255)                                                                              aCTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATgg                                                                       >  3:268447/1‑148 (MQ=255)                                                                               cTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGaa                                                                     <  4:99369/149‑1 (MQ=255)                                                                                tCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAAc                                                                    >  7:127263/1‑149 (MQ=255)                                                                                         tGCCGCTGATCAAGTGATCCAGGAAGACGTAGGGCGTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGgaga                                                           <  5:236512/149‑1 (MQ=255)                                                                                                gATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAAc                                                    >  3:35253/1‑149 (MQ=255)                                                                                                 aTCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAAcc                                                   >  6:90859/1‑149 (MQ=255)                                                                                                   cAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACtctc                                                                <  8:287056/134‑1 (MQ=255)                                                                                                   cAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACtctc                                                                >  7:287056/1‑134 (MQ=255)                                                                                                     aGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAAt                                               >  4:339633/1‑149 (MQ=255)                                                                                                     aGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAAt                                               <  2:50344/149‑1 (MQ=255)                                                                                                        gATCCAGGAAGAGTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGAGTACGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTAAACCCAAt                                               >  1:109222/1‑147 (MQ=255)                                                                                                        gATCCAGGAAGAATTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGTACGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGATTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGa                                             >  1:114040/1‑149 (MQ=255)                                                                                                             aGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATc                                       >  6:112180/1‑149 (MQ=255)                                                                                                                 aGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGc                                   <  1:129677/149‑1 (MQ=255)                                                                                                                   aCTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCgg                                 >  3:228088/1‑149 (MQ=255)                                                                                                                       aGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCaa                                                     >  5:252115/1‑125 (MQ=255)                                                                                                                       aGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCaa                                                     <  6:252115/125‑1 (MQ=255)                                                                                                                       aGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAg                              <  3:108009/148‑1 (MQ=255)                                                                                                                          gCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAg                          <  3:339633/149‑1 (MQ=255)                                                                                                                           cTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTAt                                        >  5:8165/1‑134 (MQ=255)                                                                                                                           cTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTAt                                        <  6:8165/134‑1 (MQ=255)                                                                                                                            tttCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGa                         >  4:26327/1‑148 (MQ=255)                                                                                                                            tttCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGa                         >  2:179375/1‑148 (MQ=255)                                                                                                                                 aTGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGGAATCACGCGGTGCTTGAGCCTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcg                    >  1:283685/1‑149 (MQ=255)                                                                                                                                 aTGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcgg                   >  7:351475/1‑149 (MQ=255)                                                                                                                                  tGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcggc                  <  4:116386/149‑1 (MQ=255)                                                                                                                                    cttgtGCGGCACTATTGATTGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcgg                   >  7:104472/5‑146 (MQ=255)                                                                                                                                    cttgtGCGGCACTATTGATTGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcgg                   <  8:104472/142‑1 (MQ=255)                                                                                                                                               caatcaatAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGcc     <  1:375621/149‑1 (MQ=255)                                                                                                                                                  tcaatAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGCCGGc  >  3:340754/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGCCGGC  >  NC_002947/4458604‑4458897

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