Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,468,261	G→C	8.3%	intergenic (‑25/+98)	PP_2161 ← / ← PP_2162	hypothetical protein/phosphonoacetaldehyde phosphonohydrolase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,468,261	0	G	C	8.3%	169.1 / 4.8	59	intergenic (‑25/+98)	PP_2161/PP_2162	hypothetical protein/phosphonoacetaldehyde phosphonohydrolase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base C (2/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGCCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACCAGGCAG  >  NC_002947/2468117‑2468407                                                                                                                                                 |                                                                                                                                                    cgcCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCaa                                                                                                                                                  <  5:305802/148‑1 (MQ=255) cgcCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAg                                                                                                                                                 >  6:287600/1‑149 (MQ=255)   ccTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGc                                                                                                                                                >  1:267714/1‑148 (MQ=255)     gcGATCTGGCGCGCCTCTTCGCGCAGGGCGGCCAGCTCTTCTTCCTTAAGCGTTACGTTGCGGTCCAGTTGCTGCCGTAATCTGTTCAGGCGCTCTTGCAATTCGCGGGCTGGCTTTGGATATTCTCCTTCATTGGCTTGCCaaccaat                                                                                                                                             <  4:196793/148‑6 (MQ=255)             gcgcgcCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTTCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGCCAAGCCAttg                                                                                                                                           <  5:132444/143‑3 (MQ=255)             gcgcgcCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGa                                                                                                                                          <  3:2503/144‑1 (MQ=255)             gcgcgcCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACctct                                                                                                                                     >  2:183156/1‑149 (MQ=255)             gcgcgcCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGCCAAGCCATTGa                                                                                                                                          <  2:372793/144‑1 (MQ=255)                        cGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTgg                                                                                                                            <  2:49956/147‑1 (MQ=255)                          ggCAGGGCGTCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCtaa                                                                                                                        <  2:375446/148‑1 (MQ=255)                            cAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAgg                                                                                                                      >  6:96807/1‑149 (MQ=255)                             aGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGt                                                                                                                     >  4:19838/1‑149 (MQ=255)                                   gcCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGc                                                                                                               >  1:147908/1‑149 (MQ=255)                                    ccAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTAAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGcc                                                                                                              >  5:104815/1‑149 (MQ=255)                                           ttcttcCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGATATTCTCCATCAATGGCTTGGCAAGCAATGGACCTCTGCGAAGTGGCAAAGGTTCTGGCCTGTCtt                                                                                                        >  1:186106/1‑148 (MQ=255)                                             catcCGAAAGCGCAACGTTGCGTTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAAATCCACGGCAGGCATTGGCTAATCTCCATCAATCGCTACTCAAGCCATGGACATCTGAGACGTGGCAAAGGTTCTGGCCTGTATTCAg                                                                                                     >  6:128557/3‑149 (MQ=255)                                              ttcCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAgg                                                                                                    >  5:88387/1‑149 (MQ=255)                                                ccGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAgg                                                                                                    <  2:267714/147‑1 (MQ=255)                                                ccGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGAAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAggg                                                                                                   >  2:351388/1‑148 (MQ=255)                                                 cGAAAGCGTTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTa                                                                                                 <  5:215875/149‑1 (MQ=255)                                                 cGAAAGCGGTACGGTGGGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTa                                                                                                 <  3:196909/149‑1 (MQ=255)                                                  gAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTaa                                                                                                <  3:12024/149‑1 (MQ=255)                                                  gAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTaa                                                                                                >  6:7728/1‑149 (MQ=255)                                                       cGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTAc                                                                                          <  5:96807/149‑1 (MQ=255)                                                        ggTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTAcc                                                                                         >  4:213330/1‑149 (MQ=255)                                                             gTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACctgctg                                                                                    <  5:287600/149‑1 (MQ=255)                                                                  ggTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGAt                                                                               >  5:115208/1‑149 (MQ=255)                                                                      cAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTg                                                                            <  1:164508/148‑1 (MQ=255)                                                                        gTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCAt                                                                         >  1:208676/1‑149 (MQ=255)                                                                        gTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACATCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGATGGCGATTTGCAt                                                                         >  1:176905/1‑149 (MQ=255)                                                                         ttGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATg                                                                        <  1:351388/149‑1 (MQ=255)                                                                           gCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGa                                                                       <  3:19838/148‑1 (MQ=255)                                                                           gCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGa                                                                       <  2:147908/148‑1 (MQ=255)                                                                               gcgcAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCcgc                                                                  >  4:187465/1‑149 (MQ=255)                                                                                 gcAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGG‑AA‑TCTACCTGCTGGCGATTCGCATGATCCGCGTc                                                               <  3:115180/149‑1 (MQ=255)                                                                                    aGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCACTGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCaa                                                             <  5:128557/149‑1 (MQ=255)                                                                                    aGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCaa                                                             <  3:213330/149‑1 (MQ=255)                                                                                        gTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACg                                                         >  3:163867/1‑149 (MQ=255)                                                                                                 ctctTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCa                                                <  2:96212/149‑1 (MQ=255)                                                                                                      gCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAgg                                              <  6:88387/146‑1 (MQ=255)                                                                                                             gcgGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCAcc                                    >  4:233922/1‑149 (MQ=255)                                                                                                             gcgGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCAc                                     >  4:205571/1‑148 (MQ=255)                                                                                                               gggCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGc                                                       >  1:334125/1‑128 (MQ=255)                                                                                                               gggCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGc                                                       <  2:334125/128‑1 (MQ=255)                                                                                                                      caTTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCg                           <  5:7728/149‑1 (MQ=255)                                                                                                                            cTATTCTCCATCAATGGCTTGGCAAGCAATGGACCTCTGCGACGTGGCAAAGGTTCTGGCGTGTCTTCAGGGTAA‑TCTACCTGATGGCGATTTGCATGATCAGCGTCAAGACGGCTGTATCAGGGTGTGTCAACGTTGCGCCGGCGCTg                     >  6:99689/1‑149 (MQ=255)                                                                                                                             tATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGtt                                           <  6:276482/126‑1 (MQ=255)                                                                                                                             tATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTgg                    >  2:76324/1‑149 (MQ=255)                                                                                                                               ttCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTgg                    <  6:104815/147‑1 (MQ=255)                                                                                                                                  tcCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGAt               >  7:118427/1‑149 (MQ=255)                                                                                                                                    cATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGAt               <  2:208676/147‑1 (MQ=255)                                                                                                                                     gtccaTGGCTTGCCAAGCCATTGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTCGCGAGTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTCGGCGATGTc            >  5:280059/5‑149 (MQ=255)                                                                                                                                        caTGGCTTGCCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTccac         >  5:132691/2‑149 (MQ=255)                                                                                                                                        aaTGGCTTGTCAAGGCATGGACCTATGCGACGTGGCAAAGGTTCTTGCCTTTCTTCATGGAAA‑TCTACATGCTGGCGATTTTCATTATCCGCGTCAAGACTGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCTATGTccac         <  2:186106/149‑1 (MQ=255)                                                                                                                                            gCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGAt               <  1:138491/139‑1 (MQ=255)                                                                                                                                            gCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGAt               >  2:138491/1‑139 (MQ=255)                                                                                                                                            gCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcag      <  3:233922/148‑1 (MQ=255)                                                                                                                                             cTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcagg     >  7:226889/1‑148 (MQ=255)                                                                                                                                              ttGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAAGTCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcagg     >  1:106252/1‑148 (MQ=255)                                                                                                                                                ggCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcaggcag  <  4:163867/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CGCCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAA‑TCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACCAGGCAG  >  NC_002947/2468117‑2468407

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