Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,798,704	A→T	6.0%	W66R (TGG→AGG)	PP_3358 ←	hydroxycinnamoyl‑CoA hydratase‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,798,704	0	A	T	6.0%	199.8 / 2.8	66	W66R (TGG→AGG)	PP_3358	hydroxycinnamoyl‑CoA hydratase‑lyase
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base T (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-01

GTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGAT  >  NC_002947/3798558‑3798845                                                                                                                                                    |                                                                                                                                              gTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAgg                                                                                                                                             <  5:165830/149‑1 (MQ=255) gTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAg                                                                                                                                              >  1:8925/1‑148 (MQ=255) gTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAg                                                                                                                                              >  7:219870/1‑148 (MQ=255) accgACCATCGCCATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAg                                                                                                                                              >  1:179892/4‑148 (MQ=255)  ttGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGa                                                                                                                                            <  8:165924/149‑1 (MQ=255)  ttGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGa                                                                                                                                            <  5:257410/149‑1 (MQ=255)  ttGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGa                                                                                                                                            <  4:341573/149‑1 (MQ=255)      ccATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCg                                                                                                                                        >  1:235459/1‑149 (MQ=255)       cATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGc                                                                                                                                       <  2:123630/149‑1 (MQ=255)       cATGGCGATGGTCGGGTTGGCAGACAGACGCACCACCGTCCCTTGCCCTTGCGAGGCGTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGGCCACCGCGCGGCAGTACTCCTTCAGGGCCATGCCGGCG‑GTCCAGGACTCGc                                                                                                                                       <  5:289148/149‑1 (MQ=255)             gATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAC‑CGCCGGCat                                                                                                                                 <  6:155548/149‑2 (MQ=255)                        tGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCTGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGCCTCGCCGGCACCGGTCAATag                                                                                                                       <  7:298833/148‑2 (MQ=255)                          gCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCa                                                                                                                     <  2:377817/148‑1 (MQ=255)                          gCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAg                                                                                                                    >  8:8315/1‑149 (MQ=255)                          gCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAg                                                                                                                    >  5:54193/1‑149 (MQ=255)                                aGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGcc                                                                                                              <  6:54193/149‑1 (MQ=255)                                aGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGc                                                                                                               >  2:297249/1‑148 (MQ=255)                                 gACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCca                                                                                                             <  2:235459/149‑1 (MQ=255)                                          aCTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtcagcg                                                                                                    >  6:114775/1‑149 (MQ=255)                                             tCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGc                                                                                                                   >  6:348042/1‑131 (MQ=255)                                             tCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGc                                                                                                                   <  5:348042/131‑1 (MQ=255)                                               cattgccattgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTg                                                                                               <  4:85839/149‑1 (MQ=255)                                                 ttgccattgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGc                                                                                              >  3:363004/1‑148 (MQ=255)                                                      attgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAgg                                                                                        <  8:336114/149‑1 (MQ=255)                                                      attgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAgg                                                                                        <  5:213692/149‑1 (MQ=255)                                                           gAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTc                                                                                   <  4:139918/149‑1 (MQ=255)                                                               cTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGc                                                                               <  2:179892/149‑1 (MQ=255)                                                                     gacgaATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCg                                                                         <  2:8925/149‑1 (MQ=255)                                                                        gaATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGAc                                                                       >  3:335/1‑148 (MQ=255)                                                                           tCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCAtt                                                                   >  8:204185/1‑149 (MQ=255)                                                                               ttCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCAttt                                                                  <  6:59821/146‑1 (MQ=255)                                                                                 ccTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                             >  2:91313/1‑149 (MQ=255)                                                                                     gaggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCa                                                         <  7:204185/149‑1 (MQ=255)                                                                                     gaggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCa                                                         <  1:297249/149‑1 (MQ=255)                                                                                        gaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAggg                                                      >  5:297226/1‑149 (MQ=255)                                                                                        gaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAgg                                                       <  4:335/148‑1 (MQ=255)                                                                                          tttCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGGTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTg                                                    >  3:40908/1‑149 (MQ=255)                                                                                          tttCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTcc                                                                <  4:158655/137‑1 (MQ=255)                                                                                          tttCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTcc                                                                >  3:158655/1‑137 (MQ=255)                                                                                          tttCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTg                                                    >  1:8434/1‑149 (MQ=255)                                                                                               gggCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCt                                               <  4:374271/149‑1 (MQ=255)                                                                                                   cGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATg                                           <  5:114775/149‑1 (MQ=255)                                                                                                    ggCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATgg                                          <  6:297226/149‑1 (MQ=255)                                                                                                         ccACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTg                                     >  7:253514/1‑149 (MQ=255)                                                                                                              tCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                                >  3:217094/1‑149 (MQ=255)                                                                                                                gcgGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGAAACCCCGGCACCGGTCAATACCAGCACGCCAGCGACAGCGACGTGGTAAAGGGTTTCAAGCACGTCGACCATTTCCCGGTACAGGTTGGGGCACATGGCATTGCGtttt                               >  2:344877/1‑148 (MQ=255)                                                                                                                      gTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGc                        >  3:29470/1‑149 (MQ=255)                                                                                                                          tCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACg                                                                            >  4:340707/1‑93 (MQ=255)                                                                                                                          tCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTAAAGCACGTCGACACTTACCCGGTTCATGGTGGGGCTCAAGGCATTGCGTTTTTCAGTGCGagt                    >  4:250648/1‑147 (MQ=255)                                                                                                                          tCCTTCAGGTCCATGCCGGCC‑GTCCACGACTCGCCGGCACCGGGCAATACCAGCACGCCAGCCTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGa                      <  6:193713/147‑1 (MQ=255)                                                                                                                            cTTCAGGTCCATGCCGGCT‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCTTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCTTTTTTCCgtt                         >  8:114327/1‑140 (MQ=255)                                                                                                                            cTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCggg                         <  7:114327/142‑1 (MQ=255)                                                                                                                            cTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGa                      >  5:193713/1‑145 (MQ=255)                                                                                                                            cTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGAtt                    >  3:281757/1‑147 (MQ=255)                                                                                                                               cAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCCAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAg                 <  4:29470/147‑1 (MQ=255)                                                                                                                               cAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAgg                <  2:8434/148‑1 (MQ=255)                                                                                                                               cAGCTCGATGGCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCATCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTACAGCACGTCGTCCATTTCCCGGTTCAGGTTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAg                 <  4:363004/147‑1 (MQ=255)                                                                                                                                aGGTCCATGGCGGCG‑GTCCAGGACTCGCCGGCACGGGTCAATACCAGCACGCCAGCGGCAGCGTCCTGCTCAAGGGTTTCCAGCACGGCGACCATTTCCGGGTTCCGTGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGt              <  8:219870/149‑1 (MQ=255)                                                                                                                                  gTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTggg                                                  <  3:83172/111‑1 (MQ=255)                                                                                                                                  gTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTggg                                                  >  4:83172/1‑111 (MQ=255)                                                                                                                                  gTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCa            >  3:7917/1‑149 (MQ=255)                                                                                                                                      gtGCCGGCGAGTCCTGGAC‑CGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCCCGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCAcc          >  2:179835/2‑147 (MQ=255)                                                                                                                                      gtGCCGGCGAGTCCTGGAC‑CGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAgg                                                       <  7:354193/101‑1 (MQ=38)                                                                                                                                      gtGCCGGCGAGTCCTGGAC‑CGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAgg                                                       >  8:354193/2‑102 (MQ=38)                                                                                                                                      aTGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCAcc          >  8:132309/1‑147 (MQ=255)                                                                                                                                            gcg‑gTCCAGGACTCGCCGGCACCTGTCAATACCAGCACGCCAGCGTCATCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGAt  <  4:40908/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              GTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGAT  >  NC_002947/3798558‑3798845

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