Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,715,358	A→T	11.7%	D166E (GAT→GAA)	paaJ ←	3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,715,358	0	A	T	11.7%	150.8 / 11.2	60	D166E (GAT→GAA)	paaJ	3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base T (3/4);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

TTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACC  >  NC_002947/3715213‑3715498                                                                                                                                                  |                                                                                                                                              tttCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgattc‑g                                                                                                                                           <  5:279482/149‑1 (MQ=255)     gACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTTTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGATTTCAGCCGTTTCCGGCATCGAATC‑GATGc                                                                                                                                       <  7:335196/149‑1 (MQ=255)       ccGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTC‑GATGCCg                                                                                                                                     <  5:106043/149‑1 (MQ=255)       ccGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTC‑GATGCCg                                                                                                                                     <  2:352551/149‑1 (MQ=255)         ggCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTC‑GATGCCg                                                                                                                                     <  7:165698/147‑1 (MQ=255)         ggCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTC‑GATGCCg                                                                                                                                     >  8:165698/1‑147 (MQ=255)          gCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaat                                                                                                                                              >  6:323324/1‑138 (MQ=255)          gCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaat                                                                                                                                              <  5:323324/138‑1 (MQ=255)           cACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGTTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAgg                                                                                                                                 >  7:55935/1‑149 (MQ=255)            aCGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAggc                                                                                                                                >  6:67312/1‑149 (MQ=255)              gATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAggcgg                                                                                                                              >  2:214009/1‑149 (MQ=255)                   cccGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCATGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCtt                                                                                                                          >  7:298147/1‑148 (MQ=255)                    ccGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTC‑GATGCCg                                                                                                                                     >  2:351845/1‑136 (MQ=255)                    ccGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTC‑GATGCCg                                                                                                                                     <  1:351845/136‑1 (MQ=255)                    ccGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTtca                                                                                                                        >  3:297489/1‑149 (MQ=255)                            ggcggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCggg                                                                                                                <  5:67312/149‑1 (MQ=255)                            ggcggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCggg                                                                                                                <  6:102083/149‑1 (MQ=255)                              cggcAGCGGGACTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCTGGtt                                                                                                              <  4:212498/149‑1 (MQ=255)                               ggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTg                                                                                                             <  5:223077/149‑1 (MQ=255)                                         cTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCt                                                                                                                           >  5:143544/1‑125 (MQ=255)                                         cTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCt                                                                                                                           <  6:143544/125‑1 (MQ=255)                                         cTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGc                                                                                                    >  1:90748/1‑148 (MQ=255)                                            aGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGc                                                                                                <  1:214009/149‑1 (MQ=255)                                             gctgctGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGcc                                                                                               <  2:90748/149‑1 (MQ=255)                                                 ctgcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGa                                                                                             <  5:288170/147‑1 (MQ=255)                                                 ctgcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtg                                                                                           >  1:87818/1‑149 (MQ=255)                                                 ctgcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtg                                                                                           <  4:22380/149‑1 (MQ=255)                                                     ggCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCTGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtg                                                                                        <  1:282291/149‑1 (MQ=255)                                                      gCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATggtggt                                                                                      <  2:87818/149‑1 (MQ=255)                                                             gctggctgCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGaa                                                                                >  7:232682/1‑148 (MQ=255)                                                                     cgcAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGc                                                                       >  6:370940/1‑149 (MQ=255)                                                                        aGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGcctc                                                                    <  2:50223/149‑3 (MQ=255)                                                                               aGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGGAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGccc                                                              >  1:90098/1‑148 (MQ=255)                                                                               aGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGGAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGccc                                                              <  2:90098/148‑1 (MQ=255)                                                                                     ccTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGcc                                                       >  3:172197/1‑149 (MQ=255)                                                                                      cTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGcc                                                       <  4:297489/148‑1 (MQ=255)                                                                                         gTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCt                                                   >  5:192005/1‑149 (MQ=255)                                                                                             gcgcgcGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGc                                                        <  4:112535/140‑1 (MQ=255)                                                                                             gcgcgcGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGc                                                        >  3:112535/1‑140 (MQ=255)                                                                                              cgcgcgAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGAc                                              <  4:172197/149‑1 (MQ=255)                                                                                               gcgcgAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAGTGCCTGCTCCGACt                                             >  2:122373/1‑149 (MQ=255)                                                                                               gcgcgAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAAt                                                          <  4:264066/136‑1 (MQ=255)                                                                                               gcgcgAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAAt                                                          >  3:264066/1‑136 (MQ=255)                                                                                                  cgAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTAc                                          <  8:232682/149‑1 (MQ=255)                                                                                                             aCTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAAc                               >  6:304891/1‑149 (MQ=255)                                                                                                              cTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACg                              <  5:370940/149‑1 (MQ=255)                                                                                                                 ttCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTg                           >  6:136939/1‑149 (MQ=255)                                                                                                                      ccACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAcgc                      >  7:262486/1‑149 (MQ=255)                                                                                                                      ccACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAcgc                      >  5:14681/1‑149 (MQ=255)                                                                                                                       cACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAAc                                                                               >  5:23286/1‑91 (MQ=255)                                                                                                                       cACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAAc                                                                               <  6:23286/91‑1 (MQ=255)                                                                                                                       cACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAcgcg                     >  6:297245/1‑149 (MQ=255)                                                                                                                        aCGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAAc                                                                               >  4:71597/1‑90 (MQ=255)                                                                                                                        aCGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAAc                                                                               <  3:71597/90‑1 (MQ=255)                                                                                                                           ttttCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAGTGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACAt                 <  1:122373/149‑1 (MQ=255)                                                                                                                             ttCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATg                <  8:55935/148‑1 (MQ=255)                                                                                                                                aGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGAc                                              >  2:247622/1‑115 (MQ=255)                                                                                                                                aGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGAc                                              <  1:247622/115‑1 (MQ=255)                                                                                                                                     cctacGGCATCGATTC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCcaca       >  2:11063/5‑149 (MQ=255)                                                                                                                                         cGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTc                                                  >  2:266124/1‑102 (MQ=255)                                                                                                                                         cGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTc                                                  <  1:266124/102‑1 (MQ=255)                                                                                                                                         cGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCAcaccac   >  8:33478/1‑149 (MQ=255)                                                                                                                                          ggCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACaccacc  >  4:154771/1‑149 (MQ=255)                                                                                                                                                 aatc‑gaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCAcacca    <  6:220926/140‑1 (MQ=255)                                                                                                                                                 aatc‑gaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCAcacca    >  5:220926/1‑140 (MQ=255)                                                                                                                                                  |                                                                                                                                              TTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATC‑GATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACC  >  NC_002947/3715213‑3715498

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