BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,074,827	A→G	Q78R (CAG→CGG)	yihZ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,074,827	0	A	G	99.7%	1217.9 / ‑5.4	380	Q78R (CAG→CGG)	yihZ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (170/209); total (171/209)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 407 total aligned reads displayed.

GTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTATATGACT  >  NC_000913/4074737‑4074916

gtgtGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGt                                                                                         <  1:3154699/93‑1 (MQ=255)

  gtgCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGttt                                                                                       <  1:4445024/93‑1 (MQ=255)

   tgCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTa                                                                                      <  1:1147647/93‑1 (MQ=255)

   tgCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTa                                                                                      >  1:7630979/1‑93 (MQ=255)

   tgCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTa                                                                                      >  1:8513229/1‑93 (MQ=255)

   tgCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTa                                                                                      <  1:9273104/93‑1 (MQ=255)

   tgCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGGTTa                                                                                      >  1:4905748/1‑93 (MQ=255)

   tgCTCGGCTACCGCATCATTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGTCGGGCGGCAGTGGACTGGTGGTTTCCCGGTTTa                                                                                      >  1:9956891/1‑93 (MQ=255)

     cTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTAcc                                                                                    <  2:7263919/93‑1 (MQ=255)

     cTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTAcc                                                                                    <  2:8186257/93‑1 (MQ=255)

      tCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGGTTCCCGGTTTAccc                                                                                   >  1:2337958/1‑93 (MQ=255)

            aCCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTcgccg                                                                             >  1:5526344/1‑93 (MQ=255)

              cGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCa                                                                           <  1:2463640/93‑1 (MQ=255)

              cGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCa                                                                           >  1:7795151/1‑93 (MQ=255)

               gCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAg                                                                          >  1:6591955/1‑93 (MQ=255)

               gCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAg                                                                          <  2:3455387/93‑1 (MQ=255)

               gCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAg                                                                          >  1:8081551/1‑93 (MQ=255)

               ccATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAg                                                                          <  2:2059222/92‑1 (MQ=255)

                  tCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATa                                                                       <  2:3384248/93‑1 (MQ=255)

                  tCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATa                                                                       <  2:9272731/93‑1 (MQ=255)

                  tCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATa                                                                       >  2:6436164/1‑93 (MQ=255)

                    tttAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATAcc                                                                     >  2:3324723/1‑93 (MQ=255)

                    tttAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATAcc                                                                     <  2:8277566/93‑1 (MQ=255)

                    tttAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATAcc                                                                     >  1:1692562/1‑93 (MQ=255)

                     ttAGCGATGCCGAAGACAAGATGAATCTCAACGAGCGACCGGCGGGCAGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATAcca                                                                    <  2:8975199/93‑2 (MQ=255)

                       aGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGaa                                                                  <  1:9906831/93‑1 (MQ=255)

                        gCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGTAc                                                                 >  1:3703306/1‑93 (MQ=255)

                        gCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAAc                                                                 >  2:2334328/1‑93 (MQ=255)

                         cGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACg                                                                <  2:2839670/93‑1 (MQ=255)

                          gATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACgg                                                               <  1:2094634/93‑1 (MQ=255)

                          gATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACgg                                                               >  1:7388108/1‑93 (MQ=255)

                          gATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGGGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACgg                                                               >  2:5394922/1‑93 (MQ=255)

                            tGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCTGGTTTACCCTCGCCGCAGATACCGAACgggg                                                             <  2:6929450/93‑1 (MQ=255)

                            tGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACgggg                                                             <  2:9354961/93‑1 (MQ=255)

                            tGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACgggg                                                             >  1:5050610/1‑93 (MQ=255)

                            tGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGGTTCCCGGTTTACCCTCGCCGCAGATACCGAACgggg                                                             >  1:9720274/1‑93 (MQ=255)

                             gCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACggggg                                                            >  1:8829814/1‑93 (MQ=255)

                              ccGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGa                                                           <  1:7195626/93‑1 (MQ=255)

                              ccGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGa                                                           <  1:6920960/93‑1 (MQ=255)

                               cGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGAt                                                          >  1:8769794/1‑93 (MQ=255)

                                gAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATg                                                         <  2:5944212/93‑1 (MQ=255)

                                 aaGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCTCCGCAGATACCGAACGGGGGATgc                                                        <  2:9651564/93‑1 (MQ=255)

                                 aaGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATgc                                                        >  2:6786296/1‑93 (MQ=255)

                                  aGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATgcg                                                       >  2:7715346/1‑93 (MQ=255)

                                     gaaGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGccc                                                    <  1:1865859/92‑1 (MQ=255)

                                      aaGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCa                                                   >  1:3570438/1‑93 (MQ=255)

                                       aGATGAATCTCAGCGTGCAAGAGGCGGTCGGCAGTGTGCTGGTGGTTGCGCGGTTTGCCCTCGCCGCAGTTACCGTACGGGTGATGCGCCCaa                                                  <  2:4188727/93‑1 (MQ=255)

                                       aGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCaa                                                  <  2:6101337/93‑1 (MQ=255)

                                       aGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCaa                                                  >  1:6679661/1‑93 (MQ=255)

                                       aGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCaa                                                  >  2:7944068/1‑93 (MQ=255)

                                        gATGAATCTCAACGTGCACGAGGCGGGTGGCAGTGTTCTGGTGGTTGCCGGGTTTACCCTCGCCGCAGATGCCGAACGGGCGATGCGCCCAAg                                                 <  2:4138070/93‑1 (MQ=255)

                                        gATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAg                                                 >  2:4006727/1‑93 (MQ=255)

                                          tGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGtt                                               >  2:277548/1‑93 (MQ=255)

                                            aaTCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTtc                                             >  1:7905644/1‑93 (MQ=255)

                                             aTCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTtct                                            <  2:3478572/93‑1 (MQ=255)

                                             aTCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTtct                                            <  1:1790988/93‑1 (MQ=255)

                                              tctcAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTtctc                                           <  2:4972555/93‑1 (MQ=255)

                                                tcAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGCTACCGAACGGGGGATGCGCCCAAGTTTCTCCa                                         >  1:5886314/1‑93 (MQ=255)

                                                tcAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCa                                         <  2:3737402/93‑1 (MQ=255)

                                                   aCGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAg                                      >  1:5986631/1‑93 (MQ=255)

                                                   aCGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAg                                      >  2:5083296/1‑93 (MQ=255)

                                                   aCGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAg                                      <  2:3640098/93‑1 (MQ=255)

                                                    cGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAgg                                     <  1:2604778/93‑1 (MQ=255)

                                                      tGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCTAGTTTCTCCAAAGGTg                                   <  2:9558563/93‑1 (MQ=255)

                                                       gCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGc                                  <  2:6711735/93‑1 (MQ=255)

                                                        cAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCa                                 <  1:3822/93‑1 (MQ=255)

                                                        cAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCa                                 <  2:1585414/93‑1 (MQ=255)

                                                            aGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCAc                             >  2:2977887/1‑93 (MQ=255)

                                                            aGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCAc                             <  2:7206363/93‑1 (MQ=255)

                                                            aGGCGGGCGGCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCCCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCAc                             <  1:3075880/93‑1 (MQ=255)

                                                                gggcggCAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGa                         <  1:1438003/93‑1 (MQ=255)

                                                                 ggcggcAGTGTGCTGGTGGTTTCCCGGTTTGCCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGAt                        <  1:411692/93‑1 (MQ=255)

                                                                 ggcggcAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGAt                        >  1:10009398/1‑93 (MQ=255)

                                                                  gcggcAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATc                       <  2:2017234/93‑1 (MQ=255)

                                                                   cggcAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATcg                      <  2:1300733/93‑1 (MQ=255)

                                                                    ggcAGTGTGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATcgc                     <  1:4316402/93‑1 (MQ=255)

                                                                        gtgtGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCag                 >  1:7975579/1‑93 (MQ=255)

                                                                        gtgtGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCag                 >  2:2186014/1‑93 (MQ=255)

                                                                        gtgtGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCag                 >  1:8333794/1‑93 (MQ=255)

                                                                        gtgtGCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCag                 <  2:1177700/93‑1 (MQ=255)

                                                                          gtgCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCagag               <  2:2682880/93‑1 (MQ=255)

                                                                          gtgCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCagag               >  1:916351/1‑93 (MQ=255)

                                                                           tgCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAgg              >  1:2999216/1‑93 (MQ=255)

                                                                            gCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGc             <  2:9865607/93‑1 (MQ=255)

                                                                            gCTGGTGGTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGc             >  2:4757261/1‑93 (MQ=255)

                                                                              tggtggTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCATCGGGTCGCGCCGGGGCGt           >  2:4429502/1‑93 (MQ=255)

                                                                              tggtggTTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGt           >  2:7571463/1‑93 (MQ=255)

                                                                               ggtggtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGtt          >  2:4495821/1‑93 (MQ=255)

                                                                                gtggtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTta         >  2:1049894/1‑93 (MQ=255)

                                                                                gtggtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTta         >  2:228624/1‑93 (MQ=255)

                                                                                gtggtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTta         <  2:8424483/93‑1 (MQ=255)

                                                                                 tggtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTtat        >  2:1780561/1‑93 (MQ=255)

                                                                                  ggtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTtata       <  1:6395997/93‑1 (MQ=255)

                                                                                  ggtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTtata       >  1:6189728/1‑93 (MQ=255)

                                                                                  ggtTTCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTtata       <  2:565595/93‑1 (MQ=255)

                                                                                    tttCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTATATg     >  2:8706119/1‑93 (MQ=255)

                                                                                    tttCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTATATg     >  1:4758270/1‑93 (MQ=255)

                                                                                     ttCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTATATGa    <  1:9623223/93‑1 (MQ=255)

                                                                                      tCCCGGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTATATGAc   >  2:8522303/1‑93 (MQ=255)

                                                                                       cccGGTTTACCCTCGCCGCCGATACCGAACGGGGGATGCGCCCAAGTTTCTCCACAGGTGCATCACCGGATCGCGCAGAGGCGTTATCTGACt  >  1:2707273/1‑93 (MQ=255)

GTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGTTATATGACT  >  NC_000913/4074737‑4074916

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