Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,405,314	C→T	9.5%	G83R (GGG→AGG)	PP_3886 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,405,314	0	C	T	9.5%	166.3 / 6.8	64	G83R (GGG→AGG)	PP_3886	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base T (3/3);  total (34/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 = 3.39e-01

ATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCT  >  NC_002947/4405175‑4405460                                                                                                                                            |                                                                                                                                                     aTAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGCCTGGCTCCATTTAACGAAAAAGCCCCG‑GCAtgtg                                                                                                                                            >  3:55074/1‑149 (MQ=255) aTAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCAtgtg                                                                                                                                            <  6:326975/149‑1 (MQ=255)   aGCAAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCAGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGATATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGcc                                                                                                                                          >  8:54869/1‑149 (MQ=255)    gCCAGCGAACTGCCTGGGAACAGGGCGACCCCAGATCCCCCTATGACCCGTGACGGAATTTGGAGGGACTGCGCAGAAATCTGAGCATCTAACCTCTCGACGATATACATGACTGGCTCCATTTAACCAAAAAGCCCAC‑GCATGTGcc                                                                                                                                          >  5:191112/1‑148 (MQ=255)      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)                  tGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGcccc                                                                                                                                                     <  4:43075/124‑1 (MQ=255)                  tGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGcccc                                                                                                                                                     >  3:43075/1‑124 (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                                                                                                 <  3:424406/141‑1 (MQ=255)                                                   ccGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCg                                                                                                 >  4:424406/1‑141 (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)                                                              tttGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGACTCCATTTAACGAAAAAGCCCCg                                                                                                                                                    <  7:53037/81‑1 (MQ=255)                                                              tttGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGACTCCATTTAACGAAAAAGCCCCg                                                                                                                                                    >  8:53037/1‑81 (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                                                              >  5:381120/1‑149 (MQ=255)                                                                              aGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGAc                                                              >  7:205566/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                                                                                               >  2:93474/1‑109 (MQ=255)                                                                                     tGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCt                                                                                               <  1:93474/109‑1 (MQ=255)                                                                                          aTCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGACCGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTg                                                   <  1:130917/149‑1 (MQ=255)                                                                                           tCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCa                                                 >  3:127451/1‑149 (MQ=255)                                                                                            cTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCAc                                                >  1:118736/1‑149 (MQ=255)                                                                                            cTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCAc                                                >  4:244733/1‑149 (MQ=255)                                                                                              aaCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGt                                              >  7:398530/1‑149 (MQ=255)                                                                                              aaCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGt                                              >  3:106726/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                                >  7:184580/1‑120 (MQ=255)                                                                                                                                         gCCTCG‑GCACATGCCGGGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCttt                                <  8:184580/120‑1 (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)                                                                                                                                            |                                                                                                                                                     ATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCG‑GCATGTGCCGAGGCTTGGAATGGGTGGAGA‑TGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCT  >  NC_002947/4405175‑4405460

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