BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	514,925	C→G	5.7%	L54L (CTG→CTC) ‡	PP_0426 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	514,925	0	C	G	5.7%	462.3 / 12.9	141	L54L (CTG→CTC) ‡	PP_0426	membrane protein
Reads supporting (aligned to +/- strand): ref base C (70/62); new base G (5/3); total (76/65)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 142 total aligned reads displayed.

ACGCCAGCAGCGACGGGCTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAGGTCCTTGGCGTCG  >  NC_002947/514776‑515072

gtgctagcagcGACGGGCTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATc                                                                                                                                                     >  2:5005083/6‑150 (MQ=255)

   ccAGCAGCGACGGGCTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGc                                                                                                                                                  <  2:595759/150‑1 (MQ=255)

   ccAGCAGCGACGGGCTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGc                                                                                                                                                  >  1:3284795/1‑150 (MQ=255)

           gACGGGCTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGcc                                                                                                                                          >  2:4364760/1‑150 (MQ=255)

           gACGGGCTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGcc                                                                                                                                          >  2:1841127/1‑150 (MQ=255)

              gggCTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCaaa                                                                                                                                       <  1:699009/150‑1 (MQ=255)

              gggCTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCaaa                                                                                                                                       <  1:3893180/150‑1 (MQ=255)

                 cTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTa                                                                                                                                    >  1:1484494/1‑150 (MQ=255)

                 cTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTa                                                                                                                                    <  1:459890/150‑1 (MQ=255)

                  tGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAg                                                                                                                                   >  2:920895/1‑150 (MQ=255)

                      cgcgGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAgg                                                                                                                               <  2:1346070/150‑1 (MQ=255)

                      cgcgGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAgg                                                                                                                               <  2:3925713/150‑1 (MQ=255)

                            gTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACAt                                                                                                                         <  1:1562951/150‑1 (MQ=255)

                              tGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATgg                                                                                                                       <  1:2500322/150‑1 (MQ=255)

                                 aGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGcc                                                                                                                    >  2:2460000/1‑150 (MQ=255)

                                     ggAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCCGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCt                                                                                                                >  2:2323889/1‑150 (MQ=255)

                                      gAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTg                                                                                                               >  2:2913982/1‑150 (MQ=255)

                                        gggCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGaa                                                                                                             >  2:745481/1‑150 (MQ=255)

                                         ggCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCACAGTAGAAGGCAACATGGGCCAGCTGaca                                                                                                            >  2:3714727/1‑148 (MQ=255)

                                         ggCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGaaa                                                                                                            <  2:4017522/150‑1 (MQ=255)

                                           cGAACAAGGTGATCAGCAGGATATAGAGGCCGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCAAGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGc                                                                                                          >  1:12145/1‑150 (MQ=255)

                                             aaCAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTg                                                                                                        <  2:1000482/150‑1 (MQ=255)

                                                aaGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCa                                                                                                     >  2:2524313/1‑150 (MQ=255)

                                                  ggTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACg                                                                                                   >  2:738177/1‑150 (MQ=255)

                                                  ggTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACg                                                                                                   >  1:1466248/1‑150 (MQ=255)

                                                   gTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCCAGGCGCCAAAGTAGAAGGCAACAAGGGCCAGCTGAAAGCTGCCACgt                                                                                                  >  1:1846495/1‑150 (MQ=255)

                                                         agcagGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCg                                                                                            >  2:324510/1‑150 (MQ=255)

                                                          gcagGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCgg                                                                                           >  2:1368431/1‑150 (MQ=255)

                                                          gcagGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCgg                                                                                           <  2:1181861/150‑1 (MQ=255)

                                                           cagGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCggg                                                                                          <  1:2913982/150‑1 (MQ=255)

                                                                 ataGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTg                                                                                    <  2:2991710/150‑1 (MQ=255)

                                                                 ataGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTg                                                                                    <  1:2347517/150‑1 (MQ=255)

                                                                      ggCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCa                                                                               <  1:1208894/150‑1 (MQ=255)

                                                                       gCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAg                                                                              <  2:2675828/150‑1 (MQ=255)

                                                                         tGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGACAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGa                                                                            >  2:812617/1‑150 (MQ=255)

                                                                         tGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGa                                                                            >  1:4358997/1‑150 (MQ=255)

                                                                         tGGAGATGATCAACAGCGCGCCGTCACAGAGGCTCATCCATGCCTCCGTGAGCAACCAGCCCAATTCCCCCATGAAAACCAATGCGCCAAAGAAGAAGGCAACATGGGCCCGCAGAAAGCTGCCACGAGGTGCGGGTTGCTGGCGCAGGa                                                                            >  1:4359008/1‑150 (MQ=255)

                                                                           gCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCATAGTTGAAGGCAACATGGGCCATCTTAAAGCTGCCTCGTGGTGCGGGTTGCTTGCGCATGAAt                                                                          <  2:2634298/150‑1 (MQ=255)

                                                                           gCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAAt                                                                          >  2:3496586/1‑150 (MQ=255)

                                                                              atgatCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGcc                                                                       <  1:1091687/150‑1 (MQ=255)

                                                                                  tCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACa                                                                   >  2:2332263/1‑150 (MQ=255)

                                                                                  tCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACa                                                                   >  2:1396566/1‑150 (MQ=255)

                                                                                    aaCAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAgg                                                                 >  2:4428251/1‑150 (MQ=255)

                                                                                         cgcgcCGTCACCGAGGCTCATCCATGCCTCTGTGAGGAACCAGCGCAGAGCCGCCATGATCAGCACGGCGCCAAAGCAGAAGGCTACTTGGGCCATCGTAAAGCTGCCACGTGTTTCGGGTTGGTGGCGCAGGAATGCCAACAGTCTTTg                                                            <  1:1387243/150‑1 (MQ=255)

                                                                                         cgcgcCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTg                                                            <  1:787518/150‑1 (MQ=255)

                                                                                         cgcgcCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTg                                                            <  1:1387250/150‑1 (MQ=255)

                                                                                          gcgcCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTgg                                                           <  1:5023686/150‑1 (MQ=255)

                                                                                            gcCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTc                                                         <  2:4583788/150‑1 (MQ=255)

                                                                                            gcCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAAACAGCCCATTGCCGCCATGATCAGCATTGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTTGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTc                                                         >  2:1275604/1‑150 (MQ=255)

                                                                                                tCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCtt                                                              >  2:4117790/1‑141 (MQ=255)

                                                                                                tCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCtt                                                              <  1:4117790/141‑1 (MQ=255)

                                                                                                tCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTcctgg                                                     >  1:2486074/1‑150 (MQ=255)

                                                                                                tCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTcctgg                                                     >  1:2486558/1‑150 (MQ=255)

                                                                                                  aCCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTTCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACAAGGGCCAGCTTAAAGATGCCACGTGGTGCGGGTTGCTGGCGCAGGAATTCAAACAGGCTTTGGTcctgtcc                                                   >  1:4718572/1‑150 (MQ=255)

                                                                                                  aCCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTcctggcc                                                   >  1:2550490/1‑150 (MQ=255)

                                                                                                   ccGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTcctggcct                                                  >  1:2099313/1‑150 (MQ=255)

                                                                                                    cGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTcctggcctg                                                 >  2:2609203/1‑150 (MQ=255)

                                                                                                        gCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTg                                             <  2:3326787/150‑1 (MQ=255)

                                                                                                          tCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGga                                           >  1:4852251/1‑150 (MQ=255)

                                                                                                          tCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGga                                           >  2:390092/1‑150 (MQ=255)

                                                                                                          tCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGga                                           >  2:3697881/1‑150 (MQ=255)

                                                                                                          tCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGga                                           <  1:3631271/150‑1 (MQ=255)

                                                                                                          tCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGga                                           <  2:1484494/150‑1 (MQ=255)

                                                                                                             tCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGaa                                        >  1:1915424/1‑150 (MQ=255)

                                                                                                             tCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGaa                                        >  2:1074442/1‑150 (MQ=255)

                                                                                                              atgagCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAAt                                       <  1:3065800/146‑1 (MQ=255)

                                                                                                               cATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGTTCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATg                                      <  1:1926124/150‑1 (MQ=255)

                                                                                                               cATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAAAGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATg                                      >  1:2388258/1‑150 (MQ=255)

                                                                                                                aggCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGc                                     <  2:4930119/148‑1 (MQ=255)

                                                                                                                aTGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGc                                     <  1:920895/150‑1 (MQ=255)

                                                                                                                  gCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCt                                   <  1:2665566/150‑1 (MQ=255)

                                                                                                                    cTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGc                                                                                                                              >  1:4734444/1‑57 (MQ=255)

                                                                                                                    cTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGc                                                                                                                              <  2:4734444/57‑1 (MQ=255)

                                                                                                                     tCGGTGAGCAACCAGCCCATTGCCGCCATGATCGGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTGTGGTCCTGGCCGGGGTGGAGAGTGCCTGCt                                >  2:2022114/1‑150 (MQ=255)

                                                                                                                     tCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCt                                >  1:1875402/1‑150 (MQ=255)

                                                                                                                      cGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCtt                               <  1:2329971/150‑1 (MQ=255)

                                                                                                                      cGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCtt                               <  1:2226212/150‑1 (MQ=255)

                                                                                                                         tGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCa                            >  2:4050368/1‑150 (MQ=255)

                                                                                                                         tGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCa                            >  2:4050182/1‑150 (MQ=255)

                                                                                                                         tGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCa                            >  2:1470420/1‑150 (MQ=255)

                                                                                                                           aGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACg                          >  2:22761/1‑150 (MQ=255)

                                                                                                                           aGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACg                          >  2:22982/1‑150 (MQ=255)

                                                                                                                             cAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGgc                        <  2:3563042/150‑1 (MQ=255)

                                                                                                                             cAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGgc                        <  2:3562825/150‑1 (MQ=255)

                                                                                                                               aCCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGgcgc                      >  1:2982551/1‑150 (MQ=255)

                                                                                                                                ccAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGgcgac                     >  2:1825729/1‑148 (MQ=255)

                                                                                                                                 cAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGgcgcgc                    >  2:2831840/1‑150 (MQ=255)

                                                                                                                                 cAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGgcgcgc                    >  1:4262617/1‑150 (MQ=255)

                                                                                                                                  tctactttGGCGCCATGCTGATCATGGCGGCAATGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCGTGGTTGGAGAATGCCTGCTTGCACGgcgcgcg                   >  2:4182075/7‑150 (MQ=255)

                                                                                                                                   gCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGgcgcgcgc                  >  1:3253106/1‑150 (MQ=255)

                                                                                                                                       tttGGCGCCATGCTGATCATGGCGGCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAgg              >  2:453348/2‑150 (MQ=255)

                                                                                                                                       tttGGCGCCATGCTGATCATGGCGGCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAgg              >  2:244782/2‑150 (MQ=255)

                                                                                                                                        ttGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAGGt             >  2:5119352/1‑150 (MQ=255)

                                                                                                                                           ccgccATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAAt                                                                          <  1:2904889/86‑1 (MQ=255)

                                                                                                                                             gccATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAGGTCCTTg        <  1:324510/150‑1 (MQ=255)

                                                                                                                                             gccATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAGGTCCTTg        <  2:2669776/150‑1 (MQ=255)

                                                                                                                                               cATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAGGTCCTTGGc      <  2:4333951/150‑1 (MQ=255)

                                                                                                                                                  gATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAGGTCCTTGGCGTc   <  1:687144/150‑1 (MQ=255)

                                                                                                                                                  gATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAGGTCCTTGGCGTc   <  1:80335/150‑1 (MQ=255)

                                                                                                                                                   aTCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAGGTCCTTGGCGTCg  >  1:3669359/1‑150 (MQ=255)

ACGCCAGCAGCGACGGGCTGGACGCGGCGTTGCAGATGGAGGGCGAACAAGGTGATCAGCAGGATATAGAGGCTGGCGATGATCAACAGCGCGCCGTCACCGAGGCTCATCCATGCCTCGGTGAGCAACCAGCCCATTGCCGCCATGATCAGCATGGCGCCAAAGTAGAAGGCAACATGGGCCAGCTGAAAGCTGCCACGTGGTGCGGGTTGCTGGCGCAGGAATGCCAACAGGCTTTGGTCCTGGCCTGGTTGGAGAATGCCTGCTTGCACGGCGCGCGCCAGGTCCTTGGCGTCG  >  NC_002947/514776‑515072

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`