Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,715,401	A→T	10.9%	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,401	0	A	T	10.9%	161.6 / 9.0	54	D166E (GAT→GAA)	paaJ	3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base A (21/27);  new base T (4/2);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGGATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCAT  >  NC_002947/3715252‑3715550                                                                                                                                                       |                                                                                                                                                      gggATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaa                                                                                                                                                       >  1:609277/1‑150 (MQ=255)  ggATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaat                                                                                                                                                      >  2:774235/1‑150 (MQ=255)    aTTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaatcg                                                                                                                                                    <  2:1019378/150‑1 (MQ=255)    aTTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGAGCTGGCTGCGCAGGAGAAATTCA‑TCGGCGACTTCTCTTTAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCAGGCagcgaatcg                                                                                                                                                    >  1:274665/1‑150 (MQ=255)          aCCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCg                                                                                                                                              >  1:995571/1‑150 (MQ=255)          aCCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCCGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCg                                                                                                                                              <  2:824936/150‑1 (MQ=255)           ccGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGt                                                                                                                                             <  2:1048468/150‑1 (MQ=255)              gCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATg                                                                                                                                                 <  1:443641/143‑1 (MQ=255)              gCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATg                                                                                                                                                 >  2:443641/1‑143 (MQ=255)                aCGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAggcg                                                                                                                                        >  2:967547/1‑150 (MQ=255)                  gATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGc                                                                                                                                                <  2:921970/140‑1 (MQ=255)                  gATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGc                                                                                                                                                >  1:921970/1‑140 (MQ=255)                   ctttCCCGCGACAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTCCATGGGCA‑GACAGAGCAATTCCGGGTTGGCGCGCGACATTTCGTACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGcc                                                                                                                                     <  1:875443/149‑1 (MQ=255)                   aTTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGcc                                                                                                                                     <  1:904547/150‑1 (MQ=255)                        ccGCGCCAGGCGGCCGCGGGCCTGAGCTTCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTtcat                                                                                                                                >  1:96926/1‑150 (MQ=255)                            gcCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGcc                                                                                                                                     >  2:1019795/1‑141 (MQ=255)                            gcCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGcc                                                                                                                                     <  1:1019795/141‑1 (MQ=255)                                ggcggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGt                                                                                                                        <  1:967547/150‑1 (MQ=255)                                   ggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGa                                                                                                                     <  1:433111/150‑1 (MQ=255)                                           gCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGc                                                                                                             >  1:673665/1‑150 (MQ=255)                                             cTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCa                                                                                                           >  2:1030715/1‑150 (MQ=255)                                              tGAGCAGCTGAGGCTATGTGCTGCCTGGGCAGGTCAAAGTCA‑GGCGGGGCGGCGCGCGAGATGCCGAAATGTTCGGCCCCGTTTTCAGCCGTTCCCGTCATCGATTCTAATCCGTAGGCTGCCTTCATCATCGGGTTGACGAAACGCCAg                                                                                                          >  2:1167869/1‑150 (MQ=255)                                                      tgcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGaaa                                                                                                                >  1:1001526/1‑136 (MQ=255)                                                      tgcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGaaa                                                                                                                <  2:1001526/136‑1 (MQ=255)                                                      tgcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCa                                                                                                           <  1:652485/141‑1 (MQ=255)                                                      tgcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCa                                                                                                           >  2:652485/1‑141 (MQ=255)                                                       gcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtg                                                                                                 >  1:546097/1‑150 (MQ=255)                                                       gcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtg                                                                                                 <  1:923031/150‑1 (MQ=255)                                                       gcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtg                                                                                                 <  1:1167869/150‑1 (MQ=255)                                                        cGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGt                                                                                                                                             <  1:602960/105‑1 (MQ=255)                                                        cGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGt                                                                                                                                             >  2:602960/1‑105 (MQ=255)                                                            tttGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTc                                                                                            <  1:245723/150‑1 (MQ=255)                                                                 gctggctgCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACa                                                                                       <  2:995571/150‑1 (MQ=255)                                                                     gctgCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTc                                                                                   <  1:463020/150‑1 (MQ=255)                                                                       tgCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCgg                                                                                 <  1:995582/150‑1 (MQ=255)                                                                              ggCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtgg                                                                                                <  2:520921/128‑1 (MQ=255)                                                                              ggCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTGCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtgg                                                                                                >  1:520921/1‑128 (MQ=255)                                                                                       a‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAAAAGCTCGGCAACAAGACCGAATGc                                                                 >  2:68424/1‑150 (MQ=255)                                                                                            tGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGc                                                             <  1:480595/150‑1 (MQ=255)                                                                                               tCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTcc                                                          <  2:274665/150‑1 (MQ=255)                                                                                               tCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTcc                                                          <  1:228861/150‑1 (MQ=255)                                                                                                       ttcgaTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTAcc                                                  >  2:101631/4‑150 (MQ=255)                                                                                                         agaTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCTAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTAACCa                                                <  2:417135/150‑1 (MQ=255)                                                                                                          gaTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCg                                                                                                                                              >  1:615311/1‑55 (MQ=255)                                                                                                          gaTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCg                                                                                                                                              <  2:615311/55‑1 (MQ=255)                                                                                                           aTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATc                                              >  2:581892/1‑150 (MQ=255)                                                                                                             gCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtg                                                                                                 >  2:1057004/1‑97 (MQ=255)                                                                                                             gCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtg                                                                                                 <  1:1057004/97‑1 (MQ=255)                                                                                                                         ggCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTg                                                              <  2:225566/120‑1 (MQ=255)                                                                                                                         ggCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTg                                                              >  1:225566/1‑120 (MQ=255)                                                                                                                          gCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAcgc                               <  2:627071/150‑1 (MQ=255)                                                                                                                              cGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACa                           >  1:718791/1‑150 (MQ=255)                                                                                                                               gTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACAt                          <  2:609277/150‑1 (MQ=255)                                                                                                                                     aGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCAt                                               >  2:58203/1‑123 (MQ=255)                                                                                                                                     aGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCAt                                               <  1:58203/123‑1 (MQ=255)                                                                                                                                         gTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAAAGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCcaca                >  2:879154/1‑150 (MQ=255)                                                                                                                                           ctacGGCATCGATTCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCTTCTTACCCATCACAAACGGTGCACGCGACATGGACTCCAcacc              >  1:1036680/4‑150 (MQ=255)                                                                                                                                            tacGGCATCGATTCGATGCCGTAGGAGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATTGTGGTGTCGAACAGCGCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCAcacca             >  2:1042524/3‑150 (MQ=255)                                                                                                                                                      aatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGt                                     >  1:790886/1‑116 (MQ=255)                                                                                                                                                      aatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGt                                     <  2:790886/116‑1 (MQ=255)                                                                                                                                                       atcgaTGCCGTATGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCAt  <  2:604825/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      GGGATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCA‑TCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCAT  >  NC_002947/3715252‑3715550

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