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%	1001.7 / NA	273	*617K (TAA→AAA)	yjgL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base A (157/116); total (157/116)

Only 100 of 302 total aligned reads displayed.

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910

ttAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat                                                                                          <  2:6901324/93‑1 (MQ=255)

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

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

     gAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaat                                                                                    >  2:3670222/1‑91 (MQ=255)

      aaCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaata                                                                                   >  2:6133382/1‑90 (MQ=255)

      aaCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaata                                                                                   <  1:1919424/93‑4 (MQ=255)

       aCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaataa                                                                                  <  1:2934822/93‑5 (MQ=255)

         aGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtg                                                                                 <  1:5875953/93‑1 (MQ=255)

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

          gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                <  2:5558217/93‑1 (MQ=255)

           cGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgtg                                                                               <  1:4469940/93‑1 (MQ=255)

            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              >  1:534105/1‑93 (MQ=255)

            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              <  2:3240586/93‑1 (MQ=255)

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

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

                aaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGcc                                                                          >  1:3295763/1‑93 (MQ=255)

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

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

                    aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGa                                                                      >  2:1934937/1‑93 (MQ=255)

                    aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGa                                                                      >  2:4888115/1‑93 (MQ=255)

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

                       cTAAATCCACCGGTAATTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACTTATCTCTGTATAAAATAAAGTGATGCCGGGATgg                                                                   >  2:7017408/1‑93 (MQ=255)

                       cTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATgg                                                                   <  2:2410567/93‑1 (MQ=255)

                       cTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATgg                                                                   >  2:800934/1‑93 (MQ=255)

                         aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGtt                                                                 <  1:416199/93‑1 (MQ=255)

                         aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGtt                                                                 >  2:6694918/1‑93 (MQ=255)

                          aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGttt                                                                >  1:221513/1‑93 (MQ=255)

                          aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGttt                                                                >  1:1074953/1‑93 (MQ=255)

                           aTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTg                                                               >  1:3881120‑M1/1‑92 (MQ=255)

                           aTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTg                                                               >  1:5263481‑M1/1‑92 (MQ=255)

                            tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGAGGAAGAATAGAAATATCACTGCATAAAATAATATGATGCCGGGATGGTTTgt                                                              >  1:3931058‑M1/1‑91 (MQ=255)

                             ccACCGGAAAATTTACGTATAGCAATAGACAACTTTGGGTGGAAGAATAGACCAATCACTGCATAAAATAATTTGATGCCGGGATGGttttgt                                                             >  1:1740789/1‑90 (MQ=255)

                              cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtat                                                            >  1:4686567‑M1/1‑89 (MQ=255)

                              cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtat                                                            <  2:3107580‑M1/93‑5 (MQ=255)

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

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

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

                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        <  1:1770566‑M1/93‑9 (MQ=255)

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

                                   gAAAATTTAAGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                       <  2:2277661‑M1/93‑10 (MQ=255)

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

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

                                         ttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatc                                                 >  2:2733004‑M1/1‑78 (MQ=255)

                                          tACGTATAGCAATAGAACACTTTGGGTGGAAGAATCGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatct                                                >  1:3450523‑M1/1‑77 (MQ=255)

                                          tACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatct                                                <  2:3817243‑M1/93‑17 (MQ=255)

                                           aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               <  2:6558638‑M1/93‑18 (MQ=255)

                                           aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  1:3524320‑M1/1‑76 (MQ=255)

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

                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             <  2:7204578‑M1/93‑20 (MQ=255)

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

                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           <  2:159351‑M1/93‑22 (MQ=255)

                                                taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatag                                          <  2:7695371‑M1/93‑23 (MQ=255)

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

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

                                                  gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                        <  1:4270361‑M1/93‑25 (MQ=255)

                                                    aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                      >  2:745760‑M1/1‑67 (MQ=255)

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

                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     <  2:7588091‑M1/93‑28 (MQ=255)

                                                        gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                  >  2:3755741‑M1/1‑63 (MQ=255)

                                                          acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                <  2:4629978‑M1/93‑33 (MQ=255)

                                                           cacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatt                               <  2:685526‑M1/93‑34 (MQ=255)

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

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

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

                                                             cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             <  2:5899024‑M1/93‑36 (MQ=255)

                                                              tttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                            <  2:5055108‑M1/93‑37 (MQ=255)

                                                              tttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                            <  1:6917466‑M1/93‑37 (MQ=255)

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

                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  2:6403983‑M1/93‑40 (MQ=255)

                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  2:3308021‑M1/93‑40 (MQ=255)

                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  1:3757808‑M1/93‑40 (MQ=255)

                                                                  ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                        <  2:2800044‑M1/93‑41 (MQ=255)

                                                                  ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                        <  1:1604930‑M1/93‑41 (MQ=255)

                                                                  ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                        >  1:7852622‑M1/1‑53 (MQ=255)

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

                                                                      gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                    >  2:4664829‑M1/1‑49 (MQ=255)

                                                                      gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                    <  2:5503268‑M1/93‑45 (MQ=255)

                                                                       aagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtg                   >  1:69481‑M1/1‑48 (MQ=255)

                                                                          aaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcg                >  1:5623220‑M1/1‑45 (MQ=255)

                                                                           aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgc               >  2:101717‑M1/1‑44 (MQ=255)

                                                                             aGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcga             >  2:6277243‑M1/1‑42 (MQ=255)

                                                                             aGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcga             <  2:182100‑M1/93‑52 (MQ=255)

                                                                              gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgac            <  1:844464‑M1/93‑53 (MQ=255)

                                                                              gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgac            <  1:3090962‑M1/93‑53 (MQ=255)

                                                                               aCCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacc           <  2:4737437‑M1/93‑54 (MQ=255)

                                                                                ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacaa          >  1:6378023‑M1/1‑39 (MQ=255)

                                                                                 cTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat         <  2:2968006‑M1/93‑56 (MQ=255)

                                                                                 cTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat         <  2:1325057‑M1/93‑56 (MQ=255)

                                                                                  tATCACTGCATAAAATAATGTGGTGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccata        >  2:4038822‑M1/1‑37 (MQ=255)

                                                                                  tATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc        >  1:1413457‑M1/1‑37 (MQ=255)

                                                                                   aTCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatca       >  2:713449‑M1/1‑36 (MQ=255)

                                                                                    tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac      <  2:394136‑M1/93‑59 (MQ=255)

                                                                                    tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac      <  1:7312307‑M1/93‑59 (MQ=255)

                                                                                     cACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcacc     <  1:1196380‑M1/93‑60 (MQ=255)

                                                                                     cACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcacc     <  2:3429288‑M1/93‑60 (MQ=255)

                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  >  1:4528440‑M1/1‑31 (MQ=255)

                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  <  1:4845978‑M1/93‑63 (MQ=255)

                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  >  2:4387091‑M1/1‑31 (MQ=255)

                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  >  2:4300630‑M1/1‑31 (MQ=255)

                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  >  2:1730409‑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` `< 25 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`