BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,875,934	G→C	6.3%	M359I (ATG→ATC)	PP_4284 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,875,934	0	G	C	6.3%	301.2 / 7.8	96	M359I (ATG→ATC)	PP_4284	transporter
Reads supporting (aligned to +/- strand): ref base G (38/51); new base C (4/2); total (43/53)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 102 total aligned reads displayed.

TGTTCGGCGCGGTGGTCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGA  >  NC_002947/4875786‑4876083

tgtTCGGCGCGGTGGTCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGgatga                                                                                                                                                       >  1:236504/1‑150 (MQ=255)

tgtTCGGCGCGGTGGTCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGgatga                                                                                                                                                       <  1:1616943/150‑1 (MQ=255)

 gtTCGGCGCGGTGGTCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGgatgat                                                                                                                                                      >  1:1656055/1‑150 (MQ=255)

  ttCGGCGCGGTGGTCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGAtt                                                                                                                                                     >  1:448653/1‑150 (MQ=255)

       cgcgGTGGTCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGTCGGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGc                                                                                                                                               <  2:74141/150‑1 (MQ=255)

       cgcgGTGGTCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGc                                                                                                                                               <  1:358622/150‑1 (MQ=255)

          ggtggtCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTa                                                                                                                                            >  1:570596/1‑150 (MQ=255)

           gtggtCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTAt                                                                                                                                           <  2:1413361/150‑1 (MQ=255)

              gtCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGgatga                                                                                                                                                       <  1:1139803/136‑1 (MQ=255)

              gtCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGgatga                                                                                                                                                       >  2:1139803/1‑136 (MQ=255)

                 ggTGAACCGCCGGTGACCAGCAAGGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGAca                                                                                                                                     >  2:1133632/1‑150 (MQ=255)

                  gtgtACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGAcag                                                                                                                                    <  2:218673/150‑1 (MQ=255)

                     tACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGAcagcag                                                                                                                                 <  2:75662/150‑1 (MQ=255)

                     tACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGAcagcag                                                                                                                                 <  2:1662276/150‑1 (MQ=255)

                       ccgccgGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGGAGGCCTGCTGTTCATCGCCGCGATCTTATTCGCGCCGCGGGCTTGGATGATTCC‑GGCCTATGCGACAGCAGGt                                                                                                                               >  1:1738720/1‑150 (MQ=255)

                          ccgGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCa                                                                                                                            <  1:811233/150‑1 (MQ=255)

                            ggTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACt                                                                                                                          <  1:1725095/150‑1 (MQ=255)

                             gTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTTGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTg                                                                                                                         >  2:1141315/1‑150 (MQ=255)

                                  cAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGAcagca                                                                                                                                  >  1:2270596/1‑136 (MQ=255)

                                  cAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGAcagca                                                                                                                                  <  2:2270596/136‑1 (MQ=255)

                                   aGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGATGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTAc                                                                                                                   >  1:14505/1‑150 (MQ=255)

                                     cTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGt                                                                                                                 >  2:402344/1‑150 (MQ=255)

                                          tgGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGa                                                                                                            >  2:1405992/1‑150 (MQ=255)

                                            gagaGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATg                                                                                                          >  1:1204856/1‑150 (MQ=255)

                                            gagaGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATg                                                                                                          >  2:978441/1‑150 (MQ=255)

                                             agagTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGc                                                                                                         <  1:978441/150‑1 (MQ=255)

                                                 tgctgcGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCtgat                                                                                                     <  1:1653282/150‑1 (MQ=255)

                                                      cGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATgg                                                                                                <  2:1212040/150‑1 (MQ=255)

                                                        ggggTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGgc                                                                                              <  1:777180/150‑1 (MQ=255)

                                                            ttGCCGCAGGTGGCCGAACTGGCCTGACGGCTATGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGgcagca                                                                                          >  1:1068618/1‑150 (MQ=255)

                                                             tGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCAt                                                                                         <  1:1568620/150‑1 (MQ=255)

                                                             tGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCAt                                                                                         <  1:1568667/150‑1 (MQ=255)

                                                              gCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATg                                                                                        <  1:1812142/150‑1 (MQ=255)

                                                              gCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATg                                                                                        <  1:1156159/150‑1 (MQ=255)

                                                               ccGCAGGTGGCCGAACTGGCCTGACGGCTATGGTGGTAGTCCTGCTGTTCATCTCCGCGCTGTTCTTCGCGTCGCTGGCATAG‑TACTGCCTGTGCTGAGCTACAGCAGGTGCACTGATCTACGTTGCGATGCTAATGATGGGCAGCATgg                                                                                       >  1:551234/1‑150 (MQ=255)

                                                                cGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGgc                                                                                      >  1:1753194/1‑150 (MQ=255)

                                                                cGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATG‑TCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGc                                                                                                                             <  1:937840/110‑1 (MQ=255)

                                                                cGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATG‑TCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGc                                                                                                                             >  2:937840/1‑110 (MQ=255)

                                                                 gCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGgcg                                                                                     >  1:869285/1‑150 (MQ=255)

                                                                   aGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGca                                                                                   >  2:100221/1‑150 (MQ=255)

                                                                            aaCTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTg                                                                          <  1:576819/150‑1 (MQ=255)

                                                                            aaCTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGATGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTg                                                                          <  2:14505/150‑1 (MQ=255)

                                                                            aaCTGGACTGACGGCTGTGGTATCACCCTTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGTTAGGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTg                                                                          <  1:1133632/150‑1 (MQ=255)

                                                                             aCTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCAATGCTGATGATGGGCAGCATGGCGCACATCCATTgg                                                                         <  1:1085511/150‑1 (MQ=255)

                                                                              cTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTggg                                                                        <  1:41759/150‑1 (MQ=255)

                                                                               tGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGa                                                                       <  1:402344/150‑1 (MQ=255)

                                                                                ggCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGAt                                                                      <  2:1050476/150‑1 (MQ=255)

                                                                                   cTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAg                                                                   <  2:920788/150‑1 (MQ=255)

                                                                                   cTGACGGCTGTGGTGGTAGGCCTGCTGTGCATCGCCGCGGTGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCGGATGATGGGCAGCATGGCGCACATCCATTGGGCTGAg                                                                   >  2:1631487/1‑150 (MQ=255)

                                                                                   cTGACGGCTGTGGGGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAg                                                                   >  2:478633/1‑150 (MQ=255)

                                                                                        tgCTGCGGTGGTCGGCCTGCAGTTCATTGCCGCGATGTTATGCGCGCCGCTGGCCGGAATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCGACGTTGCGATGCTGATGATCCGCAGCATGGCGCACATCCATTGGGATGAGGCCAc                                                              <  2:417840/149‑1 (MQ=255)

                                                                                           tgtgTTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCtgatga                                                                                                   >  1:865050/1‑110 (MQ=255)

                                                                                           tgtgTTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCtgatga                                                                                                   <  2:865050/110‑1 (MQ=255)

                                                                                           tgtgGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGAt                                                                                                                       >  2:1815241/1‑90 (MQ=255)

                                                                                           tgtgGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGAt                                                                                                                       <  1:1815241/90‑1 (MQ=255)

                                                                                           tgtgGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCtg                                                                                                       >  1:916636/1‑106 (MQ=255)

                                                                                           tgtgGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCtg                                                                                                       <  2:916636/106‑1 (MQ=255)

                                                                                             tggtggTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACa                                                         <  2:236504/150‑1 (MQ=255)

                                                                                             tggtggTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCAATTGGGATGAGGCCACCGaga                                                         <  2:448653/150‑3 (MQ=255)

                                                                                                tggtAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCa                                                      >  2:2381354/1‑150 (MQ=255)

                                                                                                    aGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTcc                                                  >  1:205396/1‑150 (MQ=255)

                                                                                                      gCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCAtt                                                                           <  1:707359/123‑1 (MQ=255)

                                                                                                      gCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCAtt                                                                           >  2:707359/1‑123 (MQ=255)

                                                                                                       ccTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGc                                               >  1:27078/1‑150 (MQ=255)

                                                                                                         tgctgTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGa                                             >  2:1787926/1‑150 (MQ=255)

                                                                                                           cgagtCATCGCCGCGATGCTCTTCGCGCCGCTGGTCGGGATGAGTGC‑GGCCTATGCTACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATc                                           <  2:1068618/146‑1 (MQ=255)

                                                                                                            tgTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCg                                          >  1:1585570/1‑150 (MQ=255)

                                                                                                             gTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCtgat                                                                                                     <  1:1549396/90‑1 (MQ=255)

                                                                                                             gTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCtgat                                                                                                     >  2:1549396/1‑90 (MQ=255)

                                                                                                                 aTCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCtgatg                                                                                                    >  1:667754/1‑87 (MQ=255)

                                                                                                                 aTCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCtgatg                                                                                                    <  2:667754/87‑1 (MQ=255)

                                                                                                                 aTCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCtga                                                                                                      <  1:987545/85‑1 (MQ=255)

                                                                                                                 aTCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTca                                       <  2:479164/148‑1 (MQ=255)

                                                                                                                 aTCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTca                                       >  1:479164/1‑148 (MQ=255)

                                                                                                                  tCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAg                                    >  2:1023975/1‑150 (MQ=255)

                                                                                                                      cgcgATGTTCTTGGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGAt                                <  2:156992/150‑1 (MQ=255)

                                                                                                                      cgcgATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCCGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGAt                                >  1:1378506/1‑150 (MQ=255)

                                                                                                                         gATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGGGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCat                             <  2:570596/150‑1 (MQ=255)

                                                                                                                         gATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCat                             <  1:774149/150‑1 (MQ=255)

                                                                                                                           tGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCatga                           <  1:1787926/150‑1 (MQ=255)

                                                                                                                                 tcGCGCCGCTGGCCGGGATCATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGc                     <  2:1585570/150‑1 (MQ=255)

                                                                                                                                  cgcgCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCt                    <  2:1738720/150‑1 (MQ=255)

                                                                                                                                   gcgcCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTg                   >  2:2177861/1‑150 (MQ=255)

                                                                                                                                   gcgcCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTg                   >  1:1226966/1‑150 (MQ=255)

                                                                                                                                    cgccgcTGGCCGGGATGATTCC‑GGCCTCTGCAACAATAAATACTTTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGa                  <  1:1141315/150‑1 (MQ=255)

                                                                                                                                      gcataggCCGGAATCATCCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACg                >  2:2018547/6‑150 (MQ=255)

                                                                                                                                       cgcTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGt               >  1:713513/1‑150 (MQ=255)

                                                                                                                                        ataggCCGGAATCATCCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCCTGGCGTCGACATTTTGGTGTGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTCACGtt              >  1:2222923/4‑150 (MQ=255)

                                                                                                                                         taggCCGGAATCATCCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTtc             >  2:1953637/3‑150 (MQ=255)

                                                                                                                                         cTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTtc             <  2:1204856/150‑1 (MQ=255)

                                                                                                                                          tGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTtct            <  2:27078/150‑1 (MQ=255)

                                                                                                                                          aggCCGGAATCATCCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTtct            >  1:1677236/2‑150 (MQ=255)

                                                                                                                                            gCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGAt                                >  2:321819/1‑128 (MQ=255)

                                                                                                                                            gCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGAt                                <  1:321819/128‑1 (MQ=255)

                                                                                                                                            gCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCg          <  1:214215/150‑1 (MQ=255)

                                                                                                                                            gCCGGAATCATCCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCg          <  2:1677236/150‑1 (MQ=255)

                                                                                                                                              cGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGt        >  2:1333661/1‑150 (MQ=255)

                                                                                                                                               gggATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTc       <  2:869285/150‑1 (MQ=255)

                                                                                                                                               gggATGATTCC‑GGCCTATGCGACAGCAGGTGAACTGATATAGGTTGTGAAGAACTGGATGCGCAGTATGGCGCACATCTATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTc       <  1:2342443/150‑1 (MQ=255)

                                                                                                                                                    gatTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGa  <  2:205396/150‑1 (MQ=255)

TGTTCGGCGCGGTGGTCGGTGTACCGCCGGTGACCAGCTATGTGGAGAGTGCTGCGGGGGTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCC‑GGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGA  >  NC_002947/4875786‑4876083

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