Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,405,327	A→G	9.5%	P78P (CCT→CCC)	PP_3886 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,405,327	0	A	G	9.5%	170.8 / 6.0	63	P78P (CCT→CCC)	PP_3886	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base G (3/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGC  >  NC_002947/4405180‑4405472                                                                                                                                                     |                                                                                                                                                   cAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAg                                                                                                                                                   <  7:319593/149‑1 (MQ=255)  aGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAgg                                                                                                                                                  >  6:315433/1‑149 (MQ=255)      aaCTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTg                                                                                                                                              <  7:153314/149‑1 (MQ=255)           ccTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCG‑GCACATGCCGGGGCtttttc                                                                                                                                           <  2:225868/147‑5 (MQ=255)           ccTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATg                                                                                                                                         >  3:354693/1‑149 (MQ=255)               ggAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATggg                                                                                                                                       <  5:186464/147‑1 (MQ=255)                 aaCAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGga                                                                                                                                   <  2:201435/149‑1 (MQ=255)                       gCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑Tgg                                                                                                                             >  8:305268/1‑148 (MQ=255)                         gACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCaa                                                                                                                          <  5:398443/149‑1 (MQ=255)                          aCCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCTTCTAACGTCTCGACTCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAAc                                                                                                                         <  1:352741/149‑1 (MQ=255)                           cccGAGTTCCGACGATGACCCGTGACGGAATTTGGAGGGACTGTGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGGCCCT‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAAcc                                                                                                                        <  6:191112/149‑1 (MQ=255)                            ccGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAAccc                                                                                                                       >  3:222020/1‑149 (MQ=255)                                gTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTc                                                                                                                   <  7:378094/149‑1 (MQ=255)                                gTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTc                                                                                                                   <  5:315433/149‑1 (MQ=255)                                 ttCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCa                                                                                                                  <  6:156757/149‑1 (MQ=255)                                    cgccgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCACG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGt                                                                                                               >  5:229963/1‑149 (MQ=255)                                       cgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGc                                                                                                            <  7:305268/149‑1 (MQ=255)                                             cccGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGt                                                                                                      <  1:172600/149‑1 (MQ=255)                                              ccGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCg                                                                                                             >  4:424406/1‑141 (MQ=255)                                              ccGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCg                                                                                                             <  3:424406/141‑1 (MQ=255)                                              ccGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTg                                                                                                     >  8:370340/1‑149 (MQ=255)                                               cGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTgg                                                                                                    >  2:130917/1‑149 (MQ=255)                                                         tttGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                                          <  4:222020/149‑1 (MQ=255)                                                         tttGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                                          >  6:153378/1‑149 (MQ=255)                                                            ggAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCG‑GCACATGCCGGGGCt                                                                                                                                                >  7:352348/1‑93 (MQ=255)                                                            ggAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCG‑GCACATGCCGGGGCt                                                                                                                                                <  8:352348/93‑1 (MQ=255)                                                            ggAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGTTTATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTcc                                                                                         <  1:135763/148‑1 (MQ=255)                                                            ggAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAg                                                                                       >  6:378621/1‑149 (MQ=255)                                                                 gACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTcccc                                                                                  <  7:370340/149‑1 (MQ=255)                                                                   cTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCAc                                                                                >  6:286533/1‑149 (MQ=255)                                                                         aGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGAc                                                                          >  7:205566/1‑149 (MQ=255)                                                                         aGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGAc                                                                          >  5:381120/1‑149 (MQ=255)                                                                         aGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGAc                                                                          >  5:275943/1‑149 (MQ=255)                                                                             aTCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGgcg                                                                                                  >  5:354183/1‑121 (MQ=255)                                                                             aTCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGgcg                                                                                                  <  6:354183/121‑1 (MQ=255)                                                                                tGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCt                                                                                                           <  1:93474/109‑1 (MQ=255)                                                                                tGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCt                                                                                                           >  2:93474/1‑109 (MQ=255)                                                                                     aTCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGACCGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTg                                                               <  1:130917/149‑1 (MQ=255)                                                                                      tCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCa                                                             >  3:127451/1‑149 (MQ=255)                                                                                       cTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCAc                                                            >  4:244733/1‑149 (MQ=255)                                                                                       cTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCAc                                                            >  1:118736/1‑149 (MQ=255)                                                                                         aaCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGt                                                          >  3:106726/1‑149 (MQ=255)                                                                                         aaCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGt                                                          >  7:398530/1‑149 (MQ=255)                                                                                                  aCGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCgg                                                 <  5:286533/149‑1 (MQ=255)                                                                                                         ccATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGAGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTCCCGCCGGGCtgtg                                           >  7:331319/1‑148 (MQ=255)                                                                                                           aTGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGACAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatc                                        <  7:70593/149‑1 (MQ=255)                                                                                                              aCTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatcatc                                     <  5:378621/149‑1 (MQ=255)                                                                                                                  gCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCAtctc                                   <  2:118736/147‑1 (MQ=255)                                                                                                                  gCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTcc                                  >  7:80095/1‑148 (MQ=255)                                                                                                                     ccATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCACTGCACGTTGCCGCCGGGCTTTGATCATCTCCaaa                               >  8:388704/1‑148 (MQ=255)                                                                                                                        tttAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAAc                            >  2:124767/1‑148 (MQ=255)                                                                                                                           aaCGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGAc                         <  5:153378/148‑1 (MQ=255)                                                                                                                             cGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCTCTTACCTCCAGTCCACACGCTGACTGTTACCCCTGCACGTTTCCTCCGGGCTTTGCTCATCTCCAAACACGACaaa                      >  8:145337/1‑149 (MQ=255)                                                                                                                              gAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg                     >  2:362824/1‑149 (MQ=255)                                                                                                                                aaaaGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGcc                   <  8:398530/149‑1 (MQ=255)                                                                                                                                 caaGCCTCG‑GCACATGCCGGGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGcc                   >  3:178459/2‑148 (MQ=255)                                                                                                                                    gCCTCG‑GCACATGCCGGGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCttt                                            <  8:184580/120‑1 (MQ=255)                                                                                                                                    gCCTCG‑GCACATGCCGGGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCttt                                            >  7:184580/1‑120 (MQ=255)                                                                                                                                     ccccG‑GCATGTGCCGAGGCTTGGAATGGGTTGAGA‑TGTCAACCCTGTCAGGTCGCTGTCGTGGCTCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACTTTGCCGCCGGGCTTTGATCATCTCCAAAAACGAAAAAGCCCAGCt              <  6:229963/149‑1 (MQ=255)                                                                                                                                     ccccG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCt              <  8:80095/149‑1 (MQ=255)                                                                                                                                     ccccG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGCTCGCTGTCGTGCCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCAAGCt              >  2:284472/1‑149 (MQ=255)                                                                                                                                                  ccGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGc  >  3:402244/1‑149 (MQ=255)                                                                                                                                                  ccGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGc  <  4:354693/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGC  >  NC_002947/4405180‑4405472

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