Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,531,153	A→C	7.2%	I117L (ATA→CTA)	PP_2219 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,531,153	0	A	C	7.2%	214.1 / 6.5	77	I117L (ATA→CTA)	PP_2219	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/37);  new base C (3/3);  total (37/40)
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

GCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGCGA  >  NC_002947/2531006‑2531296                                                                                                                                                     |                                                                                                                                                gctgagctgGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATa                                                                                                                                               <  1:2083418/150‑1 (MQ=255)  ctgagctgGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTATACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAg                                                                                                                                              <  1:2215187/150‑1 (MQ=255)   tgagctgGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAgg                                                                                                                                             >  2:363779/1‑150 (MQ=255)             tGCCGACCGGCGGATTGTGAGCAGGGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGc                                                                                                                                   <  1:931558/150‑1 (MQ=255)                 gACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACAt                                                                                                                               >  1:1834951/1‑150 (MQ=255)                    cggcggATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGaa                                                                                                                            >  1:873753/1‑150 (MQ=255)                      gcggATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATc                                                                                                                          <  2:1472755/150‑1 (MQ=255)                       cggATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCg                                                                                                                         <  2:1682540/150‑1 (MQ=255)                        ggATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCg                                                                                                                         >  2:1443509/1‑150 (MQ=255)                        ggATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCg                                                                                                                         >  2:1441551/1‑150 (MQ=255)                                aGCAGTTGAAGCAAATTCGCACCCATCTGCAGTTGCTCCACGGGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATCCCGCTGGCTATGCATAGGCATCGACGGCACATGAATCCCCCAAATGc                                                                                                                <  2:995399/150‑1 (MQ=255)                                aGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAATAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGc                                                                                                                <  1:2237773/150‑1 (MQ=255)                                aGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCCGTCGATGCCTATGCATCgccagc                                                                                                                                   >  1:919236/1‑126 (MQ=255)                                aGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCCGTCGATGCCTATGCATCgccagc                                                                                                                                   <  2:919236/131‑6 (MQ=255)                                 gCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGAGCAGCTTGTTCTTAGCAGAAAGCCGCCTGCGGTGCATAGGCCTCGACGGCACATGAATCGCCCCAATGcc                                                                                                               >  1:877193/1‑149 (MQ=255)                                        aaGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACAttttt                                                                                                        >  1:2012567/1‑150 (MQ=255)                                        aaGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACAttttt                                                                                                        >  2:165959/1‑150 (MQ=255)                                            aaaTTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCAt                                                                                                    >  2:3029/1‑150 (MQ=255)                                             aaTTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATa                                                                                                   <  2:1834951/150‑1 (MQ=255)                                             aaTTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATa                                                                                                   <  1:770453/150‑1 (MQ=255)                                               ttCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATaga                                                                                                 >  1:776048/1‑150 (MQ=255)                                                    aCCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAAt                                                                                            >  2:2418040/1‑150 (MQ=255)                                                    aCCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAAt                                                                                            >  1:1582592/1‑150 (MQ=255)                                                      ccATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCg                                                                                          >  2:486874/1‑150 (MQ=255)                                                      ccATTTGCAGTCGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCg                                                                                          >  2:292652/1‑150 (MQ=255)                                                        aTTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGAt                                                                                        <  1:1716570/150‑1 (MQ=255)                                                          ttGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGa                                                                                      <  2:776048/150‑1 (MQ=255)                                                                    tCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACa                                                                            >  2:858042/1‑150 (MQ=255)                                                                        cgcgGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTg                                                                        >  1:985523/1‑150 (MQ=255)                                                                        cgcgGCGTGGCGCCGTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCCCAGCAGTCGACGGCACCTGAATCGCCCAAATGCACCTTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTg                                                                        <  1:1993741/150‑1 (MQ=255)                                                                            gcgTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTgg                                                                    <  2:2012567/150‑1 (MQ=255)                                                                            gcgTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTgg                                                                    <  1:2418040/150‑1 (MQ=255)                                                                                ggCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCAt                                                                >  2:1016088/1‑150 (MQ=255)                                                                                 gcgcCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATg                                                               <  2:985523/150‑1 (MQ=255)                                                                                   gcCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGGCAATGAATCGCCCACATGCACATTTTTTCATAGAGACATCGATGACCAAGGGACAGTTGCTGGCCATGTc                                                             >  2:790245/1‑150 (MQ=255)                                                                                    ccTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACCAATTTTCATAGTGAGAGCTATTATCAAGTAACAGTTTCTGGCCATGTCg                                                            >  2:1592029/1‑150 (MQ=255)                                                                                     cTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGa                                                                                      <  1:1215132/123‑1 (MQ=255)                                                                                     cTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGa                                                                                      >  2:1732272/1‑123 (MQ=255)                                                                                     cTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGa                                                                                      <  1:1732272/123‑1 (MQ=255)                                                                                     cTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGa                                                                                      >  2:1215132/1‑123 (MQ=255)                                                                                        aGACTGCCGCCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGcc                                                         <  1:1441551/150‑1 (MQ=255)                                                                                        aGACTGCCGCCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGcc                                                         <  1:1443509/150‑1 (MQ=255)                                                                                         gACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGa                                                       >  1:1045056/1‑150 (MQ=255)                                                                                           cTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGc                                                                       >  1:1034962/1‑132 (MQ=255)                                                                                           cTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGc                                                                       <  2:1034962/132‑1 (MQ=255)                                                                                            ggccg‑ccCCCTGTAAAAGAGTGCCGGGTGGCGCTTGAACAATTCCGCTGGCGATGCATAGGCAGCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGa                                                    <  2:877193/149‑1 (MQ=255)                                                                                                g‑ccCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTAc                                                >  1:540495/1‑150 (MQ=255)                                                                                                    cccTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGaa                                                                                                                            >  1:1825164/1‑71 (MQ=255)                                                                                                    cccTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGaa                                                                                                                            <  2:1825164/71‑1 (MQ=255)                                                                                                    cccTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATg                                             >  1:577694/1‑150 (MQ=255)                                                                                                     ccTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGa                                            >  2:302098/1‑150 (MQ=255)                                                                                                       tGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAAc                                          <  1:858042/150‑1 (MQ=255)                                                                                                       tGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAAGCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAAc                                          <  1:1022785/150‑1 (MQ=255)                                                                                                         tAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGc                                        <  1:3029/150‑1 (MQ=255)                                                                                                            aaGAGTGCCTGTTTGCGCTTGAAGAATTCCGCTTGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGTAACACTTGCTGACCATATCGGCCCATGACTACATTAACGCCGa                                     <  2:873753/150‑1 (MQ=255)                                                                                                            aaGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAAt                                                                                                                           <  1:534048/64‑1 (MQ=255)                                                                                                            aaGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAAt                                                                                                                           >  2:534048/1‑64 (MQ=255)                                                                                                             agagTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCAGCTACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGAc                                    <  1:854694/150‑1 (MQ=255)                                                                                                                 tGCCGGTTGGCGCTTGAAGAATTTCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGc                                >  1:1400676/1‑150 (MQ=255)                                                                                                                 tGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGc                                <  2:540495/150‑1 (MQ=255)                                                                                                                   ccGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACCTTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCtt                              <  1:165959/150‑1 (MQ=255)                                                                                                                          gcgcTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCcttct                       >  1:1894647/1‑150 (MQ=255)                                                                                                                           cgcTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCcttctt                      <  1:1735404/150‑1 (MQ=255)                                                                                                                           cgcTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCcttctt                      <  2:1111615/150‑1 (MQ=255)                                                                                                                                gaagaaTTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCCAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTg                 >  2:802335/1‑150 (MQ=255)                                                                                                                                    aaTTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTTAgct             >  1:638407/1‑150 (MQ=255)                                                                                                                                       tCCGCTGGCGATGCATAGGCATCGACGGCATATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTg                                                                     >  2:1791236/1‑91 (MQ=255)                                                                                                                                       tCCGCTGGCGATGCATAGGCATCGACGGCATATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTg                                                                     <  1:1791236/91‑1 (MQ=255)                                                                                                                                        ccGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGctgctg         <  2:1136567/150‑1 (MQ=255)                                                                                                                                          gCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCa       <  2:244242/150‑1 (MQ=255)                                                                                                                                           ccgtcGATGCCTATGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAg      >  2:2118331/5‑150 (MQ=255)                                                                                                                                             gtcGATGCCTATGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGc    >  1:1668607/3‑150 (MQ=255)                                                                                                                                             gtcGATGCCTATGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGc    <  2:1668607/148‑1 (MQ=255)                                                                                                                                             ggCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGc    <  2:1894647/150‑1 (MQ=255)                                                                                                                                               cGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGCGa  <  1:486874/150‑1 (MQ=255)                                                                                                                                                gATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGt                                                              >  2:1286163/1‑89 (MQ=255)                                                                                                                                                gATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGt                                                              <  1:1286163/89‑1 (MQ=255)                                                                                                                                                     aTAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAgct             >  2:1634697/1‑133 (MQ=255)                                                                                                                                                     aTAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAgct             <  1:1634697/133‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                GCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGCGA  >  NC_002947/2531006‑2531296

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