BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,474,826	+A	intergenic (+1/+44)	yjgL → / ← argI	hypothetical protein/ornithine carbamoyltransferase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,474,823	1	.	A	100.0%	741.1 / NA	193	*617K (TAA→AAA)	yjgL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base A (135/58); total (135/58)

Only 100 of 208 total aligned reads displayed.

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910

ttAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat                                                                                          >  1:3263506/1‑93 (MQ=255)

 tAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑a                                                                                        >  1:1490250/1‑93 (MQ=255)

   aTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaa                                                                                      <  1:3556224/93‑1 (MQ=255)

   aTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaa                                                                                      <  1:2394921/93‑1 (MQ=255)

   aTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaa                                                                                      >  2:4241987/1‑93 (MQ=255)

    tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaa                                                                                     >  2:4341900/1‑92 (MQ=255)

          gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                >  1:4737541/1‑93 (MQ=255)

           cGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgtg                                                                               >  1:4520817/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  1:2136849/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  1:5320851/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             <  2:6135925/93‑1 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  2:2957303/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  2:5784015/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  1:2606263/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  1:2207111/1‑93 (MQ=255)

              ttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATg                                                                            >  2:4072774/1‑93 (MQ=255)

               taaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGc                                                                           >  1:1444263/1‑93 (MQ=255)

                 aataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCg                                                                         >  1:3860900/1‑93 (MQ=255)

                   taaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCggg                                                                       >  2:4418337/1‑93 (MQ=255)

                    aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGa                                                                      <  2:5756829/93‑1 (MQ=255)

                     aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGAt                                                                     >  1:201078/1‑93 (MQ=255)

                     aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGAt                                                                     >  1:3993157/1‑93 (MQ=255)

                        tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGt                                                                  >  1:7079966/1‑93 (MQ=255)

                        tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGt                                                                  >  2:515769/1‑93 (MQ=255)

                            tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgt                                                              >  2:2147035‑M1/1‑91 (MQ=255)

                             ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgta                                                             <  2:7314603‑M1/93‑4 (MQ=255)

                               aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatt                                                           >  2:248046‑M1/1‑88 (MQ=255)

                                ccGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtattt                                                          >  2:1739715‑M1/1‑87 (MQ=255)

                                 cGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttc                                                         >  1:7132869‑M1/1‑86 (MQ=255)

                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  2:6737789‑M1/1‑85 (MQ=255)

                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        <  2:6073090‑M1/93‑9 (MQ=255)

                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  1:3358673‑M1/1‑85 (MQ=255)

                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  1:3587289‑M1/1‑85 (MQ=255)

                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  1:2348674‑M1/1‑85 (MQ=255)

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttctc                                                       >  1:4567417‑M1/1‑84 (MQ=255)

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                       >  1:5142988‑M1/1‑84 (MQ=255)

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                       >  2:3653704‑M1/1‑84 (MQ=255)

                                     aaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccgg                                                     <  2:5585901‑M1/93‑12 (MQ=255)

                                     aaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccgg                                                     >  1:751716‑M1/1‑82 (MQ=255)

                                     aaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccgg                                                     <  2:7098042‑M1/93‑12 (MQ=255)

                                        tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                  >  2:1922416‑M1/1‑79 (MQ=255)

                                        tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                  >  2:3770524‑M1/1‑79 (MQ=255)

                                        tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                  >  2:7431505‑M1/1‑79 (MQ=255)

                                        tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTACATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                  >  1:7301603‑M1/1‑79 (MQ=255)

                                            cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                              >  2:2815978‑M1/1‑75 (MQ=255)

                                            cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                              >  1:1705338‑M1/1‑75 (MQ=255)

                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             >  1:4589434‑M1/1‑74 (MQ=255)

                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           <  1:997705‑M1/93‑22 (MQ=255)

                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           >  2:7392502‑M1/1‑72 (MQ=255)

                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           >  2:5814744‑M1/1‑72 (MQ=255)

                                                taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatag                                          >  1:5529158‑M1/1‑71 (MQ=255)

                                                 aGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                                         <  1:488848‑M1/93‑24 (MQ=255)

                                                  gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                        >  2:4753686‑M1/1‑69 (MQ=255)

                                                  gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                        >  2:3361461‑M1/1‑69 (MQ=255)

                                                   cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                                       <  1:6419971‑M1/93‑26 (MQ=255)

                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     <  1:3005365‑M1/93‑28 (MQ=255)

                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  1:6835085‑M1/1‑66 (MQ=255)

                                                      tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                    >  2:5515771‑M1/1‑65 (MQ=255)

                                                         aaCACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                 >  1:3192471‑M1/1‑62 (MQ=255)

                                                          acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                >  2:1160931‑M1/1‑61 (MQ=255)

                                                           cacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatt                               >  1:6632477‑M1/1‑60 (MQ=255)

                                                             cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             <  1:5331345‑M1/93‑36 (MQ=255)

                                                             cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  2:1029154‑M1/1‑58 (MQ=255)

                                                             cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  1:5362104‑M1/1‑58 (MQ=255)

                                                             cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  1:5230342‑M1/1‑58 (MQ=255)

                                                               ttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattt                           <  1:1000066‑M1/93‑38 (MQ=255)

                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTCgtatttcccggcatctttatagcgatagcaattatttac                         >  1:1947927‑M1/1‑54 (MQ=255)

                                                                   gTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactg                       <  2:2284152‑M1/93‑42 (MQ=255)

                                                                   gTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactg                       <  1:6136106‑M1/93‑42 (MQ=255)

                                                                   gTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactg                       >  2:6167240‑M1/1‑52 (MQ=255)

                                                                   gTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactg                       <  2:445263‑M1/93‑42 (MQ=255)

                                                                    tGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                      >  2:6554367‑M1/1‑51 (MQ=255)

                                                                     ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     >  1:416350‑M1/1‑50 (MQ=255)

                                                                     ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     <  1:3482403‑M1/93‑44 (MQ=255)

                                                                         gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtc                 <  1:131265‑M1/93‑48 (MQ=255)

                                                                         gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtc                 >  2:3332746‑M1/1‑46 (MQ=255)

                                                                           aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc               >  2:3985611‑M1/1‑44 (MQ=255)

                                                                           aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc               >  2:3206957‑M1/1‑44 (MQ=255)

                                                                           aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc               >  2:5401494‑M1/1‑44 (MQ=255)

                                                                           aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc               >  1:3615920‑M1/1‑44 (MQ=255)

                                                                           aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc               >  2:5781443‑M1/1‑44 (MQ=255)

                                                                           aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc               >  1:6579290‑M1/1‑44 (MQ=255)

                                                                            tAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcg              >  1:2871891‑M1/1‑43 (MQ=255)

                                                                            tAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcg              <  1:4409347‑M1/93‑51 (MQ=255)

                                                                              gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgac            >  1:5449071‑M1/1‑41 (MQ=255)

                                                                              gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgac            >  1:5598764‑M1/1‑41 (MQ=255)

                                                                               aCCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacc           >  2:4636259‑M1/1‑40 (MQ=255)

                                                                               aCCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacc           <  1:5286509‑M1/93‑54 (MQ=255)

                                                                                ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacca          >  1:1384933‑M1/1‑39 (MQ=255)

                                                                                ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacca          <  2:1608475‑M1/93‑55 (MQ=255)

                                                                                 cTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccat         >  2:261714‑M1/1‑38 (MQ=255)

                                                                                 cTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccat         <  2:452460‑M1/93‑56 (MQ=255)

                                                                                   aTCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatca       <  1:3107051‑M1/93‑58 (MQ=255)

                                                                                    tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcac      <  1:7355399‑M1/93‑59 (MQ=255)

                                                                                     cACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcacc     <  1:7196633‑M1/93‑60 (MQ=255)

                                                                                      aCTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccg    >  2:7374578‑M1/1‑33 (MQ=255)

                                                                                       cTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgc   >  1:2286075‑M1/1‑32 (MQ=255)

                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgct  <  2:916789‑M1/93‑63 (MQ=255)

                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgct  >  2:7496255‑M1/1‑31 (MQ=255)

                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgct  >  2:3901841‑M1/1‑31 (MQ=255)

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910

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