Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	649,625	A→G	7.2%	L18P (CTG→CCG)	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,625	0	A	G	7.2%	199.7 / 3.6	69	L18P (CTG→CCG)	accC	acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (35/28);  new base G (3/2);  total (38/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.91e-01

GATATGCAGGTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGCC  >  NC_002947/649480‑649772                                                                                                                                                    |                                                                                                                                                    gATATGCAGGTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACGGACTCGTGTGCCAGGCCCAGGTGCATCAGTTGACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAggg                                                                                                                                                  <  3:370255/149‑1 (MQ=255)        aGGTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGACTTGATACACAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTc                                                                                                                                           >  5:283274/1‑149 (MQ=255)         ggTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTcc                                                                                                                                          >  7:180868/1‑149 (MQ=255)          gTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTccc                                                                                                                                         >  3:264126/1‑149 (MQ=255)          gTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTccc                                                                                                                                         >  4:137050/1‑149 (MQ=255)          gTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTccc                                                                                                                                         >  5:355332/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:268422/146‑1 (MQ=255)                                                         ccAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGtt                                                                                             <  8:116139/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                                                                           <  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                                     >  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)                                                                                                                                                    |                                                                                                                                                    GATATGCAGGTAGGAATCCTTGGACGATGCAGGGCCGATGCAGACCGACTCGTCTGCCAGGCCCAGGTGCATCAGTTCACGGTCGGCAGTGGAGTGCACG‑‑GCGACGGTCTTGATACCCAGCTCTTTGCAGGCACGCAGGATCCGCAGGGCAATTTCCCCGCGGTTGGCGATCAGGACTTTTTCGAGCTTCCCAGACATCGTTGGCTCTCCGCGATTCAAACGATGGTGAACAGCGGCTGGTCGAACTCAACCGGCTGGCCGTCTTCTACCAGGATGGCGTCGATGACACCGCC  >  NC_002947/649480‑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: ‑