Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,586,387	A→T	6.2%	L300Q (CTG→CAG)	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,387	0	A	T	6.2%	188.2 / 2.7	64	L300Q (CTG→CAG)	PP_4914	phytase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base T (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

CCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAA  >  NC_002947/5586247‑5586526                                                                                                                                             |                                                                                                                                             cccGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACTGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTc                                                                                                                                       <  7:266432/148‑1 (MQ=255)       gcCTGGCTCATGCTTTGGCGCTCCAGCCGCCTTCCGGTGAGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACCGGCCCCCTTGTAGTTGTCCAGGCCAATCCAGGCGCCCATGGCATCGATAACCAGTGCTTCggccag‑g                                                                                                                               >  6:233915/1‑149 (MQ=255)          tGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCggccag‑gcc                                                                                                                             >  3:214128/1‑148 (MQ=255)               cATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACgg                                                                                                                       >  8:277867/1‑149 (MQ=255)                aTGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCggc                                                                                                                                    <  4:220921/136‑1 (MQ=255)                aTGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCggc                                                                                                                                    >  3:220921/1‑136 (MQ=255)                          gcTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAgc                                                                                                            <  1:83802/149‑1 (MQ=255)                            tCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGcgcg                                                                                                                <  5:171783/143‑1 (MQ=255)                            tCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGcgcg                                                                                                                >  6:171783/1‑143 (MQ=255)                              cAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAAc                                                                                                                         >  7:26807/1‑132 (MQ=255)                              cAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAAc                                                                                                                         <  8:26807/132‑1 (MQ=255)                              cAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCg                                                                                                        >  6:306536/1‑149 (MQ=255)                               aGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGt                                                                                                       >  6:175071/1‑149 (MQ=255)                                  cgccTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCggccag                                                                                                                                 >  7:300734/1‑121 (MQ=255)                                  cgccTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCggccag                                                                                                                                 <  8:300734/121‑1 (MQ=255)                                  cgccTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGAc                                                                                                    >  2:89627/1‑149 (MQ=255)                                     cTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTc                                                                                                  >  7:80973/1‑148 (MQ=255)                                     cTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCg                                                                                                 >  5:200310/1‑149 (MQ=255)                                           gTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAgc                                                                                           >  8:16885/1‑149 (MQ=255)                                           gTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAgc                                                                                           >  1:211776/1‑149 (MQ=255)                                                 cGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCAt                                                                                       >  2:100480/1‑147 (MQ=255)                                                  gAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATcag                                                                                    >  4:220162/1‑149 (MQ=255)                                                    aGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATcagca                                                                                  <  6:129279/149‑1 (MQ=255)                                                      cgcCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCagcagc                                                                                <  4:210187/149‑1 (MQ=255)                                                         cAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGAGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGaa                                                                             >  4:276286/1‑149 (MQ=255)                                                          aGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAGCGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGaa                                                                             <  1:89627/149‑1 (MQ=255)                                                          aGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAAc                                                                            <  4:214128/149‑1 (MQ=255)                                                              gacCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGAcc                                                                        <  3:6445/149‑1 (MQ=255)                                                               agcggccggGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCCGCCAG‑GCCTAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCa                                                                       <  5:233915/147‑1 (MQ=255)                                                               acCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCa                                                                       <  7:16885/149‑1 (MQ=255)                                                                   gccggGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCaa                                                                   >  2:223993/1‑149 (MQ=255)                                                                    ccggGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAAc                                                                  <  3:220162/149‑1 (MQ=255)                                                                    ccggGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAAc                                                                  <  5:306536/149‑1 (MQ=255)                                                                       ggTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCggc                                                                                                                                    <  4:13749/81‑1 (MQ=255)                                                                       ggTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCggc                                                                                                                                    >  3:13749/1‑81 (MQ=255)                                                                       ggTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCt                                                               <  1:100480/149‑1 (MQ=255)                                                                            cACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGAc                                                           <  6:200310/148‑1 (MQ=255)                                                                                    gCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACt                                                  >  2:257754/1‑149 (MQ=255)                                                                                      aCGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGt                                                >  2:199180/1‑149 (MQ=255)                                                                                        gggCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTc                                               >  4:115090/1‑148 (MQ=255)                                                                                                ttgttgttgTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGa                                                            >  4:37570/1‑127 (MQ=255)                                                                                                ttgttgttgTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGa                                                            <  3:37570/127‑1 (MQ=255)                                                                                                  gttgttgtCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGt                                     <  1:12550/148‑1 (MQ=255)                                                                                                   ttgttgtCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGCCCTTGCCACTGTCGATATCACGTCg                                   >  4:213941/1‑149 (MQ=255)                                                                                                           cAGGCCAATCCAGGCGCCCTTGGCATCGATCACCTGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGCCCTTGCCACTGTCGATATCCCGTCGGCACAccc                           >  5:298467/1‑149 (MQ=255)                                                                                                           cAGGCCAATCCAGGCGCCCTTGGCATCGATCACCTGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACAccc                           >  5:303664/1‑149 (MQ=255)                                                                                                           cAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACAccc                           >  1:44014/1‑149 (MQ=255)                                                                                                           cAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACAcc                            <  1:176451/148‑1 (MQ=255)                                                                                                              gCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGCACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCgg                         >  3:283257/1‑148 (MQ=255)                                                                                                              gCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCgg                         >  7:80142/1‑148 (MQ=255)                                                                                                                cAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAg                      >  2:30567/1‑149 (MQ=255)                                                                                                                cAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAg                      <  5:175071/149‑1 (MQ=255)                                                                                                                  acaacggCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGc                    <  7:251/144‑1 (MQ=255)                                                                                                                      aGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTa                <  3:115090/149‑1 (MQ=255)                                                                                                                             cTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGTCTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGa         <  8:80973/149‑1 (MQ=255)                                                                                                                             cTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGa         <  7:277867/149‑1 (MQ=255)                                                                                                                              ttGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAg        <  8:80142/149‑1 (MQ=255)                                                                                                                                ggCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAgg       <  5:144430/148‑1 (MQ=255)                                                                                                                                  cATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTa    <  1:199180/149‑1 (MQ=255)                                                                                                                                   aTCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTa    >  3:156878/1‑148 (MQ=255)                                                                                                                                    tCGATCACCTGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTaaa  <  6:303664/149‑1 (MQ=255)                                                                                                                                    tCGATCACCTGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTaaa  <  6:298467/149‑1 (MQ=255)                                                                                                                                       aTCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCa                                        <  6:254924/108‑1 (MQ=255)                                                                                                                                       aTCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCa                                        >  5:254924/1‑108 (MQ=255)                                                                                                                                             |                                                                                                                                             CCCGGCGCCTGGCTCATGGTTTGGCGCTCCAGCCGCCTTCCGGTGCGTCGAAGCGCCAGACGACCGGCCGGGTTTCACCGTCGGCACGGGCGCCGTTGTTGTTGTCCAGGCCAATCCAGGCGCCCTTGGCATCGATCACCAGTGCTTCGGCCAG‑GCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAA  >  NC_002947/5586247‑5586526

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