Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,383,254	A→T	18.2%	L16Q (CTG→CAG)	PP_1204 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,383,254	0	A	T	18.2%	102.5 / 17.7	59	L16Q (CTG→CAG)	PP_1204	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (19/27);  new base T (4/7);  total (24/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 = 4.13e-01

CACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTAC  >  NC_002947/1383106‑1383399                                                                                                                                                      |                                                                                                                                                  cacaCGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGcaca                                                                                                                                                   <  5:133645/149‑1 (MQ=255)  acacGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAg                                                                                                                                                  <  4:209162/149‑1 (MQ=255)  acacGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAg                                                                                                                                                  <  4:295599/149‑1 (MQ=255)           tGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGcc                                                                                                                                         >  1:322899/1‑149 (MQ=255)             cGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGcc                                                                                                                                         <  7:84128/147‑1 (MQ=255)             cGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGcc                                                                                                                                         <  1:289313/147‑1 (MQ=255)             cGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGcc                                                                                                                                         <  2:365240/147‑1 (MQ=255)             cGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGcccc                                                                                                                                       <  4:271570/149‑1 (MQ=255)                   cTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑ATCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAggg                                                                                                                                 <  1:70296/149‑1 (MQ=255)                    tCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGCCGAGCAGGGCGCCGCACCGTGCGGCGCCCCCAAggg                                                                                                                                 >  3:93980/1‑148 (MQ=255)                     ccATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACTGTGCGgcgc                                                                                                                                          >  3:344703/1‑138 (MQ=255)                     ccATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACTGTGCGgcgc                                                                                                                                          <  4:344703/138‑1 (MQ=255)                      cATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTa                                                                                                                               >  7:23800/1‑148 (MQ=255)                      cATCACCTTCGGCGCGCCGACATCGATCTCCAGGGTTTGCTCGAGCAGTTTGCGCTGGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTGGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTaa                                                                                                                              <  6:211233/149‑1 (MQ=255)                        tCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAg                                                                                                                            <  5:299758/149‑1 (MQ=255)                        tCACCTTCGGCCCGCCGACAGCGATCTGCAGGGTATGCGCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAg                                                                                                                            <  8:98460/149‑1 (MQ=255)                          aCCTTCGGGTGGCCGACATCGTTGTGCTGGGTATGCTCGAGCAGTTTGCGCTCGACGGGGTCCTCATGGTCGCTGCGGTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGt                                                                                                                          <  1:328566/149‑1 (MQ=255)                            cTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTAc                                                                                                                        <  2:322899/149‑1 (MQ=255)                             ttCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACt                                                                                                                       >  1:203391/1‑149 (MQ=255)                                        gACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACa                                                                                                            >  8:2765/1‑149 (MQ=255)                                         aCATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACAt                                                                                                           <  2:203391/149‑1 (MQ=255)                                         aCATCGATCTGCAGGGTATGCGCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACAt                                                                                                           <  4:55554/149‑1 (MQ=255)                                                   gCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGt                                                                                                 >  5:187809/1‑149 (MQ=255)                                                    cagtttttGCTCGAGCAGTTTTGGCGCGACGCTCTCCTCTTGGTCGCTCCGGTGGGACATCTGTTTGGCGC‑AGCCGGTGATCTGCTGGGCGCCGCACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTc                                                                                                   <  8:143647/139‑1 (MQ=255)                                                    cAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCAc                                                                                                     <  5:1821/144‑1 (MQ=255)                                                            gCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAACGTACTTCGCTTGTACATGATCACTCGTGCATTGGtt                                                                                        >  4:319051/1‑149 (MQ=255)                                                               cGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTg                                                                                                                >  8:268555/1‑122 (MQ=255)                                                               cGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTg                                                                                                                <  7:268555/122‑1 (MQ=255)                                                                  gCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGt                                                                                                 <  5:213985/134‑1 (MQ=255)                                                                   cAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAgcg                                                                                 <  8:23800/149‑1 (MQ=255)                                                                           gcTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAAcgc                                                                         <  6:126731/149‑1 (MQ=255)                                                                               gACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCgg                                                                       <  8:24300/147‑1 (MQ=255)                                                                               gACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCgg                                                                       <  3:319051/147‑1 (MQ=255)                                                                                   cgcCCCTCCTGGCCGCGGCGCGGGGCCAGCGGGTGGGCGC‑CGCCGCTCACCAGCCGGGCCCCGCACCGGGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGc                                                                 <  6:255232/149‑1 (MQ=255)                                                                                    gcTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCa                                                                <  6:187809/149‑1 (MQ=255)                                                                                        ctcATGGTCGCTGCGCTGGGACATCTGGTTGGCGCTAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAgg                                                              >  1:54638/1‑148 (MQ=255)                                                                                        ctcATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGccct                                                                                                                                       <  4:268583/74‑2 (MQ=255)                                                                                        ctcATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGccct                                                                                                                                       >  3:268583/1‑73 (MQ=255)                                                                                        ctcATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAgg                                                              >  1:120930/1‑147 (MQ=255)                                                                                           aTGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAg                                                           <  2:120930/147‑1 (MQ=255)                                                                                           aTGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGt                                                         >  1:23790/1‑149 (MQ=255)                                                                                                cgctgcgctgGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAg                                                      >  3:82101/1‑147 (MQ=255)                                                                                                cgctgcgctgGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAg                                                      <  4:82101/147‑1 (MQ=255)                                                                                                          ggACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAg                                          >  6:344744/1‑149 (MQ=255)                                                                                                              aTCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCgg                                       >  2:338167/1‑148 (MQ=255)                                                                                                                      ggCGC‑AGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGc                                                                 <  1:162944/114‑1 (MQ=255)                                                                                                                       gcgc‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTg                              <  1:338167/148‑1 (MQ=255)                                                                                                                         gc‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACGTCGCTTGTACATGATCACTCGTGCATTGGGTATCACCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACGTCATTGAGCGGCACCCCCTGGGc                           >  8:224317/1‑149 (MQ=255)                                                                                                                             gCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTc                        <  3:5790/149‑1 (MQ=255)                                                                                                                             gCCGCTGACGAGCAGGGCGCCGCACAGGGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTc                        <  1:169470/149‑1 (MQ=255)                                                                                                                                  tGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCt                   >  6:93192/1‑149 (MQ=255)                                                                                                                                    aCGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcg                 <  5:344744/149‑1 (MQ=255)                                                                                                                                        gCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggtgc             >  5:68970/1‑149 (MQ=255)                                                                                                                                        gCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggtgc             >  1:140319/1‑149 (MQ=255)                                                                                                                                          aGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggtgcgg           <  5:93192/149‑1 (MQ=255)                                                                                                                                             gcgcCGCACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCggtgcggtgc        >  7:45194/1‑149 (MQ=255)                                                                                                                                             gcgcCGCACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCggtgcggtgc        >  5:106196/1‑149 (MQ=255)                                                                                                                                                ccgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTg     <  2:54638/149‑1 (MQ=255)                                                                                                                                                 cgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTg     >  2:142733/1‑148 (MQ=255)                                                                                                                                                  gcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTa   >  3:215147/1‑149 (MQ=255)                                                                                                                                                   cacaGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGGACTCCTCGGAGCGGTGCTTGTAc  >  6:258535/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGC‑AGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTAC  >  NC_002947/1383106‑1383399

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