Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	649,638	C→T	8.6%	G14R (GGG→AGG)	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,638	0	C	T	8.6%	162.1 / 3.9	58	G14R (GGG→AGG)	accC	acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base T (3/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

AGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGCCAATATCGGC  >  NC_002947/649491‑649781                                                                                                                                                      |                                                                                                                                                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:36222/1‑149 (MQ=255)                                                    ccAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCcg                                                                                             >  4:262204/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                                                                                    <  3:285380/149‑1 (MQ=255)                                                             tCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCaaa                                                                                    <  7:2823/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                                              >  1:27365/1‑149 (MQ=255)                                                                                                     ttGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTgg                                              >  5:5451/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                            <  3:137050/149‑1 (MQ=255)                                                                                                                       cAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATg                            <  4:329640/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                                                              >  8:68740/1‑101 (MQ=255)                                                                                                                                     ccGCGGGGAAATTGCCCTGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTc                                                              <  7:68740/101‑1 (MQ=255)                                                                                                                                     ccGCGGGGAAATTGCCCTGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAg                                <  1:288232/131‑1 (MQ=255)                                                                                                                                     ccGCGGGGAAATTGCCCTGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAg                                >  2:288232/1‑131 (MQ=255)                                                                                                                                       gCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGc            <  6:153087/149‑1 (MQ=255)                                                                                                                                       gCAGGGCAATTTACCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGCCATCGTTGGCTCTCCGCGATTCAACCGATGGTGAACAGCGGATGTGCGACCTAAACCGGATGGCCGTCTTCTACCAGGATGGCGTCGCTGACAcctc            >  5:4743/1‑147 (MQ=255)                                                                                                                                         aGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGcc           >  1:317492/1‑148 (MQ=255)                                                                                                                                                  tCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGCCAATATCGGc  <  3:123598/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                AGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGCCAATATCGGC  >  NC_002947/649491‑649781

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