BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,572,526	T→C	D234G (GAC→GGC)	pqqL ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,572,526	0	T	C	98.7%	476.9 / ‑4.0	149	D234G (GAC→GGC)	pqqL	putative peptidase
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (72/75); total (73/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

GCTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGTCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCG  >  NC_000913/1572436‑1572615

gCTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCGTTTATCAGCGCCAGCGCTTCTTTACTGCCg                                                                                         >  2:4636992/1‑93 (MQ=255)

gCTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCg                                                                                         >  1:1868026/1‑93 (MQ=255)

  ttttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGat                                                                                       >  2:2910016/1‑93 (MQ=255)

   tttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGGTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGata                                                                                      <  2:2319882/93‑1 (MQ=255)

   tttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGata                                                                                      <  2:2894815/93‑1 (MQ=255)

   tttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGata                                                                                      <  1:2723106/93‑1 (MQ=255)

    ttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGatat                                                                                     <  1:3971828/93‑1 (MQ=255)

         ggCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCg                                                                                <  1:3992423/93‑1 (MQ=255)

          gCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGa                                                                               >  1:2631951/1‑93 (MQ=255)

           ccAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGAc                                                                              <  2:2313490/93‑1 (MQ=255)

           ccAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGAc                                                                              <  1:1364037/93‑1 (MQ=255)

             aGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCa                                                                            <  2:1031518/93‑1 (MQ=255)

              gACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCAc                                                                           >  1:544804/1‑93 (MQ=255)

              gACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCAc                                                                           >  1:6715/1‑93 (MQ=255)

                cgcgATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGa                                                                         <  2:480959/93‑1 (MQ=255)

                cgcgATTTTCAGCTGCTTTGTTAGCCGGAACCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGa                                                                         <  1:604012/93‑1 (MQ=255)

                 gcgATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGAt                                                                        >  1:2838082/1‑93 (MQ=255)

                  cgATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTCCTGTCGATATCGCCGCCCACGATa                                                                       <  2:656147/93‑1 (MQ=255)

                   gATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATaa                                                                      >  2:4612748/1‑93 (MQ=255)

                   gATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATaa                                                                      >  1:4693269/1‑93 (MQ=255)

                    aTTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATaaa                                                                     <  2:161235/93‑1 (MQ=255)

                     ttttCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAg                                                                    <  1:105124/93‑1 (MQ=255)

                     ttttCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAg                                                                    >  2:585700/1‑93 (MQ=255)

                      tttCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAgg                                                                   <  1:1875453/93‑1 (MQ=255)

                      tttCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAgg                                                                   <  2:2555194/93‑1 (MQ=255)

                       ttCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGt                                                                  >  1:4093447/1‑93 (MQ=255)

                       ttCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGt                                                                  <  1:3607195/93‑1 (MQ=255)

                        tCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTc                                                                 <  2:2592899/93‑1 (MQ=255)

                         cAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCa                                                                <  2:3352415/93‑1 (MQ=255)

                         cAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCa                                                                >  2:656060/1‑93 (MQ=255)

                           gctgctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCata                                                              >  2:1228829/1‑93 (MQ=255)

                            ctgctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCatat                                                             <  2:2529988/93‑1 (MQ=255)

                             tgctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCAtatt                                                            >  2:3676369/1‑93 (MQ=255)

                                tttGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATattatt                                                         <  2:3012990/93‑1 (MQ=255)

                                tttGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATattatt                                                         <  1:982477/93‑1 (MQ=255)

                                    ttAGCCGGAATCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGt                                                     <  1:3359284/93‑1 (MQ=255)

                                    ttAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGt                                                     <  2:4536623/93‑1 (MQ=255)

                                    ttAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGt                                                     <  1:79387/93‑1 (MQ=255)

                                     tAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGtt                                                    <  2:1869423/93‑1 (MQ=255)

                                      aGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTg                                                   >  2:653500/1‑93 (MQ=255)

                                      aGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTg                                                   >  2:653501/1‑93 (MQ=255)

                                        ccGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGAt                                                 <  2:1612014/93‑1 (MQ=255)

                                        ccGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGAt                                                 >  2:2280044/1‑93 (MQ=255)

                                         cGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATa                                                >  1:1063517/1‑93 (MQ=255)

                                         cGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATa                                                >  2:3849200/1‑93 (MQ=255)

                                          ggAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATAc                                               >  2:1210013/1‑93 (MQ=255)

                                          ggAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATAc                                               <  1:4660/93‑1 (MQ=255)

                                          ggAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTAATAc                                               >  2:2668242/1‑93 (MQ=255)

                                           gAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATAcc                                              >  1:1822569/1‑93 (MQ=255)

                                           gAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATAcc                                              >  2:2917056/1‑93 (MQ=255)

                                             aGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAg                                            <  1:4089926/93‑1 (MQ=255)

                                             aGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAg                                            >  2:3511973/1‑93 (MQ=255)

                                             aGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAg                                            >  2:3511974/1‑93 (MQ=255)

                                              gCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGc                                           <  2:4673415/93‑1 (MQ=255)

                                                ttacttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCAGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGt                                         >  1:384979/1‑93 (MQ=255)

                                                 tacttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATATCAGCGtt                                        <  2:1933889/93‑1 (MQ=255)

                                                    ttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGAt                                     >  2:3026500/1‑93 (MQ=255)

                                                    ttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGAt                                     >  1:2481727/1‑93 (MQ=255)

                                                    ttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGAt                                     >  1:795292/1‑93 (MQ=255)

                                                     taAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATa                                    <  2:3324333/93‑1 (MQ=255)

                                                        aTTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATaaaa                                 <  2:706466/93‑1 (MQ=255)

                                                          tATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAAtt                               <  2:895463/93‑1 (MQ=255)

                                                          tATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAAtt                               >  2:2417338/1‑93 (MQ=255)

                                                          tATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAAtt                               <  1:1296008/93‑1 (MQ=255)

                                                          tATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAAtt                               <  1:3676049/93‑1 (MQ=255)

                                                          tATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAAtt                               <  1:3729456/93‑1 (MQ=255)

                                                          tATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAAtt                               <  2:1440647/93‑1 (MQ=255)

                                                           aTCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTg                              <  2:4022927/93‑1 (MQ=255)

                                                           aTCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTg                              <  2:289738/93‑1 (MQ=255)

                                                           aTCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTg                              <  2:3299646/93‑1 (MQ=255)

                                                             ccTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgc                            <  1:1721811/93‑1 (MQ=255)

                                                             ccTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgc                            >  2:1974914/1‑93 (MQ=255)

                                                              cTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgcg                           >  1:4285703/1‑93 (MQ=255)

                                                                ttATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCa                         >  1:3318105/1‑93 (MQ=255)

                                                                    tagcgccagctcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATtg                     <  1:4712526/92‑1 (MQ=255)

                                                                     agcgccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATtgt                    <  1:1562298/93‑1 (MQ=255)

                                                                     agcgccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATtgt                    >  2:740743/1‑93 (MQ=255)

                                                                     agcgccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATtgt                    >  2:1454117/1‑93 (MQ=255)

                                                                      gcgccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATtgtg                   >  1:1107127/1‑93 (MQ=255)

                                                                        gccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGcc                 <  1:2003192/93‑1 (MQ=255)

                                                                         ccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCg                <  2:2540420/93‑1 (MQ=255)

                                                                           agcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgt              <  1:2910668/93‑1 (MQ=255)

                                                                            gcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgtg             >  2:3604096/1‑93 (MQ=255)

                                                                            gcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgtg             >  2:2116584/1‑93 (MQ=255)

                                                                             cgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgtgt            >  1:3795396/1‑93 (MQ=255)

                                                                             cgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgtgt            >  2:1780219/1‑93 (MQ=255)

                                                                              gcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGtgtg           >  1:4558671/1‑93 (MQ=255)

                                                                               cttcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGa          >  1:3802602/1‑93 (MQ=255)

                                                                               cttcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGa          >  2:4277900/1‑93 (MQ=255)

                                                                               cttcttTACTGCCGATATCGCCGACCACGATAAAGGCCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGa          >  2:4475428/1‑93 (MQ=255)

                                                                                ttcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGAc         <  1:575161/93‑1 (MQ=255)

                                                                                 tcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGAcc        >  1:1095135/1‑93 (MQ=255)

                                                                                 tcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGAcc        <  1:1964189/93‑1 (MQ=255)

                                                                                 tcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGAcc        >  1:2226570/1‑93 (MQ=255)

                                                                                   tttACTGCCGATATCGCCGACCACGCTAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCgg      >  2:1468921/1‑92 (MQ=255)

                                                                                   tttACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCgg      >  2:4356782/1‑92 (MQ=255)

                                                                                   tttACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGt      <  1:1541782/93‑1 (MQ=255)

                                                                                     tACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTTGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTgg    <  1:2229411/93‑1 (MQ=255)

                                                                                      aCTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGc   <  2:3284879/93‑1 (MQ=255)

                                                                                       cTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCg  <  2:2379370/93‑1 (MQ=255)

GCTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGTCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCG  >  NC_000913/1572436‑1572615

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