BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	286,834	T→C	S16P (TCC→CCC)	yagH →	CP4‑6 prophage; putative xylosidase/arabinosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	286,834	0	T	C	100.0%	1243.1 / NA	304	S16P (TCC→CCC)	yagH	CP4‑6 prophage; putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (200/104); total (200/104)

Only 100 of 309 total aligned reads displayed.

GGCGAGGCCGTTCCCGACGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTG  >  NC_000913/286736‑286931

   gAGGCCGTTCCCGACGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCtg                                                                                              >  2:558993/1‑101 (MQ=255)

     ggCCGTTCCCGACGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCtgtg                                                                                            >  1:558998/1‑101 (MQ=255)

     ggCCGTTCCCGACGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCtgtg                                                                                            >  2:558997/1‑101 (MQ=255)

       ccGTTCCCGACGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTgcc                                                                                          >  2:559003/1‑101 (MQ=255)

         gTTCCCGACGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTgccgc                                                                                        >  1:559007/1‑101 (MQ=255)

           tCCCGACGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCa                                                                                      >  2:559010/1‑101 (MQ=255)

              cGACGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGcccc                                                                                                >  1:559014/1‑88 (MQ=255)

              cGACGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAggg                                                                                   >  2:559013/1‑101 (MQ=255)

                aCGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCg                                                                                 >  1:559017/1‑101 (MQ=255)

                 cgccgcGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGa                                                                                >  2:559021/1‑101 (MQ=255)

                   ccgcGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAgg                                                                              >  1:559022/1‑101 (MQ=255)

                    cgcgACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGa                                                                             >  1:559027/1‑101 (MQ=255)

                     gcgACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGac                                                                            >  2:559029/1‑101 (MQ=255)

                       gACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGacta                                                                          >  1:559033/1‑101 (MQ=255)

                         cAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGactact                                                                        >  2:559037/1‑101 (MQ=255)

                          aGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGActacta                                                                       >  1:559038/1‑101 (MQ=255)

                           gCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCACCAGGGCGAGGACtactac                                                                      >  2:559041/1‑101 (MQ=255)

                               cATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCg                                                                  >  2:559050/1‑101 (MQ=255)

                               cATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCg                                                                  >  1:559049/1‑101 (MQ=255)

                                aTCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGc                                                                 >  1:559053/1‑101 (MQ=255)

                                 tCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGcc                                                                <  1:558605/101‑1 (MQ=255)

                                  cccATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCa                                                               >  2:559058/1‑101 (MQ=255)

                                   ccATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCAc                                                              >  1:559059/1‑101 (MQ=255)

                                     aTTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCt                                                            >  1:559065/1‑101 (MQ=255)

                                     aTTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCt                                                            >  2:559064/1‑101 (MQ=255)

                                       tAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCg                                                          >  2:559070/1‑101 (MQ=255)

                                       tAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCg                                                          >  1:559071/1‑101 (MQ=255)

                                          ccGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGAcc                                                       <  1:558561/101‑1 (MQ=255)

                                          ccGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGAcc                                                       >  2:559076/1‑101 (MQ=255)

                                          ccGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGAcc                                                       >  1:559077/1‑101 (MQ=255)

                                           cGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCt                                                      >  2:559082/1‑101 (MQ=255)

                                           cGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCt                                                      >  1:559084/1‑101 (MQ=255)

                                            ggAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCtt                                                     >  2:559086/1‑101 (MQ=255)

                                              aGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCg                                                   >  2:559090/1‑101 (MQ=255)

                                                gCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAg                                                 >  1:559093/1‑101 (MQ=255)

                                                gCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAg                                                 <  1:558467/101‑1 (MQ=255)

                                                 cAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGt                                                >  1:559096/1‑101 (MQ=255)

                                                  aaTATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTg                                               >  2:559099/1‑101 (MQ=255)

                                                   atatGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTgg                                              >  2:559103/1‑101 (MQ=255)

                                                    tatGGAAATCACTACCCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGt                                             <  1:558780/101‑1 (MQ=255)

                                                     attgAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGtt                                            >  1:559111/4‑101 (MQ=255)

                                                     atGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGtt                                            >  1:559106/1‑101 (MQ=255)

                                                     atGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGtt                                            >  2:559110/1‑101 (MQ=255)

                                                         aaaTCACTACCCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTcctg                                        <  1:558669/101‑3 (MQ=255)

                                                         aaaTCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCg                                        >  1:559116/1‑101 (MQ=255)

                                                            tCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGc                                     >  2:559121/1‑101 (MQ=255)

                                                            tCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGc                                     >  1:559123/1‑101 (MQ=255)

                                                             cACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCg                                    >  2:559126/1‑101 (MQ=255)

                                                             cACTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCg                                    >  1:559128/1‑101 (MQ=255)

                                                              aCTAACCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGt                                   >  2:559131/1‑101 (MQ=255)

                                                                 aaCCCGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTgcg                                >  1:559136/1‑101 (MQ=255)

                                                                   cccGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCa                              >  1:559139/1‑101 (MQ=255)

                                                                   cccGATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCa                              >  2:559141/1‑101 (MQ=255)

                                                                      gATACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCt                           >  1:559144/1‑101 (MQ=255)

                                                                        tACTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTAc                         >  2:559148/1‑101 (MQ=255)

                                                                         aCTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTAcc                        >  1:559149/1‑101 (MQ=255)

                                                                         aCTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTAcc                        >  1:559153/1‑101 (MQ=255)

                                                                          cTCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTTGAGTGGTTCCCGGGCGTGCGCATCTACCa                       <  1:558958/101‑1 (MQ=255)

                                                                           tCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGa                                                                                >  1:559158/1‑43 (MQ=255)

                                                                           tCACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCAc                      >  2:559156/1‑101 (MQ=255)

                                                                            cACCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGa                                                                                <  1:558546/42‑1 (MQ=255)

                                                                             aCCGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCg                                    >  2:559163/1‑85 (MQ=255)

                                                                              ccGGCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTcc                   >  2:559166/1‑101 (MQ=255)

                                                                                ggCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCg                 >  1:559170/1‑101 (MQ=255)

                                                                                 gCTTCACCCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGt                <  1:558677/101‑1 (MQ=255)

                                                                                 gCTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGt                >  1:559173/1‑101 (MQ=255)

                                                                                  cTTCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTg               >  2:559176/1‑101 (MQ=255)

                                                                                   ttCAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGa              >  1:559177/1‑101 (MQ=255)

                                                                                     caccccGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCggg                                      <  2:558634/72‑1 (MQ=255)

                                                                                     cAACCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGact                                                                           <  2:558527/38‑1 (MQ=38)

                                                                                      aaCCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGact                                                                           <  2:558651/37‑1 (MQ=38)

                                                                                      aaCCCGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCt           >  1:559180/1‑101 (MQ=255)

                                                                                        cccGGCCCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTg               <  2:558574/95‑1 (MQ=255)

                                                                                        cccGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGAc             <  1:558319/97‑1 (MQ=255)

                                                                                        cccGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGa         >  2:559183/1‑101 (MQ=255)

                                                                                         ccGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGaa        >  2:559189/1‑101 (MQ=255)

                                                                                         ccGGACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGaa        >  1:559186/1‑101 (MQ=255)

                                                                                           ggccccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCa                       <  1:558579/81‑1 (MQ=255)

                                                                                           ggACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGt                                   <  2:558465/72‑1 (MQ=255)

                                                                                           ggACCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTg                                  <  2:558697/73‑1 (MQ=255)

                                                                                             aCCCGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAAc    >  1:559192/1‑101 (MQ=255)

                                                                                              cccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTgg                                              <  2:558482/58‑1 (MQ=255)

                                                                                              cccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGtt                                            <  2:558554/60‑1 (MQ=255)

                                                                                              cccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCgg                                       <  1:558722/65‑1 (MQ=255)

                                                                                              cccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGt                <  2:558673/88‑1 (MQ=255)

                                                                                              cccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCt           <  2:558495/93‑1 (MQ=255)

                                                                                              cccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGaaaaa     <  2:558773/99‑1 (MQ=255)

                                                                                              cccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGaaaa      <  2:558625/98‑1 (MQ=255)

                                                                                              cccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAAc    <  2:558843/100‑1 (MQ=255)

                                                                                              cccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACt   <  2:558714/101‑1 (MQ=255)

                                                                                               ccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATc                                                                   <  1:558556/36‑1 (MQ=255)

                                                                                               ccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCt                                                            <  2:558367/43‑1 (MQ=255)

                                                                                               ccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTcc                                          <  1:558515/61‑1 (MQ=255)

                                                                                               ccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTTCCGGGCGTgcgc                               <  2:558400/72‑1 (MQ=255)

                                                                                               ccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGc                                     <  1:558626/66‑1 (MQ=255)

                                                                                               ccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTg                                  <  2:558580/69‑1 (MQ=255)

                                                                                               ccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTAc                         <  1:558391/78‑1 (MQ=255)

                                                                                               ccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACt                     <  2:558612/82‑1 (MQ=255)

                                                                                               ccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGaaaaa     <  1:558712/98‑1 (MQ=255)

                                                                                               ccGCCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTg  <  1:558735/101‑1 (MQ=255)

GGCGAGGCCGTTCCCGACGCCGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTG  >  NC_000913/286736‑286931

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