BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,862,881	T→C	L208P (CTA→CCA)	ygbM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,862,881	0	T	C	99.8%	1304.1 / ‑5.2	405	L208P (CTA→CCA)	ygbM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (172/232); total (172/233)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

Only 100 of 430 total aligned reads displayed.

ATACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACTACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCG  >  NC_000913/2862791‑2862970

aTACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCAc                                                                                         <  1:449139/93‑1 (MQ=255)

 tACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGaccacc                                                                                        <  2:9868531/93‑1 (MQ=255)

 tACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGaccacc                                                                                        >  2:2439476/1‑93 (MQ=255)

  aCTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCt                                                                                       >  2:8516063/1‑93 (MQ=255)

  aCTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCt                                                                                       <  1:4791383/93‑1 (MQ=255)

   cTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTg                                                                                      <  2:5248727/93‑1 (MQ=255)

    ctttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTGAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGa                                                                                     <  1:4156703/92‑1 (MQ=255)

     tttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACTAATTGCCGGACCACCTGAt                                                                                    >  1:8933057/1‑93 (MQ=255)

     tttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGCCCACCTGAt                                                                                    >  1:254875/1‑93 (MQ=255)

     tttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACTACCTGAt                                                                                    <  1:1487026/93‑1 (MQ=255)

     tttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGAt                                                                                    >  2:2052424/1‑93 (MQ=255)

          tGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACa                                                                               <  2:4250909/93‑1 (MQ=255)

          tGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACa                                                                               >  2:3438315/1‑93 (MQ=255)

          tGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACa                                                                               <  1:2143270/93‑1 (MQ=255)

          tGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACa                                                                               <  2:3253517/93‑1 (MQ=255)

           gcgcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACAt                                                                              <  2:7822644/93‑1 (MQ=255)

           gcgcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACAt                                                                              <  2:8296063/93‑1 (MQ=255)

           gcgcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACAt                                                                              <  1:1040982/93‑1 (MQ=255)

              cAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGaa                                                                           <  2:2415917/93‑1 (MQ=255)

               aaaaaCTAGATGGCAACCTGAAGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAAc                                                                          <  1:2784305/93‑1 (MQ=255)

                aaaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAAcc                                                                         >  1:3780895/1‑93 (MQ=255)

                 aaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCg                                                                        <  2:485551/93‑1 (MQ=255)

                  aaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCgg                                                                       >  2:5571449/1‑93 (MQ=255)

                  aaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCgg                                                                       <  1:5390713/93‑1 (MQ=255)

                  aaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCgg                                                                       >  1:7587191/1‑93 (MQ=255)

                   aGTAGATGGCAACCTGACGCATTTAATCCGCGACTAAGCCGGGAAATATGCGCACGTACAAATTGCCGAACCACCTGATCGACATGAACCGGa                                                                      >  1:6896080/1‑93 (MQ=255)

                    gTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGAc                                                                     <  1:4193683/93‑1 (MQ=255)

                    gTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGAc                                                                     >  1:7276689/1‑93 (MQ=255)

                      aGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGa                                                                   <  1:5101804/93‑1 (MQ=255)

                      aGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGa                                                                   <  1:6737453/93‑1 (MQ=255)

                        agtgCAACCGGACGCATTGAATCCGCGTCTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATg                                                                 <  1:4955560/90‑1 (MQ=255)

                        aTGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATg                                                                 <  2:5926657/93‑1 (MQ=255)

                         tGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATgg                                                                <  2:3764196/93‑1 (MQ=255)

                         tGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATgg                                                                >  1:5959480/1‑93 (MQ=255)

                         tGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATgg                                                                >  1:7382974/1‑93 (MQ=255)

                         tGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGAACAAATTGCCGGACCACCTGATCGACATGAACCGGACCATgg                                                                >  1:7112585/1‑93 (MQ=255)

                          ggCACCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGga                                                               <  2:7974173/93‑1 (MQ=255)

                           gCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGgag                                                              >  1:9432851/1‑93 (MQ=255)

                             aaCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaa                                                            >  2:667974/1‑93 (MQ=255)

                             aaCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaa                                                            <  1:8480423/93‑1 (MQ=255)

                             aaCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaa                                                            >  2:9551089/1‑93 (MQ=255)

                              aCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaac                                                           >  1:2032426/1‑92 (MQ=255)

                              aCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaaa                                                           >  2:5104041/1‑93 (MQ=255)

                               ccTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAAt                                                          <  1:1640222/93‑1 (MQ=255)

                                cTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCTGGACCACCTGATCGACATGAACCGGACGATGGAGAAATc                                                         >  1:3988293/1‑93 (MQ=255)

                                cTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATc                                                         <  2:5975140/93‑1 (MQ=255)

                                     gCATTTACTCCGCGACTACGCCGGGCAATATGCGCACGTACAAATTGGCGGACCACCTGATCGACATGAACCGGACGATGGAGGCATCAACTa                                                    >  1:4896490/1‑93 (MQ=255)

                                     gCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTa                                                    <  1:6632090/93‑1 (MQ=255)

                                      cATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTAt                                                   >  2:7097673/1‑93 (MQ=255)

                                        tttAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATcc                                                 >  2:9699640/1‑93 (MQ=255)

                                        tttAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATcc                                                 >  2:6620940/1‑93 (MQ=255)

                                        tttAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATcc                                                 <  2:2959430/93‑1 (MQ=255)

                                         ttAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCg                                                <  2:412492/93‑1 (MQ=255)

                                           aaTCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTg                                              <  2:6954895/93‑1 (MQ=255)

                                           aaTCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTg                                              >  1:1880502/1‑93 (MQ=255)

                                           aaTCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTg                                              <  2:2779515/93‑1 (MQ=255)

                                            aTCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTgg                                             <  2:8745180/93‑1 (MQ=255)

                                             tCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGc                                            >  2:7651718/1‑93 (MQ=255)

                                             tCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGc                                            >  2:4415923/1‑93 (MQ=255)

                                             tCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGc                                            >  1:1100438/1‑93 (MQ=255)

                                             cccGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGc                                            <  2:6491631/92‑1 (MQ=255)

                                               cgcgACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTg                                          >  1:645631/1‑93 (MQ=255)

                                                gcgACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGt                                         >  2:2884012/1‑93 (MQ=255)

                                                 cgACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGtt                                        <  1:5202439/93‑1 (MQ=255)

                                                  gACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTc                                       >  1:7160690/1‑93 (MQ=255)

                                                   aCTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTcc                                      >  2:3839399/1‑93 (MQ=255)

                                                   aCTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTcc                                      >  1:8574154/1‑93 (MQ=255)

                                                    cTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCg                                     >  1:7889498/1‑93 (MQ=255)

                                                     tACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGc                                    >  2:3570082/1‑93 (MQ=255)

                                                      aCGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGcc                                   <  2:9116542/93‑1 (MQ=255)

                                                      aCGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGcc                                   <  1:4282048/93‑1 (MQ=255)

                                                       cGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCt                                  >  1:10041742/1‑93 (MQ=255)

                                                         ccGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGt                                <  2:6289290/93‑1 (MQ=255)

                                                            ggAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCg                             >  1:5537308/1‑93 (MQ=255)

                                                            ggAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCg                             >  1:5725976/1‑93 (MQ=255)

                                                             gAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGa                            <  2:4500862/93‑1 (MQ=255)

                                                              aaaTATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGAt                           >  1:1339464/1‑93 (MQ=255)

                                                               aaTATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATg                          >  2:5824688/1‑93 (MQ=255)

                                                                atatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGa                         >  1:2858109/1‑93 (MQ=255)

                                                                 tatGCGCACGTACAAATTGCCTGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAg                        <  1:2319074/93‑1 (MQ=255)

                                                                 aatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAg                        <  1:1567124/92‑1 (MQ=255)

                                                                  atGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAgg                       <  1:3215003/93‑1 (MQ=255)

                                                                   tGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGt                      >  2:1416636/1‑93 (MQ=255)

                                                                   tGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGt                      <  1:974939/93‑1 (MQ=255)

                                                                    gcgcACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTg                     >  2:1048656/1‑93 (MQ=255)

                                                                       cACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGa                  <  2:7277065/93‑1 (MQ=255)

                                                                         cGTACAAATTGCCGGACCACCTGATCGACCTGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGata                >  1:7219340/1‑93 (MQ=255)

                                                                         cGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGata                <  1:9905512/93‑1 (MQ=255)

                                                                          gTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGGTGAGGTGGGgtat               >  2:4661684/1‑93 (MQ=255)

                                                                            aCAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCa             <  1:6183665/93‑1 (MQ=255)

                                                                            aCAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATAAGGTGGGATATCa             <  1:341586/93‑1 (MQ=255)

                                                                                aTTGCCCGACCCCCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGc         <  1:8033004/93‑1 (MQ=255)

                                                                                 ttGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCt        >  2:8893377/1‑93 (MQ=255)

                                                                                  tGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTg       <  2:1131342/93‑1 (MQ=255)

                                                                                  tGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTg       <  1:3601702/93‑1 (MQ=255)

                                                                                     cGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGAt    >  2:7781835/1‑93 (MQ=255)

                                                                                     cGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGAt    <  1:877727/93‑1 (MQ=255)

                                                                                     cGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGAt    <  1:5882456/93‑1 (MQ=255)

                                                                                      ggACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATc   >  2:8163595/1‑93 (MQ=255)

                                                                                       gACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCg  <  1:9547225/93‑1 (MQ=255)

ATACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACTACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCG  >  NC_000913/2862791‑2862970

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