BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,166,934	Δ1 bp	5.5%	coding (586/1212 nt)	creA ←	creatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,166,934	0	G	.	5.5%	338.2 / 14.1	109	coding (586/1212 nt)	creA	creatinase
Reads supporting (aligned to +/- strand): ref base G (63/40); new base . (2/4); total (65/44)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

Only 100 of 114 total aligned reads displayed.

CTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCA  >  NC_002947/4166788‑4167083

ctcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGc                                                                                                                                                    <  1:1981730/150‑1 (MQ=255)

     ccTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCt                                                                                                                                               >  2:1894918/1‑150 (MQ=255)

      cTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCACTGCA‑TGCCAccca                                                                                                                                              <  2:1435587/150‑3 (MQ=255)

      cTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGCCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCACTGCA‑TGCCAccca                                                                                                                                              <  1:1094122/150‑3 (MQ=255)

       ttGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCg                                                                                                                                             >  2:301777/1‑150 (MQ=255)

       ttGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGCTGCACGACGCTCTTCCGATCCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCACTGCA‑TGCCAcccag                                                                                                                                             <  2:1242370/150‑4 (MQ=255)

         gCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTa                                                                                                                                           <  1:2026350/150‑1 (MQ=255)

             acacGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCg                                                                                                                                       >  1:2098765/1‑150 (MQ=255)

              cacGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCgg                                                                                                                                      <  2:1724386/150‑1 (MQ=255)

                          tGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCAc                                                                                                                          >  1:850179/1‑150 (MQ=255)

                           gACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCCCG‑AGCATTTCTCGGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCCCCTGCTCACg                                                                                                                         >  2:663824/1‑150 (MQ=255)

                           gACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCAACAGCGCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCGCGTACTCGGGCACCTGCTcgcg                                                                                                                         >  1:424911/1‑150 (MQ=255)

                              cGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGt                                                                                                                                                      >  2:152474/1‑118 (MQ=255)

                              cGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAg                                                                                                                      <  2:975875/150‑1 (MQ=255)

                              cGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAg                                                                                                                      >  2:160085/1‑150 (MQ=255)

                              cGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAg                                                                                                                      >  2:2241524/1‑150 (MQ=255)

                              cGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAg                                                                                                                      >  1:1794003/1‑150 (MQ=255)

                              cGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAg                                                                                                                      >  1:1381906/1‑150 (MQ=255)

                                   tgtgtgCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCtt                                                                                                                 <  1:2413965/150‑1 (MQ=255)

                                    gtgtgCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTc                                                                                                                <  1:2421755/150‑1 (MQ=255)

                                    gtgtgCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTc                                                                                                                <  1:2248088/150‑1 (MQ=255)

                                       tgCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGAc                                                                                                             >  1:1665972/1‑150 (MQ=255)

                                         cgcCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCa                                                                                                           <  2:475111/150‑1 (MQ=255)

                                            cATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCGTCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTg                                                                                                        >  2:2180482/1‑150 (MQ=255)

                                                        gCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGa                                                                                            <  2:1794003/150‑1 (MQ=255)

                                                        gCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGa                                                                                            >  2:847712/1‑150 (MQ=255)

                                                        gCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGa                                                                                            >  1:2327539/1‑150 (MQ=255)

                                                         cccGGACTGGAACCAGGTCCAGGTGTCTATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGAt                                                                                           >  1:917092/1‑150 (MQ=255)

                                                           cGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGa                                                                                                              >  1:1080278/1‑129 (MQ=255)

                                                           cGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGa                                                                                                              <  2:1080278/129‑1 (MQ=255)

                                                            ggACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTc                                                                                                                >  1:1447862/1‑126 (MQ=255)

                                                            ggACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTc                                                                                                                >  1:1442571/1‑126 (MQ=255)

                                                            ggACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTc                                                                                                                <  2:1442571/126‑1 (MQ=255)

                                                            ggACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTc                                                                                                                <  2:1447862/126‑1 (MQ=255)

                                                             gACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCg                                                                                       <  2:1792773/150‑1 (MQ=255)

                                                               cTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGc                                                                                     >  2:1113209/1‑150 (MQ=255)

                                                                tGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGcc                                                                                    >  2:554274/1‑150 (MQ=255)

                                                                tGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGcc                                                                                    >  2:558209/1‑150 (MQ=255)

                                                                  gAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGNGTTTGGNAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCAc                                                                                  >  1:2173974/1‑150 (MQ=255)

                                                                   aaCCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCa                                                                                                           >  2:27436/1‑124 (MQ=255)

                                                                   aaCCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCa                                                                                                           <  1:27436/124‑1 (MQ=255)

                                                                      caggtccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCtgt                                                                              >  2:1823068/1‑150 (MQ=255)

                                                                         gtccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCACTGCA‑TGCCAccca                                                                                                                                              >  1:1221140/1‑81 (MQ=255)

                                                                               gtgtCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGa                                                                     >  1:1855454/1‑150 (MQ=255)

                                                                               gtgtCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGa                                                                     >  2:2452572/1‑150 (MQ=255)

                                                                                tgtCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGAc                                                                    <  1:1894918/150‑1 (MQ=255)

                                                                                tgtCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGAc                                                                    >  2:811544/1‑150 (MQ=255)

                                                                                tgtCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGGGCACCGTGAc                                                                    <  2:2435782/150‑1 (MQ=255)

                                                                                 gtCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGGTTCGACCACTGCCGCCCCCCCGAAGGCGGCGCGCTGTGCACCGTGCCg                                                                   >  2:552608/1‑150 (MQ=255)

                                                                                 gtCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACg                                                                   >  2:1193581/1‑150 (MQ=255)

                                                                                  tCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACgg                                                                  >  1:640301/1‑150 (MQ=255)

                                                                                   ccATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCNA‑ATGAGAAGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCGCTCCCGCCCAGCCGATGTCGGCCACCGGTGCACCGTGACGGa                                                                 >  1:300700/1‑150 (MQ=255)

                                                                                    cATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGAt                                                                >  2:1921282/1‑150 (MQ=255)

                                                                                    cATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCAACTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGAt                                                                >  1:604599/1‑150 (MQ=255)

                                                                                      tCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATaa                                                              >  2:1307373/1‑150 (MQ=255)

                                                                                        aGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCAGCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGt                                                            <  1:178273/150‑1 (MQ=255)

                                                                                             ggCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTg                                                       >  2:945758/1‑150 (MQ=255)

                                                                                              gCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGc                                                      <  2:2310167/150‑1 (MQ=255)

                                                                                                   ggggAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGctcctc                                                 <  1:2293054/150‑1 (MQ=255)

                                                                                                   ggggAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGctcctc                                                 >  2:908824/1‑150 (MQ=255)

                                                                                                   ggggAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGctcctc                                                 >  1:176649/1‑150 (MQ=255)

                                                                                                   ggggAAGCGTTTGGCAATTTCGCACACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGcccctc                                                 >  2:1818088/1‑150 (MQ=255)

                                                                                                    gggAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAg                                                                                                                                                       >  2:1130724/1‑47 (MQ=255)

                                                                                                    gggAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAg                                                                                                                                                       <  1:1130724/47‑1 (MQ=255)

                                                                                                    gggAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCg                                                <  1:2452572/150‑1 (MQ=255)

                                                                                                    gggAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCg                                                <  1:1193581/150‑1 (MQ=255)

                                                                                                     ggAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTCGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCgg                                               <  1:1307373/150‑1 (MQ=255)

                                                                                                            tttGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACtt                                        <  1:1113209/150‑1 (MQ=255)

                                                                                                              tGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTCCTCGGCCCCCTGCCCACGCACGGCGGGGTGCACTGCCGCTCCGCTGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGa                                      <  1:2120485/150‑1 (MQ=255)

                                                                                                               ggCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGAt                                     >  1:2366992/1‑150 (MQ=255)

                                                                                                                gCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATc                                    >  1:459900/1‑150 (MQ=255)

                                                                                                                  aaTTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCAt                                  >  1:1369857/1‑150 (MQ=255)

                                                                                                                  aaTTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCAt                                  >  2:2368594/1‑150 (MQ=255)

                                                                                                                   aTTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATg                                 <  2:459900/150‑1 (MQ=255)

                                                                                                                    tttCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATgc                                >  1:506706/1‑150 (MQ=255)

                                                                                                                    tttCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACACAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGAt                                                                >  1:835316/1‑118 (MQ=255)

                                                                                                                    tttCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACACAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGAt                                                                <  2:835316/118‑1 (MQ=255)

                                                                                                                          gcACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTACTCAAGCTTGACGTCGACCAATGGCACCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATAc                          >  1:1486275/1‑150 (MQ=255)

                                                                                                                           cACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCTGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCAt                                  <  1:99489/141‑1 (MQ=255)

                                                                                                                           cACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACg                         >  1:2083014/1‑150 (MQ=255)

                                                                                                                              cATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCAt                      <  2:2083014/150‑1 (MQ=255)

                                                                                                                                tCGC‑CTGGGTGGCG‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATAc                          <  2:1591180/144‑1 (MQ=255)

                                                                                                                                tCGC‑CTGGGTGGCG‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATAc                          >  1:1591180/1‑144 (MQ=255)

                                                                                                                                tCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGc                    >  1:1568411/1‑150 (MQ=255)

                                                                                                                                  gC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGa                                                                     >  2:2098016/1‑99 (MQ=255)

                                                                                                                                     cgaggTGGCACTGCA‑TGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACgg                >  2:844659/4‑150 (MQ=255)

                                                                                                                                     cTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCg                                                >  2:2225788/1‑118 (MQ=255)

                                                                                                                                     cTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCg                                                <  1:2225788/118‑1 (MQ=255)

                                                                                                                                     cTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACgg                <  1:1679156/150‑1 (MQ=255)

                                                                                                                                          tGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAg           >  1:1926921/1‑150 (MQ=255)

                                                                                                                                            gCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCg         >  1:1090907/1‑150 (MQ=255)

                                                                                                                                             cA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACAGAAAAGTGCCTGCTGCTTGGAACACGTGAGCTTCGGCACAAGCATGGACGGCCCGGCGa        >  2:519862/1‑150 (MQ=255)

                                                                                                                                              a‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGTATGCGCATgatcaagt              <  2:128532/143‑8 (MQ=255)

                                                                                                                                              a‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGTATGCGCATgatcaagt              >  1:128532/1‑136 (MQ=255)

                                                                                                                                                   aGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGAc                                          <  1:2078100/111‑1 (MQ=255)

                                                                                                                                                   aGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGAc                                          >  2:2078100/1‑111 (MQ=255)

                                                                                                                                                   aGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGc          >  1:1111484/1‑143 (MQ=255)

                                                                                                                                                   aGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGc          <  2:1111484/143‑1 (MQ=255)

                                                                                                                                                    gTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGt     <  1:2227849/147‑1 (MQ=255)

                                                                                                                                                    gTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTcca  <  1:589087/150‑1 (MQ=255)

CTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCA  >  NC_002947/4166788‑4167083

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