BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,166,929:1	+C	5.7%	coding (591/1212 nt)	creA ←	creatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,166,929	1	.	C	5.7%	435.3 / 18.6	106	coding (591/1212 nt)	creA	creatinase
Reads supporting (aligned to +/- strand): ref base . (61/39); new base C (2/4); total (63/43)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 109 total aligned reads displayed.

CTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGA  >  NC_002947/4166780‑4167077

cTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCa                                                                                                                                                       <  1:1581158/150‑1 (MQ=255)

    ggATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCa                                                                                                                                                  <  1:1361342/150‑1 (MQ=255)

        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)

              cTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTc                                                                                                                                        <  1:1276614/150‑1 (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)

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

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

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

                                      cGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAg                                                                                                                >  2:2241524/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                                                                                      >  1:2327539/1‑150 (MQ=255)

                                                                gCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGa                                                                                      >  2:847712/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:1442571/1‑126 (MQ=255)

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

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

                                                                    ggACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTc                                                                                                          <  2:1442571/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                                                                              >  1:1837136/1‑150 (MQ=255)

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

                                                                        tGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGcc                                                                              >  2:554274/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)

                                                                             ccaggtccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCtg                                                                         >  2:1576388/1‑150 (MQ=255)

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

                                                                                 gtccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCACTGCA‑TGCCAccca                                                                                                                                        <  2:1221140/83‑3 (MQ=255)

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

                                                                                  tccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCAc                                                                    >  2:380146/1‑150 (MQ=255)

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

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

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

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

                                                                                        tgtCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGGGCACCGTGAc                                                              <  2:2435782/150‑1 (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)

                                                                                            cATCAGCTCGGCGTCGGGGAAGCGTTGGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGCCGGAt                                                          >  2:701592/1‑150 (MQ=255)

                                                                                              tCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATaa                                                        >  1:2310167/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                                                                                                                                                 <  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)

                                                                                                                   gTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACt                                   <  1:301777/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)

                                                                                                                              tCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATgcgc                        <  1:847712/150‑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)

                                                                                                                                      cATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCAt                <  2:1837136/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                                                               <  1:2098016/99‑1 (MQ=255)

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

                                                                                                                                             cgaggTGGCACTGCA‑TGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACgg          >  2:844659/4‑150 (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)

CTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGC‑CTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGA  >  NC_002947/4166780‑4167077

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