BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,765,234	G→A	G129G (GGG→GGA)	yibJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,765,234	0	G	A	99.5%	707.3 / ‑5.6	208	G129G (GGG→GGA)	yibJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (105/102); total (105/103)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 213 total aligned reads displayed.

CGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATACCC  >  NC_000913/3765184‑3765323

cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                                                 <  1:6073452/93‑1 (MQ=9)

cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                                                 <  2:1046273/93‑1 (MQ=9)

cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                                                 <  2:2604443/93‑1 (MQ=9)

 gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                                                <  1:1217314/93‑1 (MQ=12)

 gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                                                >  1:6820862/1‑93 (MQ=12)

 gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                                                >  2:1704163/1‑93 (MQ=12)

 gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                                                >  2:5075685/1‑93 (MQ=12)

  cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                                               >  2:5375293/1‑93 (MQ=255)

  cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                                               <  1:7474054/93‑1 (MQ=16)

   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGGTGGAATTTTTATCAGTATCCGTTGAATCCGGTc                                              >  2:7773312/1‑93 (MQ=255)

   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc                                              >  1:2796592/1‑93 (MQ=255)

    attatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTGTCAGTATCCGTTGAATCCGGTCa                                             >  2:4359280/1‑93 (MQ=11)

    attatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCa                                             >  2:2268079/1‑93 (MQ=255)

    attatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCa                                             >  1:386633/1‑93 (MQ=255)

     ttatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTAGCAGTATCCGTTGAATCCGGTCAt                                            >  2:4301765/1‑93 (MQ=11)

       atGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaa                                          >  1:620520/1‑93 (MQ=255)

       atGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaa                                          <  2:3019787/93‑1 (MQ=255)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                         <  1:4040128/93‑1 (MQ=255)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                         <  1:1931944/93‑1 (MQ=255)

          aCCCGTTGCTGGGGCGGTAGATCACTCAGTATCCGATTGGACTGAAGTGGGGATGGAATTTTTATCAGTAGCCGTTGAATCCGGTCATAAATg                                       >  1:7019772/1‑93 (MQ=14)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                                    >  1:3457640/1‑93 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                                    >  2:3503463/1‑93 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                                    <  1:2026284/93‑1 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                                    >  1:2228633/1‑93 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                                    >  2:5558204/1‑93 (MQ=255)

              gctgctGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa                                   <  2:4042918/93‑1 (MQ=255)

              gctgctGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa                                   <  1:3381836/93‑1 (MQ=255)

               ctgctgGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt                                  >  2:1344903/1‑93 (MQ=255)

               ctgctgGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt                                  >  2:7740911/1‑93 (MQ=255)

                tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                                 <  1:3666745/93‑1 (MQ=255)

                tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                                 <  1:3513674/93‑1 (MQ=255)

                tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                                 >  2:3759289/1‑93 (MQ=255)

                tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                                 <  2:5411555/93‑1 (MQ=255)

                 gctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATcc                                <  1:4959683/90‑1 (MQ=255)

                 gctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATcc                                <  2:4981011/90‑1 (MQ=255)

                   tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCCTAAATGTAGATCCGc                              >  2:3209395/2‑93 (MQ=255)

                   tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGc                              <  2:1474791/92‑1 (MQ=255)

                   tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGc                              <  2:4239957/92‑1 (MQ=255)

                   tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGc                              >  2:4205153/2‑93 (MQ=255)

                   tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGc                              >  2:1408530/2‑93 (MQ=255)

                     gggCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                            <  2:6910576/93‑1 (MQ=255)

                     gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                            >  2:3794392/1‑93 (MQ=255)

                     gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                            <  1:2460419/93‑1 (MQ=255)

                     gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                            <  1:7604875/93‑1 (MQ=255)

                     gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                            <  1:1852188/93‑1 (MQ=255)

                     gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                            <  2:6520793/93‑1 (MQ=255)

                      ggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAg                           <  2:3690020/93‑1 (MQ=255)

                      ggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAg                           >  2:3630443/1‑93 (MQ=255)

                      ggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAg                           >  1:3569859/1‑93 (MQ=255)

                       gCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAgg                          >  1:6473672/1‑93 (MQ=255)

                       gCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTAATAAATGTAGATCCGCAAgg                          >  1:5932433/1‑93 (MQ=255)

                         ggTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGtt                        >  2:2205251/1‑93 (MQ=255)

                         ggTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGtt                        >  2:3589343/1‑93 (MQ=255)

                         ggTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGtt                        >  2:293285/1‑93 (MQ=255)

                          gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                       >  1:5569336/1‑93 (MQ=255)

                          gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                       <  1:91798/93‑1 (MQ=255)

                          gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                       <  1:5486887/93‑1 (MQ=255)

                          gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                       >  1:1888532/1‑93 (MQ=255)

                          gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                       <  1:536309/93‑1 (MQ=255)

                          gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                       >  1:4385133/1‑93 (MQ=255)

                           tataTCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg                      >  2:1605811/1‑93 (MQ=255)

                            atatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTgg                     <  2:5304212/93‑1 (MQ=255)

                            atatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTgg                     >  1:2939210/1‑93 (MQ=255)

                            atatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCAGCAAGGTTTgg                     >  1:6115680/1‑93 (MQ=255)

                             tatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt                    >  2:2331439/1‑93 (MQ=255)

                             tatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt                    >  1:4475438/1‑93 (MQ=255)

                              atCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGtt                   >  1:7535621/1‑93 (MQ=255)

                               tCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTg                  >  2:3884980/1‑93 (MQ=255)

                                 aCTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGat                >  2:2702570/1‑93 (MQ=255)

                                  cTCAGGATCCGATTGGCCTAAAGGGGGGATGGAATTTTTATCAGTATACGTTGAATCCGGTCATAAATGTAGATCCGCACGGATTGGTTGata               <  2:7031939/93‑1 (MQ=18)

                                  cTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGata               <  2:7806038/93‑1 (MQ=255)

                                  cTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGata               <  1:122674/93‑1 (MQ=255)

                                  cTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGata               <  1:1783278/93‑1 (MQ=255)

                                  cTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGata               >  2:1659348/1‑93 (MQ=255)

                                   tCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGatat              <  1:4007779/93‑1 (MQ=255)

                                   tCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGatat              >  1:4782902/1‑93 (MQ=255)

                                   tCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGatat              <  1:4299020/93‑1 (MQ=255)

                                   tCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGatat              <  2:7308878/93‑1 (MQ=255)

                                   tCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGatat              >  2:2023607/1‑93 (MQ=255)

                                    cAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGAtata             <  2:6157171/93‑1 (MQ=255)

                                      ggATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATaaa           <  1:5812714/93‑1 (MQ=255)

                                       gATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAt          >  1:4936336/1‑93 (MQ=255)

                                       gATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAt          <  1:7383310/93‑1 (MQ=255)

                                        aTCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAtt         <  2:4747807/93‑1 (MQ=255)

                                         tCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAttt        >  1:2747279/1‑93 (MQ=255)

                                          ccGATTGGACTGAAGGGGGTATGGAATTTTTATCAGTATCCGTGGAATCCGGTCATAAATGGAGATCCGCAAGGTTTGGTTGATATAAATTta       <  2:72399/93‑1 (MQ=25)

                                          ccGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTta       <  1:61565/93‑1 (MQ=255)

                                          ccGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTta       <  2:5342541/93‑1 (MQ=255)

                                          ccGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTta       <  1:2574228/93‑1 (MQ=255)

                                          ccGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCAGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTta       >  1:1972994/1‑93 (MQ=255)

                                           cGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTtat      <  2:145758/93‑1 (MQ=255)

                                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTtata     >  2:5930198/1‑93 (MQ=255)

                                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTtata     <  1:1073819/93‑1 (MQ=255)

                                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTtata     >  2:6744233/1‑93 (MQ=255)

                                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTtata     <  2:1758012/93‑1 (MQ=255)

                                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTAGCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTGTGGTTGAAATAAATTtata     <  1:1677580/93‑1 (MQ=25)

                                             aTTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATAc    >  2:7562075/1‑93 (MQ=255)

                                             aTTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATAc    <  2:769239/93‑1 (MQ=255)

                                              ttGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATAcc   <  2:2383786/93‑1 (MQ=255)

                                              ttGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATAcc   >  1:7176127/1‑93 (MQ=255)

CGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATACCC  >  NC_000913/3765184‑3765323

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: `‑`