BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,503,183	C→A	9.7%	L29M (CTG→ATG) ‡	oprC →	copper receptor OprC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,503,183	0	C	A	9.7%	337.0 / 17.4	113	L29M (CTG→ATG) ‡	oprC	copper receptor OprC
Reads supporting (aligned to +/- strand): ref base C (53/49); new base A (5/6); total (58/55)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

Only 100 of 121 total aligned reads displayed.

AACATTTTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCC  >  NC_002947/5503036‑5503332

aacaTTTTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGTATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGTGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGACCATGGCCCTg                                                                                                                                                     >  2:1258689/1‑150 (MQ=255)

aacaTTTTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTg                                                                                                                                                     >  2:2159633/1‑150 (MQ=255)

aacaTTTTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTg                                                                                                                                                     <  2:1811849/150‑1 (MQ=255)

 acaTTTTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTgg                                                                                                                                                    <  2:1703775/150‑1 (MQ=255)

       tGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGccgccg                                                                                                                                              <  1:2159633/150‑1 (MQ=255)

          gcANGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGcc                                                                                                                                           >  2:1500763/1‑150 (MQ=255)

            aCGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCt                                                                                                                                         >  1:1005170/1‑150 (MQ=255)

            aCGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGccag                                                                                                                                         >  1:184859/1‑148 (MQ=255)

               aaGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGc                                                                                                                                      <  1:1058517/150‑1 (MQ=255)

                aGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGcc                                                                                                                                     <  1:883617/150‑1 (MQ=255)

                aGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGcc                                                                                                                                     <  1:2186707/150‑1 (MQ=255)

                    aggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgg                                                                                                                                 <  2:1892396/150‑2 (MQ=255)

                    aggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgg                                                                                                                                 <  2:184859/150‑2 (MQ=255)

                    aggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgg                                                                                                                                 <  1:1924581/150‑2 (MQ=255)

                     ggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCCCGCTCCCATGGCCCTGGCCGCCGGCCCTGGCCATGaa                                                                                                                                <  1:789324/150‑1 (MQ=255)

                         gCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGAcc                                                                                                                            >  2:1941715/1‑150 (MQ=255)

                         gCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGAcc                                                                                                                            >  2:561773/1‑150 (MQ=255)

                           ttCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCAGCCGGCCCTGGCCATGAAGACCAt                                                                                                                          >  1:1188915/1‑150 (MQ=255)

                            tCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtg                                                                                                                         >  1:266892/1‑150 (MQ=255)

                             cgcgCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgt                                                                                                                        <  1:1367559/150‑1 (MQ=255)

                              gcgcTAACTCATGACTGATCGAGTGGAATCATAATGTCCGGCTGCACCCCTGTGTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCACGCTGATCGCGCCCAAGACCCTAACCGCCGCACCAGGCCGTGAGGAACAtgtg                                                                                                                       >  2:1828585/1‑150 (MQ=255)

                               cgcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGc                                                                                                                      <  2:324348/150‑1 (MQ=255)

                                gcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGAGCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCa                                                                                                                     >  1:291612/1‑150 (MQ=255)

                                gcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCa                                                                                                                     >  2:2325434/1‑150 (MQ=255)

                                gcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCa                                                                                                                     >  1:290179/1‑150 (MQ=255)

                                  tAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCacg                                                                                                                   >  2:85800/1‑150 (MQ=255)

                                   aaCTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCacga                                                                                                                  >  1:1219548/1‑150 (MQ=255)

                                    aCTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCacgac                                                                                                                 <  2:105087/150‑1 (MQ=255)

                                     cTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCacgacg                                                                                                                <  2:1831291/150‑1 (MQ=255)

                                      tCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCa                                                                                                                                    >  1:2362734/1‑129 (MQ=255)

                                       cATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGcc                                                                                                              >  1:899061/1‑149 (MQ=255)

                                          gACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCAccc                                                                                                           >  2:1375827/1‑150 (MQ=255)

                                            cTGAACGACTGGAATCATTATGGCCCGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGCGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGACGACCATGTGCACGACGCACCCGa                                                                                                         <  2:1192876/150‑1 (MQ=255)

                                               aTCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGtt                                                                                                      <  1:411128/150‑1 (MQ=255)

                                               aTCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGtt                                                                                                      <  2:557239/150‑1 (MQ=255)

                                                   aCTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTTTTTTTTTTTTTTGTGGGACCATCGCCTCCCTGTTCGGCTCGCTGCTCGCGCCCATTGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGCGTTGTGc                                                                                                  >  2:652588/1‑150 (MQ=255)

                                                     tGGAATCATTATGTCCGGCTGCACGCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGccc                                                                                                >  2:1624479/1‑150 (MQ=255)

                                                          tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGgcc                                                                                                                                                  >  2:1429668/1‑95 (MQ=255)

                                                          tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGgcc                                                                                                                                                  <  1:1429668/95‑1 (MQ=255)

                                                          tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACgg                                                                                           <  2:291612/150‑1 (MQ=255)

                                                          tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACgg                                                                                           <  2:290179/150‑1 (MQ=255)

                                                          tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATg                                                                                                                                  <  1:1226894/111‑1 (MQ=38)

                                                          tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATg                                                                                                                                  >  2:1226894/1‑111 (MQ=38)

                                                           cATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGt                                                                                          <  1:495363/150‑1 (MQ=255)

                                                            aTTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTg                                                                                         >  1:1017448/1‑150 (MQ=255)

                                                             ttATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGa                                                                                        >  1:1447248/1‑150 (MQ=255)

                                                             ttATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGa                                                                                        >  1:1642820/1‑150 (MQ=255)

                                                                tGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCa                                                                                     <  2:266892/150‑1 (MQ=255)

                                                                 gTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGATGATCAc                                                                                    <  1:922377/150‑1 (MQ=255)

                                                                  tCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCac                                                                                                                    >  2:2401924/1‑117 (MQ=255)

                                                                  tCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCac                                                                                                                    <  1:2401924/117‑1 (MQ=255)

                                                                  tCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCAcc                                                                                   <  2:2383851/150‑1 (MQ=255)

                                                                        tGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgt                                                                                                                        >  1:446668/1‑107 (MQ=255)

                                                                        tGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgt                                                                                                                        <  2:446668/107‑1 (MQ=255)

                                                                        tGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGt                                                                              >  1:2078062/1‑149 (MQ=255)

                                                                          cACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGc                                                                           <  1:1375827/150‑1 (MQ=255)

                                                                           aCCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGcc                                                                          >  1:1438110/1‑150 (MQ=255)

                                                                            ccccTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACAGCCGTGGccc                                                                         >  1:702946/1‑150 (MQ=255)

                                                                              ccTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCg                                                                       >  2:2204965/1‑150 (MQ=255)

                                                                               cTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCAt                                                                                                                                   <  2:319756/89‑1 (MQ=255)

                                                                               cTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCAt                                                                                                                                   >  1:319756/1‑89 (MQ=255)

                                                                               cTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGa                                                                      <  1:921611/150‑1 (MQ=255)

                                                                                tGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAg                                                                     <  1:901075/150‑1 (MQ=255)

                                                                                   ttttGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTc                                                                  <  2:2083488/150‑1 (MQ=255)

                                                                                      tGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgcc                                                               <  1:1452212/150‑1 (MQ=255)

                                                                                        ccTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgccgc                                                             >  2:2195891/1‑150 (MQ=255)

                                                                                        ccTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgggcg                                                                                                                              <  1:2217036/85‑5 (MQ=21)

                                                                                         cTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGTCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCt                                                            >  2:461292/1‑150 (MQ=255)

                                                                                         cTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCt                                                            >  1:1067321/1‑150 (MQ=255)

                                                                                         cTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCt                                                            >  2:227987/1‑150 (MQ=255)

                                                                                             tCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACt                                                        <  2:1112476/150‑1 (MQ=255)

                                                                                             tCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACt                                                        >  2:1733142/1‑150 (MQ=255)

                                                                                             tCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACt                                                        <  1:2325434/150‑1 (MQ=255)

                                                                                                cTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACtgtg                                                     >  2:593441/1‑150 (MQ=255)

                                                                                                    gTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCa                                                 >  1:2436804/1‑150 (MQ=255)

                                                                                                     tGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCAc                                                <  2:2238661/150‑1 (MQ=255)

                                                                                                       ggACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCAccc                                              >  2:2368363/1‑149 (MQ=255)

                                                                                                        gACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGCCGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCaa                                             <  1:2369883/150‑1 (MQ=255)

                                                                                                        gACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCaa                                             <  2:201585/150‑1 (MQ=255)

                                                                                                        gACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCaa                                             >  1:190863/1‑150 (MQ=255)

                                                                                                         aCCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCGCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGCGTTGAGCCCGACGGTGATCACCGCCGTGGCCTCGAGTTCGCCGCTGACTGTGGTCACCAAc                                            <  1:1828585/150‑1 (MQ=255)

                                                                                                           cATCGCGGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAAccc                                          >  2:1912186/1‑150 (MQ=255)

                                                                                                               gccgccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGaaa                                      <  1:1290590/150‑1 (MQ=255)

                                                                                                                   cccTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGAcc                                  <  1:1500763/150‑1 (MQ=255)

                                                                                                                   cccTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGAcc                                  >  2:8643/1‑150 (MQ=255)

                                                                                                                              tCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCgccagcca                       <  2:1438110/150‑1 (MQ=255)

                                                                                                                                gctgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCACGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGt                     >  2:40401/1‑150 (MQ=255)

                                                                                                                                    ctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCg                 <  1:593441/150‑1 (MQ=255)

                                                                                                                                    ctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCACCGAGTTCGCCGCTGACTAAAATCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCg                 >  2:1145915/1‑150 (MQ=255)

                                                                                                                                       gcgcCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGcc              >  2:2145499/1‑150 (MQ=255)

                                                                                                                                        tcttcATGGCCAGGGCCGGCGGCCAGGGCCATGGAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACtgtg                                                     >  1:308078/5‑110 (MQ=25)

                                                                                                                                        cgcCCATGGCCCTGGCCGCCGGCCCTGGCCATGCTGTCGAGGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCa             <  2:190863/150‑1 (MQ=255)

                                                                                                                                             gtttcccTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGAc        >  1:83081/5‑150 (MQ=255)

                                                                                                                                              tGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACg       >  2:1109839/1‑150 (MQ=255)

                                                                                                                                              tGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACg       <  1:2145499/150‑1 (MQ=255)

                                                                                                                                                 cccTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgccgc                                                             >  1:434340/1‑93 (MQ=255)

                                                                                                                                                 cccTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgccgc                                                             <  2:434340/93‑1 (MQ=255)

                                                                                                                                                 cccTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgcg    >  2:913380/1‑150 (MQ=255)

                                                                                                                                                  ccTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgcgc   >  2:1998383/1‑150 (MQ=255)

                                                                                                                                                   cTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGcc  >  2:746478/1‑150 (MQ=255)

AACATTTTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCC  >  NC_002947/5503036‑5503332

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: `‑`