BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,248,398	C→T	7.1%	R558H (CGC→CAC)	PP_3724 ←	acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,248,398	0	C	T	7.1%	285.9 / 8.8	98	R558H (CGC→CAC)	PP_3724	acyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base C (57/34); new base T (4/3); total (61/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 101 total aligned reads displayed.

GCACGCCGCCGGCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGC  >  NC_002947/4248249‑4248542

gCACGCCGCCGGCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGgc                                                                                                                                                  <  2:2216808/150‑1 (MQ=255)

          ggCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCAc                                                                                                                                        <  2:194770/150‑1 (MQ=255)

             aCCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAg                                                                                                                                     >  1:1268218/1‑150 (MQ=255)

              ccGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGt                                                                                                                                    <  1:613874/150‑1 (MQ=255)

              ccGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGt                                                                                                                                    <  1:838183/150‑1 (MQ=255)

               cGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTg                                                                                                                                   >  2:816683/1‑150 (MQ=255)

               cGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTg                                                                                                                                   >  1:1501761/1‑150 (MQ=255)

                    aCCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGTGCat                                                                                                                                >  1:1142890/1‑147 (MQ=255)

                    aCCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGTGCat                                                                                                                                <  2:1142890/148‑2 (MQ=255)

                      cAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGTCGGCCCGGCCATGTTCGCGAGCAGCACCTTATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTc                                                                                                                            >  2:2103781/1‑150 (MQ=255)

                      cAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTc                                                                                                                            >  2:1948419/1‑150 (MQ=255)

                      cAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTc                                                                                                                            >  1:1523237/1‑150 (MQ=255)

                       aaGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGgc                                                                                                                                                  >  1:1070614/1‑127 (MQ=255)

                       aaGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGgc                                                                                                                                                  <  2:1070614/127‑1 (MQ=255)

                        aGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTAGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGTGCatcgctga                                                                                                                          <  1:2036771/150‑8 (MQ=255)

                        aGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGAGCACCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCACCTGCACCAGTGCAGCTTCCa                                                                                                                          <  1:1206815/150‑1 (MQ=255)

                           ccAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGGCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAggg                                                                                                                       >  1:469690/1‑150 (MQ=255)

                           ccAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCGCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAggg                                                                                                                       <  2:469690/150‑1 (MQ=255)

                                       tgttTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCANcgc                                                                                                           <  2:1501761/150‑1 (MQ=255)

                                         ttttCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCgcgc                                                                                                         >  1:1169707/1‑150 (MQ=255)

                                         ttttCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCgcgc                                                                                                         <  2:418875/150‑1 (MQ=255)

                                           ttCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTACGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAg                                                                                                       >  1:1256162/1‑150 (MQ=255)

                                             cACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGt                                                                                                     >  2:181401/1‑150 (MQ=255)

                                              aCGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGtt                                                                                                    >  1:320055/1‑150 (MQ=255)

                                                gCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGCGCATCGCGAAGGTTCa                                                                                                  >  2:1224190/1‑150 (MQ=255)

                                                   tGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACg                                                                                               <  1:1510612/150‑1 (MQ=255)

                                                     aTGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCt                                                                                                                                            <  2:1556914/103‑1 (MQ=255)

                                                     aTGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCt                                                                                                                                            >  1:1556914/1‑103 (MQ=255)

                                                     aTGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGcc                                                                                             <  2:1129445/150‑1 (MQ=255)

                                                         ccACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGc                                                                                                                                  <  2:2273676/109‑1 (MQ=255)

                                                         ccACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGc                                                                                                                                  >  1:2273676/1‑109 (MQ=255)

                                                         ccACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAg                                                                                         >  1:113958/1‑150 (MQ=255)

                                                               tGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAg                                                                                                       >  1:692326/1‑130 (MQ=255)

                                                               tGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAg                                                                                                       <  2:692326/130‑1 (MQ=255)

                                                                  gCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCtt                                                                                <  2:1680593/150‑1 (MQ=255)

                                                                  gCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCtt                                                                                <  2:1594241/150‑1 (MQ=255)

                                                                  gCTCGGCCAAGGGCGGCTCGGCCATGGACTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCtt                                                                                >  2:2211943/1‑150 (MQ=255)

                                                                   cTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCTCAGGTTCAACGCCACAGCACCAGCTTc                                                                               <  1:1242661/150‑1 (MQ=255)

                                                                   cTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTc                                                                               >  2:119068/1‑150 (MQ=255)

                                                                             ggcggcTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGc                                                                     <  1:181401/150‑1 (MQ=255)

                                                                                ggcTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCCTCCCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGc                                                                  >  1:959231/1‑150 (MQ=255)

                                                                                ggcTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGc                                                                  >  1:2198937/1‑150 (MQ=255)

                                                                                 gcTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCa                                                                 >  1:855841/1‑150 (MQ=255)

                                                                                   tCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGc                                                               >  2:1777474/1‑150 (MQ=255)

                                                                                   tCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGc                                                               >  1:1777856/1‑150 (MQ=255)

                                                                                      gCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTGGTGCGGCAGCTgc                                                            >  2:1292976/1‑150 (MQ=255)

                                                                                       ccATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACCGCACCAGCTTCGGTCTTGTGCGGCAGCTgcg                                                           >  2:2135344/1‑150 (MQ=255)

                                                                                       ccATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTgcg                                                           >  2:1379151/1‑150 (MQ=255)

                                                                                       ccATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTgcg                                                           <  2:1169707/150‑1 (MQ=255)

                                                                                        cATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGAGCGTAGACAACGGTGCACAAGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCAACCAGCAGATGCAACGCCACAGCACCAGCTTCGGACCTGTACGGCATCTgcgc                                                          >  2:1046614/1‑150 (MQ=255)

                                                                                             gCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGCCTCAACGCCACAGCACCAGCTGCGGTCTTTTGCGGCAGCTGCGCGCTGa                                                     >  2:1104227/1‑150 (MQ=255)

                                                                                               tCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTGGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGAcc                                                   <  1:1948419/150‑1 (MQ=255)

                                                                                               tCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCa                                                                                                                                 <  2:187313/72‑1 (MQ=255)

                                                                                               tCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCa                                                                                                                                 >  1:187313/1‑72 (MQ=255)

                                                                                                  agcagcACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCCCACTCGCGTCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCGCAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGcc                                                >  1:1240912/1‑150 (MQ=255)

                                                                                                  agcagcACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGcc                                                <  1:1292976/150‑1 (MQ=255)

                                                                                                   gcagcaCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCt                                               >  1:232932/1‑150 (MQ=255)

                                                                                                    cagcaCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCtt                                              <  1:2211943/150‑1 (MQ=255)

                                                                                                       caCCTGATCAAACGCAACCTACGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAg                                           <  2:1256162/150‑1 (MQ=255)

                                                                                                        aCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACGGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGt                                          >  1:1368462/1‑150 (MQ=255)

                                                                                                        aCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCa                                            <  2:680873/148‑1 (MQ=255)

                                                                                                        aCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCa                                            >  1:680873/1‑148 (MQ=255)

                                                                                                        aCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGt                                          >  2:2187441/1‑150 (MQ=255)

                                                                                                        aCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGAATCCTCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGt                                          >  1:2252829/1‑150 (MQ=255)

                                                                                                         ccTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTa                                         >  2:651292/1‑150 (MQ=255)

                                                                                                          cTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTAc                                        <  2:113958/150‑1 (MQ=255)

                                                                                                          cTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTAc                                        <  1:651292/150‑1 (MQ=255)

                                                                                                             aTCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTTCGGCAGCTGCGCGCTGACCGCCTTCAGTACCaa                                     <  1:1379151/150‑1 (MQ=255)

                                                                                                             aTCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTgc                                                            >  2:1175883/1‑127 (MQ=255)

                                                                                                             aTCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTgc                                                            <  1:1175883/127‑1 (MQ=255)

                                                                                                              tCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCATGGCGCGGGCATCGCGCAGGTTCACCGCCACAGCCCCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAg                                    >  2:718157/1‑150 (MQ=255)

                                                                                                              tCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAg                                    >  2:718145/1‑150 (MQ=255)

                                                                                                              tCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAg                                    >  2:718146/1‑150 (MQ=255)

                                                                                                               cAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAgg                                   <  2:464469/150‑1 (MQ=255)

                                                                                                                aaaCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                  >  1:1821101/1‑150 (MQ=255)

                                                                                                                       aCCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCtt                                              >  2:1801190/1‑131 (MQ=255)

                                                                                                                       aCCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCtt                                              <  1:1801190/131‑1 (MQ=255)

                                                                                                                              gggCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACgcggg                    >  1:371576/1‑150 (MQ=255)

                                                                                                                              gggCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACgcggg                    >  1:753242/1‑150 (MQ=255)

                                                                                                                              gggCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACgcggg                    >  1:1477384/1‑150 (MQ=255)

                                                                                                                               ggCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTAAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACgcgggc                   >  1:714130/1‑150 (MQ=255)

                                                                                                                                gCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACgcgggcg                  <  1:2109511/150‑1 (MQ=255)

                                                                                                                                  gTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGggcggc                <  2:1523237/150‑1 (MQ=255)

                                                                                                                                     gTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGAcc                                                   <  1:445790/112‑1 (MQ=255)

                                                                                                                                     gTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGAcc                                                   >  2:445790/1‑112 (MQ=255)

                                                                                                                                      tCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCt            >  1:2070496/1‑150 (MQ=255)

                                                                                                                                      tCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGAAACTCTCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCt            >  1:1968987/1‑150 (MQ=255)

                                                                                                                                        aGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGc          >  1:710911/1‑150 (MQ=255)

                                                                                                                                           gATGCACTGTCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCca                                                                                            <  1:1485098/65‑1 (MQ=255)

                                                                                                                                           gATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCca                                                                                            >  2:1485098/1‑65 (MQ=255)

                                                                                                                                            aTGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACg      >  1:2141775/1‑150 (MQ=255)

                                                                                                                                              gCACTGGTGCAGATGCGCCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGgc    >  1:1054174/1‑150 (MQ=255)

                                                                                                                                                aCTGGTGCAGATGCGCCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGgcgc  >  2:602936/1‑150 (MQ=255)

                                                                                                                                                aCTGGTGCAGATGCGCCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGCCTTGGGCGTCAGCTTCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGgcgc  >  1:288657/1‑150 (MQ=255)

                                                                                                                                                aCTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCAc                                                                                      <  2:762604/66‑1 (MQ=255)

                                                                                                                                                aCTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCAc                                                                                      >  1:762604/1‑66 (MQ=255)

                                                                                                                                                aCTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCa        <  1:943192/144‑1 (MQ=255)

                                                                                                                                                aCTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCa        >  2:943192/1‑144 (MQ=255)

                                                                                                                                                aCTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGCAACCGACGCGGGCGGCGTCTGCCACGgcgc  >  1:603512/1‑150 (MQ=255)

                                                                                                                                                aCTGGCGCATCTGCACCAGTGCAGCTACAAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGgcgc  >  1:437362/1‑150 (MQ=255)

GCACGCCGCCGGCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGC  >  NC_002947/4248249‑4248542

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