Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,213,650	T→A	6.1%	A1105A (GCA→GCT)	PP_1061 ←	ATP‑dependent DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,213,650	0	T	A	6.1%	212.2 / 4.1	69	A1105A (GCA→GCT)	PP_1061	ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base A (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

ACACCAGCGGGGTACTGGTCAACGTGCTGCTCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCTTGCAGGCTGACCGGTTCGATCTCGCGGCGCA  >  NC_002947/1213502‑1213795                                                                                                                                                      |                                                                                                                                                  acacCAGCGGGGTACTGGTCAACGTGCTACTCGCCACNGTTGCCGCCNGCCGGCTCGAGGCAAACTGGCCGCTTCGGGAGGCTTCGGCGAGCCAGTGCGGACTGTAGTCTTTGATACGCCCCGGCGGCAGTTCCGCCTCCCAG‑GCCCCTg                                                                                                                                                  >  1:351906/1‑150 (MQ=255)    ccAGCGGGGTACTGGTCAACGTGCTGCTCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTgcgg                                                                                                                                               <  2:320616/150‑1 (MQ=255)          gggTACTGGTCAACGTGCTGCTCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTgg                                                                                                                                         <  1:995005/150‑1 (MQ=255)          gggTACTGGTCAACGTGCTGCTCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTgg                                                                                                                                         <  1:1181900/150‑1 (MQ=255)            gTACTGGTCAACGTGCTGCTCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAGCGCCGCAGCGGC‑CTggga                                                                                                                                       <  2:1023665/150‑2 (MQ=255)                tGGTCAACGTGCTGCTCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAGCGCCGCAGCGGC‑CTgggaggc                                                                                                                                    >  1:611152/1‑145 (MQ=255)                tGGTCAACGTGCTGCTCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAGCGCCGCAGCGGC‑CTgggaggc                                                                                                                                    <  2:611152/149‑5 (MQ=255)                         tgctgcTCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACt                                                                                                                          <  2:89572/150‑1 (MQ=255)                          gctgctCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTg                                                                                                                         <  1:335389/150‑1 (MQ=255)                            tgctCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGc                                                                                                                       >  1:362788/1‑150 (MQ=255)                             gctCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGcc                                                                                                                      >  2:680930/1‑150 (MQ=255)                                 gCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTgg                                                                                                                                         <  1:464816/127‑1 (MQ=255)                                 gCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTgg                                                                                                                                         >  2:464816/1‑127 (MQ=255)                                 gCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGc                                                                                                                  >  1:143443/1‑150 (MQ=255)                                 gCCACTCTTGCCGCCAGCCGGCTCCATGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGTACTGTAGTCTTGGAGACGCGCCGGCAGCAGTCCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACAGCCCGAgg                                                                                                                  <  2:351906/150‑2 (MQ=255)                                  ccACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCa                                                                                                                 >  2:916878/1‑150 (MQ=255)                                       gTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCg                                                                                                            <  2:120483/150‑1 (MQ=255)                                                  ccGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTg                                                                                                 >  2:809879/1‑150 (MQ=255)                                                      cTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGa                                                                                             <  2:1030425/150‑1 (MQ=255)                                                          aGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTgg                                                                                                                                         >  2:878508/1‑102 (MQ=255)                                                          aGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTgg                                                                                                                                         <  1:878508/102‑1 (MQ=255)                                                          aGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCAc                                                                                         >  1:917719/1‑150 (MQ=255)                                                          aGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCCGATAACTGGCCGAGCACTTCGGCCACCGCCTGCGgacgac                                                                                         >  1:918600/1‑150 (MQ=255)                                                                     ccGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCg                                                                              >  1:475482/1‑150 (MQ=255)                                                                        cTTCGGCAGGCTTCGTCGAGCCGGTGCGGAGGGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACGGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCa                                                                           <  2:78381/150‑1 (MQ=255)                                                                        cTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCa                                                                           <  1:395895/150‑1 (MQ=255)                                                                        cTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCa                                                                           <  2:362788/150‑1 (MQ=255)                                                                         ttCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGgcgc                                                                                 >  1:53357/1‑143 (MQ=255)                                                                         ttCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGgcgc                                                                                 <  2:53357/143‑1 (MQ=255)                                                                                ggCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAg                                                                   <  2:475482/150‑1 (MQ=255)                                                                                ggCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTACGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGAACTTCGGCCACCGCCTGCGGACCACGCAGGCGCCCGTCAGGCGCCAg                                                                   >  2:506268/1‑150 (MQ=255)                                                                                 gCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAgg                                                                  <  1:809879/150‑1 (MQ=255)                                                                                   ttCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTg                                                                                                                               >  1:790374/1‑87 (MQ=255)                                                                                   ttCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTg                                                                                                                               <  2:790374/87‑1 (MQ=255)                                                                                    tcttcgATCCAGTGCGGCCTGTCATGTTTCATTCGCGCTGTCTGGTGTCGCGCTTGCCAG‑TCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGtgc                                                               <  2:918600/147‑1 (MQ=255)                                                                                     cgtcgAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGt                                                                 >  1:1127146/1‑147 (MQ=255)                                                                                     cgtcgAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGt                                                                 <  2:1127146/147‑1 (MQ=255)                                                                                         gAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCa                                                          >  1:455147/1‑150 (MQ=255)                                                                                            ccAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATc                                                       <  2:455147/150‑1 (MQ=255)                                                                                               gTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTc                                                                            <  1:1105347/126‑1 (MQ=255)                                                                                               gTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTc                                                                            >  2:1105347/1‑126 (MQ=255)                                                                                               gTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGaa                                                    >  1:958330/1‑150 (MQ=255)                                                                                               gTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCATGCTCG‑GCCGGCGTGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGaa                                                    <  2:634450/150‑1 (MQ=255)                                                                                               gTGCGGACTGTAGTCTTTGATACGCGCCGACAGCAGTTCCGCCTACCAG‑GCCGCGGCGGCGCTGGCGAACCCTTGTAACTGGCCGACCGCTCCGGCCACCGCCTGCGGACCGCGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGaa                                                    <  2:345444/150‑1 (MQ=255)                                                                                                    gACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGaa                                               >  1:1116466/1‑150 (MQ=255)                                                                                                         tAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGcc                                                                                                                      <  2:1198540/74‑1 (MQ=255)                                                                                                         tAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGcc                                                                                                                      >  1:1198540/1‑74 (MQ=255)                                                                                                            tCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATg                                       >  1:779352/1‑150 (MQ=255)                                                                                                              tttGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGaa                                     >  2:509323/1‑150 (MQ=255)                                                                                                               ttGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAg                                    <  1:361399/150‑1 (MQ=255)                                                                                                                 gATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTc                                  <  2:779352/150‑1 (MQ=255)                                                                                                                   tACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCtt                                <  1:903474/150‑1 (MQ=255)                                                                                                                      gcgcCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGgcgc                                                                                 <  2:305625/98‑1 (MQ=255)                                                                                                                      gcgcCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGgcgc                                                                                 >  1:305625/1‑98 (MQ=255)                                                                                                                        gcCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCTTGCAgg                           <  1:240594/150‑1 (MQ=255)                                                                                                                          cGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCTTGCAGGCt                         >  1:1070910/1‑150 (MQ=255)                                                                                                                               gcagTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACCCAGGCGCTCGTCCGGCGCCAGGTGCTGCCACTCGAACAGACAGCGCCTGAAGTCTTGCAGGCTGACCg                    >  2:602349/1‑150 (MQ=255)                                                                                                                                   ttCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGgc                                                                        >  1:230987/1‑94 (MQ=255)                                                                                                                                   ttCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGgc                                                                        <  2:230987/94‑1 (MQ=255)                                                                                                                                     ccGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTc                                  <  1:511329/130‑1 (MQ=255)                                                                                                                                     ccGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTc                                  >  2:511329/1‑130 (MQ=255)                                                                                                                                      cGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTTGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCTTGCAGGCTGACCGGtt                 <  2:7659/146‑1 (MQ=255)                                                                                                                                      cGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTTGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCTTGCAGGCTGACCGGtt                 >  1:7659/1‑146 (MQ=255)                                                                                                                                           cccagcgCCGCAGCGGC‑CTGGGGAACCCTTGTAACTGGCCGGGCACTTGGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCTTGCAGGCTGACCGGTTCGATCTcgc        >  1:881297/7‑150 (MQ=255)                                                                                                                                            ccAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCTTGCAGGCTGACCGGTTCGATCTcgcg       >  1:701853/1‑150 (MQ=255)                                                                                                                                             cAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCTTGCAGGCTGACCGGTTCGATCTCgcgg      >  1:985152/1‑150 (MQ=255)                                                                                                                                                  ccGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGaa                                     <  2:550857/115‑1 (MQ=255)                                                                                                                                                  ccGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGaa                                     >  1:550857/1‑115 (MQ=255)                                                                                                                                                  ccGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCAACGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCTTGCAGGCTGACCGGTTCGATCTCGCGGCGCa  >  1:336509/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ACACCAGCGGGGTACTGGTCAACGTGCTGCTCGCCACTGTTGCCGCCAGCCGGCTCCAGGCAAACTGGCCGCTTCGGCAGGCTTCGTCGAGCCAGTGCGGACTGTAGTCTTTGATACGCGCCGGCAGCAGTTCCGCCTCCCAG‑GCCGCTGCGGCGCTGGGGAACCCTTGTAACTGGCCGAGCACTTCGGCCACCGCCTGCGGACCACGCAGGCGCTCGTCAGGCGCCAGGTGCTGCCAATCGAACAGAAAGCGCATGAAGTCTTGCAGGCTGACCGGTTCGATCTCGCGGCGCA  >  NC_002947/1213502‑1213795

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑