Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,586,526	A→G	7.1%	F254L (TTT→CTT)	PP_4914 ←	phytase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,586,526	0	A	G	7.1%	200.8 / 4.5	70	F254L (TTT→CTT)	PP_4914	phytase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base G (3/2);  total (40/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 = 9.65e-01

GATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTC  >  NC_002947/5586380‑5586670                                                                                                                                                   |                                                                                                                                                 gATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTaaa                                                                                                                                                  >  2:310175/1‑147 (MQ=255)  aTCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAg                                                                                                                                               >  1:270600/1‑149 (MQ=255)      ccAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTg                                                                                                                                           >  2:183818/1‑149 (MQ=255)         gTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCACCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGccc                                                                                                                                        >  7:254824/1‑149 (MQ=255)           gCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCt                                                                                                                                       <  2:270600/148‑1 (MQ=255)             ttCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCg                                                                                                                                     >  3:214694/1‑148 (MQ=255)              tCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGaa                                                                                                                                   >  7:8449/1‑149 (MQ=255)               cGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGa                                                                                                                                    <  6:316592/147‑1 (MQ=255)                  ccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCtt                                                                                                                                            >  6:290640/1‑136 (MQ=255)                  ccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCtt                                                                                                                                            <  5:290640/136‑1 (MQ=255)                  ccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccag                                                                                                                               >  4:298053/1‑149 (MQ=255)                      gccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggc                                                                                                                             >  4:116931/1‑147 (MQ=255)                       ccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggccag                                                                                                                          <  4:243502/149‑1 (MQ=255)                          aaCGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTc                                                                                                                       <  3:116931/149‑1 (MQ=255)                           aCGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGACATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCg                                                                                                                                     >  2:2331/1‑134 (MQ=255)                           aCGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGACATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCg                                                                                                                                     <  1:2331/134‑1 (MQ=255)                             ggCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGTGTAAACAGCTTGCCCTCGAACCAGGCCAGTTCAGc                                                                                                                    <  8:8449/149‑1 (MQ=255)                               cTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGGCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGa                                                                                                                  <  7:340644/149‑1 (MQ=255)                                       aGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTcgcgc                                                                                                          >  2:321852/1‑149 (MQ=255)                                          accGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGAc                                                                                                       <  8:254824/148‑1 (MQ=255)                                                aCTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCTGGCCAGGTCAGCGaa                                                                                                                 <  7:135638/133‑1 (MQ=255)                                                aCTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTAGAGGGTAAACAGCTTGCCCTCGAACCTGGCCAGGTCAGCGaa                                                                                                                 >  8:135638/1‑133 (MQ=255)                                                  tCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCt                                                                                               >  7:300622/1‑149 (MQ=255)                                                  tCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCt                                                                                               >  7:46855/1‑149 (MQ=255)                                                         gcgcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTg                                                                                        <  8:46855/149‑1 (MQ=255)                                                         gcgcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTg                                                                                        <  1:183818/149‑1 (MQ=255)                                                          cgcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTgg                                                                                       <  3:298053/149‑1 (MQ=255)                                                             aTCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCa                                                                                    >  1:322715/1‑149 (MQ=255)                                                                      gAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                           <  1:66007/149‑1 (MQ=255)                                                                            ccAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCt                                                                                 >  1:276968/1‑137 (MQ=255)                                                                            ccAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCt                                                                                 <  2:276968/137‑1 (MQ=255)                                                                             cAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAt                                                                    >  6:143198/1‑149 (MQ=255)                                                                               gCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAttt                                                                  <  1:310175/149‑1 (MQ=255)                                                                                   cGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGAc                                                              >  7:192777/1‑149 (MQ=255)                                                                                    gCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGAcc                                                             >  7:151579/1‑149 (MQ=255)                                                                                              gCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTc                                                     <  4:214694/147‑1 (MQ=255)                                                                                              gCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGc                                                   >  6:103703/1‑149 (MQ=255)                                                                                                cACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTGAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAttt                                                                  <  7:340383/132‑1 (MQ=255)                                                                                                cACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAttt                                                                  >  8:340383/1‑132 (MQ=255)                                                                                                cACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTc                                                 <  6:89981/149‑1 (MQ=255)                                                                                                     tCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagcag                                            >  2:188267/1‑149 (MQ=255)                                                                                                     tCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagc                                              <  5:143198/147‑1 (MQ=255)                                                                                                       gATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCTCCCCCCGCCTCGCTCAgcagca                                          >  6:195565/1‑149 (MQ=255)                                                                                                       gATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAgcagca                                          >  6:358616/1‑149 (MQ=255)                                                                                                       gATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAgcagca                                          >  6:363788/1‑149 (MQ=255)                                                                                                          atCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGAcc                                                             >  8:162658/1‑127 (MQ=255)                                                                                                          atCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGAcc                                                             <  7:162658/127‑1 (MQ=255)                                                                                                             aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCaa                                     <  1:321852/148‑1 (MQ=255)                                                                                                             aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAAc                                    >  8:121256/1‑149 (MQ=255)                                                                                                             aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAAc                                    <  2:322715/149‑1 (MQ=255)                                                                                                                   gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGcc                              >  4:14544/1‑149 (MQ=255)                                                                                                                    cacaCCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGccc                             <  3:49594/149‑1 (MQ=255)                                                                                                                     acacCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCg                            >  5:179422/1‑149 (MQ=255)                                                                                                                     acacCCGGTACGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCg                            >  5:85484/1‑149 (MQ=255)                                                                                                                         ccGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAgg                        >  1:269792/1‑149 (MQ=255)                                                                                                                            gTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTcc                     <  3:14544/149‑1 (MQ=255)                                                                                                                                  gTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCg               <  5:103703/149‑1 (MQ=255)                                                                                                                                   ttACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGa              >  1:215308/1‑149 (MQ=255)                                                                                                                                     tggTTCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGa                  <  7:243140/141‑1 (MQ=255)                                                                                                                                     tggTTCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGa                  >  8:243140/3‑143 (MQ=255)                                                                                                                                      cGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGActg           >  6:14315/1‑149 (MQ=255)                                                                                                                                       gTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGActgc          >  2:235809/1‑149 (MQ=255)                                                                                                                                        ttCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGActgct         >  8:210399/1‑149 (MQ=255)                                                                                                                                        ttCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCT‑GCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagca                                             >  4:357987/1‑112 (MQ=38)                                                                                                                                        ttCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCT‑GCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagca                                             <  3:357987/112‑1 (MQ=38)                                                                                                                                          cGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTgg       >  3:37614/1‑149 (MQ=255)                                                                                                                                           gAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTgg       <  2:215308/148‑1 (MQ=255)                                                                                                                                           gAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGgc      >  3:346172/1‑149 (MQ=255)                                                                                                                                               gTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTc  <  5:14315/149‑1 (MQ=255)                                                                                                                                               gTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACACGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCACCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTc  >  1:290025/1‑149 (MQ=255)                                                                                                                                                tAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTc  >  5:38598/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTC  >  NC_002947/5586380‑5586670

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