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	99.6%	735.3 / NA	226	D234G (GAC→GGC)	pqqL	putative peptidase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (123/102); total (124/102)
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 230 total aligned reads displayed.

GCTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGTCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCG  >  NC_000913/1572436‑1572615

gCTTTTGTTGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCg                                                                                         >  1:1065577/1‑93 (MQ=255)

 cTTTTGTTGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGa                                                                                        <  2:2425362/93‑1 (MQ=255)

 cTTTTGTTGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGa                                                                                        <  2:4996444/93‑1 (MQ=255)

  ttttGTTGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGat                                                                                       >  1:5669410/1‑93 (MQ=255)

  ttttGTTGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGat                                                                                       >  1:3460350/1‑93 (MQ=255)

  ttttGTTGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGat                                                                                       >  2:1142279/1‑93 (MQ=255)

  ttttGTTGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGat                                                                                       >  1:2976289/1‑93 (MQ=255)

    ttgttgGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGatat                                                                                     >  1:2398829/1‑93 (MQ=255)

       ttggCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGc                                                                                  <  2:2309443/91‑1 (MQ=255)

         ggCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCg                                                                                >  1:4797724/1‑93 (MQ=255)

         ggCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCg                                                                                >  2:1646352/1‑93 (MQ=255)

           ccAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGAc                                                                              >  1:6254891/1‑93 (MQ=255)

              gACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCAc                                                                           <  2:37687/93‑1 (MQ=255)

               aCGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACg                                                                          >  1:5744248/1‑93 (MQ=255)

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

                cgcgATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGa                                                                         >  2:2265260/1‑93 (MQ=255)

                 ggggttttCAGCTACTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGAt                                                                        >  2:6063715/5‑93 (MQ=255)

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

                  cgATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATa                                                                       >  2:2648463/1‑93 (MQ=255)

                  cgATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATa                                                                       <  1:1938658/93‑1 (MQ=255)

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

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

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

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

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

                         cAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCGAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCa                                                                <  2:5305337/93‑1 (MQ=255)

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

                           gctgctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCata                                                              <  2:846768/93‑1 (MQ=255)

                            ctgctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCatat                                                             >  2:2090797/1‑93 (MQ=255)

                             tgctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGCTGAAGGTCAtatt                                                            >  1:2008105/1‑93 (MQ=255)

                             tgctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCAtatt                                                            <  2:4586735/93‑1 (MQ=255)

                             tgctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCAtatt                                                            <  1:4625861/93‑1 (MQ=255)

                             tgctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCAtatt                                                            <  2:2221341/93‑1 (MQ=255)

                              gctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCAtatta                                                           >  1:5060413/1‑93 (MQ=255)

                              gctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCAtatta                                                           >  2:1399241/1‑93 (MQ=255)

                                tttGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATattatt                                                         >  1:5830033/1‑93 (MQ=255)

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

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

                                 ttGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTAttt                                                        <  2:3595707/93‑1 (MQ=255)

                                   gTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTgg                                                      <  2:254042/93‑1 (MQ=255)

                                       gCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGa                                                  >  1:4005101/1‑93 (MQ=255)

                                       gCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGa                                                  <  2:2139620/93‑1 (MQ=255)

                                       gCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGa                                                  <  1:3837903/93‑1 (MQ=255)

                                         cGGAAGCTTACTTAAATTATCCTTTATAAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATa                                                >  1:1876129/1‑93 (MQ=255)

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

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

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

                                           gAAGCTTACTTAAATTATCCTTTATAAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATAcc                                              >  1:4555762/1‑93 (MQ=255)

                                            aaGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCa                                             <  2:485225/93‑1 (MQ=255)

                                            aaGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCa                                             >  2:2998215/1‑93 (MQ=255)

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

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

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

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

                                                ttacttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGt                                         >  1:6118635/1‑93 (MQ=255)

                                                ttacttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGt                                         >  1:597948/1‑93 (MQ=255)

                                                ttacttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGt                                         <  1:2805003/93‑1 (MQ=255)

                                                 tacttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGtt                                        <  1:3927835/93‑1 (MQ=255)

                                                  tcttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGCCCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTg                                       <  2:3932170/92‑1 (MQ=255)

                                                   cttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGa                                      >  1:4278022/1‑93 (MQ=255)

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

                                                    ttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGAt                                     <  1:1365774/93‑1 (MQ=255)

                                                       aaTTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATaaa                                  >  2:2495350/1‑93 (MQ=255)

                                                       aaTTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATaaa                                  <  2:5608569/93‑1 (MQ=255)

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

                                                        aTTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATaaaa                                 >  1:5302538/1‑93 (MQ=255)

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

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

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

                                                            tCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTgg                             >  1:3345873/1‑93 (MQ=255)

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

                                                               tttATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgcgc                          >  1:4439441/1‑93 (MQ=255)

                                                               tttATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgcgc                          <  1:1000016/93‑1 (MQ=255)

                                                                    cagcgccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATtg                     <  1:862626/93‑1 (MQ=255)

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

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

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

                                                                         ccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGcct                <  2:1866831/93‑2 (MQ=255)

                                                                          cagcgcTTCTTTACTGCCGATATAGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCgg               <  1:5376467/93‑1 (MQ=255)

                                                                           agcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCggg              >  2:1096930/1‑92 (MQ=255)

                                                                           agcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgt              >  2:2498272/1‑93 (MQ=255)

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

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

                                                                             cgcTTCTCTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgtgt            <  2:3331854/93‑1 (MQ=255)

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

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

                                                                                ttcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGAc         >  1:2572239/1‑93 (MQ=255)

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

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

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

                                                                                 tcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTAGCGCAATTGTGCCGGTGTGAcc        >  1:2481421/1‑93 (MQ=255)

                                                                                  cttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCg       >  1:4712837/1‑93 (MQ=255)

                                                                                    ttACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTg     >  1:2304261/1‑93 (MQ=255)

                                                                                    ttACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTg     <  1:1462483/93‑1 (MQ=255)

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

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

                                                                                      aCTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGc   >  1:6231845/1‑93 (MQ=255)

                                                                                       cTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCg  >  1:2346916/1‑93 (MQ=255)

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

                                                                                       cTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCg  >  2:5785262/1‑93 (MQ=255)

GCTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGTCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCG  >  NC_000913/1572436‑1572615

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