Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	649,629	C→A	8.3%	A17S (GCC→TCC)	accC ←	acetyl‑CoA carboxylase biotin carboxylase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	649,629	0	C	A	8.3%	207.8 / 6.6	69	A17S (GCC→TCC)	accC	acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base C (35/28);  new base A (4/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

AGGTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGCC  >  NC_002947/649487‑649772                                                                                                                                                 |                                                                                                                                                aGGTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGACTTGATACACAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTc                                                                                                                                           >  5:283274/1‑149 (MQ=255)  ggTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTcc                                                                                                                                          >  7:180868/1‑149 (MQ=255)   gTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTccc                                                                                                                                         >  5:355332/1‑149 (MQ=255)   gTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTccc                                                                                                                                         >  4:137050/1‑149 (MQ=255)   gTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTccc                                                                                                                                         >  3:264126/1‑149 (MQ=255)   gTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTcc                                                                                                                                          <  5:348832/148‑1 (MQ=255)   gTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTcc                                                                                                                                          >  6:348832/1‑148 (MQ=255)   gTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGAAGGCACGCAGGATCCGCAGGGCAATTTccc                                                                                                                                         >  2:90014/1‑149 (MQ=255)   gTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGAACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCAGCAGGGAAATTTccc                                                                                                                                         >  6:206957/1‑149 (MQ=255)     aGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCcg                                                                                                                                       >  6:297401/1‑149 (MQ=255)     aGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATAAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCcg                                                                                                                                       >  8:198547/1‑149 (MQ=255)       gAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCcgcg                                                                                                                                     >  2:327320/1‑149 (MQ=255)            cTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTgg                                                                                                                                <  1:90014/149‑1 (MQ=255)               tgACGATGCAGGGCCGATGCAGACCGACTCCGCTGCCATGCCCATGTGGATCAGTTGACCGTCGGCATTGGAGTGGACGCCACGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTg                                                                                                                                 <  1:311214/146‑1 (MQ=255)                    aTGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAgg                                                                                                                        <  3:50270/149‑1 (MQ=255)                        aGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCGGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACtt                                                                                                                    >  5:178194/1‑149 (MQ=255)                        aGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACtt                                                                                                                    >  7:116139/1‑149 (MQ=255)                            ccGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTc                                                                                                                <  6:283274/149‑1 (MQ=255)                                         aCTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACACAGCTCTTTGCAGGCACGCAGGATCAGCAGGGCAATTACCCCGCGGTTGGCGATCAGGACGTTTTCGAGCTTCCCAGa                                                                                                    >  6:42815/1‑148 (MQ=255)                                             gTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACTCAGGATCCGAAGGGCAATTTACCCGCGGTTGGCTATCAGGACTTTTTCGAGCTTCCCAGACATc                                                                                                <  5:206957/148‑1 (MQ=255)                                                  ccAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCTGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTgg                                                                                           <  6:178194/148‑1 (MQ=255)                                                  ccAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGtt                                                                                             <  8:116139/146‑1 (MQ=255)                                                  ccAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGtt                                                                                             <  8:268422/146‑1 (MQ=255)                                                  ccAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTg                                                                                            <  8:180868/147‑1 (MQ=255)                                                      gCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCtctc                                                                                      <  1:327320/149‑1 (MQ=255)                                                        ccAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCcg                                                                                    >  4:262204/1‑149 (MQ=255)                                                        ccAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCcg                                                                                    >  4:36222/1‑149 (MQ=255)                                                           ggTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGa                                                                                 <  4:279635/149‑1 (MQ=255)                                                                aTCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCa                                                                             >  4:285380/1‑148 (MQ=255)                                                                aTCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GAGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCa                                                                             >  3:329640/1‑148 (MQ=255)                                                                 tCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCaaa                                                                           <  7:2823/149‑1 (MQ=255)                                                                 tCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCaaa                                                                           <  3:285380/149‑1 (MQ=255)                                                                    gTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTACCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGa                                                                        >  7:333004/1‑149 (MQ=255)                                                                       cACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATg                                                                      <  3:262204/148‑1 (MQ=255)                                                                          ggTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTg                                                                   >  4:123598/1‑148 (MQ=255)                                                                                       gTGCACG‑‑GCGACGGGCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTc                                                     <  4:264126/149‑1 (MQ=255)                                                                                            cG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGaa                                                                 >  1:19686/1‑132 (MQ=255)                                                                                            cG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGaa                                                                 <  2:19686/132‑1 (MQ=255)                                                                                             g‑‑gCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTc                                               <  6:355332/149‑1 (MQ=255)                                                                                                 cGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCaa                                             >  5:153087/1‑149 (MQ=255)                                                                                                      gTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGc                                        >  7:32871/1‑149 (MQ=255)                                                                                                         ttGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTgg                                     >  5:5451/1‑149 (MQ=255)                                                                                                         ttGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTgg                                     >  1:27365/1‑149 (MQ=255)                                                                                                           gATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTg                                      <  5:297401/146‑1 (MQ=255)                                                                                                            aTACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCg                                  >  5:42594/1‑149 (MQ=255)                                                                                                                ccAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGc                                                             <  5:98318/118‑1 (MQ=255)                                                                                                                ccAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGc                                                             >  6:98318/1‑118 (MQ=255)                                                                                                                ccAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGtct                               >  7:365934/1‑148 (MQ=255)                                                                                                                ccAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAAACGGCTGGCCGtcg                               >  7:336769/1‑147 (MQ=255)                                                                                                                ccAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGAGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGtct                               >  1:197186/1‑148 (MQ=255)                                                                                                                 cAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGtcttc                             >  1:86947/1‑149 (MQ=255)                                                                                                                   gCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTa                           >  1:104346/1‑149 (MQ=255)                                                                                                                      ctTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCa                        <  8:32871/149‑1 (MQ=255)                                                                                                                         tGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGTTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGa                     >  3:243767/1‑149 (MQ=255)                                                                                                                           cAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATg                   <  4:329640/149‑1 (MQ=255)                                                                                                                           cAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATg                   <  3:137050/149‑1 (MQ=255)                                                                                                                             ggCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGc                 <  2:197186/149‑1 (MQ=255)                                                                                                                                aCGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTc              >  5:50980/1‑149 (MQ=255)                                                                                                                                 cGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCg             <  3:36222/149‑1 (MQ=255)                                                                                                                                  gCAGTATCCGCCTGTCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCGCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGa            <  8:336769/149‑1 (MQ=255)                                                                                                                                      caaccGCGGGGAAATTGCCCTGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGac        >  7:335255/4‑149 (MQ=255)                                                                                                                                        tCCTCAGGGCAATTTCCCCGCTGTTGGCGTTCAGGACTTTTTCGATCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAAGAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGacac      <  3:288618/149‑1 (MQ=255)                                                                                                                                         ccGCGGGGAAATTGCCCTGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTc                                                     <  7:68740/101‑1 (MQ=255)                                                                                                                                         ccGCGGGGAAATTGCCCTGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTc                                                     >  8:68740/1‑101 (MQ=255)                                                                                                                                         ccGCGGGGAAATTGCCCTGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAg                       >  2:288232/1‑131 (MQ=255)                                                                                                                                         ccGCGGGGAAATTGCCCTGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAg                       <  1:288232/131‑1 (MQ=255)                                                                                                                                           gCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGc   <  6:153087/149‑1 (MQ=255)                                                                                                                                           gCAGGGCAATTTACCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGCCATCGTTGGCTCTCCGCGATTCAACCGATGGTGAACAGCGGATGTGCGACCTAAACCGGATGGCCGTCTTCTACCAGGATGGCGTCGCTGACAcctc   >  5:4743/1‑147 (MQ=255)                                                                                                                                             aGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGcc  >  1:317492/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                AGGTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGCC  >  NC_002947/649487‑649772

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