BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,857,508	T→C	E79G (GAA→GGA)	yidK ←	putative cotransporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,857,508	0	T	C	99.6%	1508.3 / ‑3.5	467	E79G (GAA→GGA)	yidK	putative cotransporter
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (234/231); total (235/232)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 481 total aligned reads displayed.

AGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA  >  NC_000913/3857418‑3857597

aGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTccc                                                                                         >  2:748472/1‑93 (MQ=255)

  aaaaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCa                                                                                       <  2:1441302/93‑1 (MQ=255)

  aaaaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCa                                                                                       >  1:7748259/1‑93 (MQ=255)

   aactCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAg                                                                                      <  1:1731796/90‑1 (MQ=255)

    aaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGc                                                                                     >  1:7901691/1‑93 (MQ=255)

    aaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGc                                                                                     <  1:3899973/93‑1 (MQ=255)

      tCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGccc                                                                                   >  2:1201715/1‑93 (MQ=255)

        ggggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCAt                                                                                 >  1:6161808/1‑93 (MQ=255)

         gggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATc                                                                                >  2:695957/1‑93 (MQ=255)

          ggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATcc                                                                               >  2:7247050/1‑92 (MQ=255)

          ggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCa                                                                               <  1:9727537/93‑1 (MQ=255)

           gATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCAc                                                                              <  2:4476645/93‑1 (MQ=255)

            aTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCAcc                                                                             <  1:1043281/93‑1 (MQ=255)

            aTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCAcc                                                                             <  2:9290349/93‑1 (MQ=255)

            aTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGCCACCGCTGAAGTCACTCCCCAGCCCATCAcc                                                                             <  1:7280346/93‑1 (MQ=255)

             tggtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCg                                                                            >  1:6597811/1‑93 (MQ=255)

             tggtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCg                                                                            >  2:6110444/1‑93 (MQ=255)

              ggtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGa                                                                           >  2:8404877/1‑93 (MQ=255)

               gtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGAc                                                                          <  2:6268397/93‑1 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:3407078/93‑1 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  2:5047159/93‑1 (MQ=255)

                tgacAATTCCCCGCTTCAGATAGCCCCGTAAAAAGATTAGCGCGAGGAAGATCAGCTTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:2460520/90‑1 (MQ=255)

                   gaaTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGc                                                                      <  2:7545456/92‑1 (MQ=255)

                    aaTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGcc                                                                     <  1:8327477/93‑1 (MQ=255)

                     aTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCg                                                                    >  1:4278371/1‑93 (MQ=255)

                      tGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGc                                                                   >  1:9870915/1‑93 (MQ=255)

                         cgcgCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCttt                                                                >  1:9268976/1‑93 (MQ=255)

                         cgcgCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCttt                                                                >  2:9403253/1‑93 (MQ=255)

                          gcgcTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCtttt                                                               >  2:6537241/1‑93 (MQ=255)

                             cTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTa                                                            >  1:4902934/1‑93 (MQ=255)

                             cTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTa                                                            <  2:2582041/93‑1 (MQ=255)

                              ttCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAg                                                           <  2:4534258/93‑1 (MQ=255)

                               tCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAgg                                                          >  2:3006512/1‑93 (MQ=255)

                                cAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCCCTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGc                                                         <  2:9825421/93‑1 (MQ=255)

                                 agagagCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGcc                                                        <  2:4653319/93‑1 (MQ=255)

                                  gATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCt                                                       >  1:8875910/1‑93 (MQ=255)

                                  gATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCt                                                       >  2:7148080/1‑93 (MQ=255)

                                     aGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGcc                                                    <  2:6395309/93‑1 (MQ=255)

                                        gcgGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCgg                                                 >  1:6385618/1‑93 (MQ=255)

                                          ggTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaa                                               <  1:3667990/93‑1 (MQ=255)

                                          ggTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaa                                               <  1:4671025/93‑1 (MQ=255)

                                          ggTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGACCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaa                                               >  1:5336761/1‑93 (MQ=255)

                                           gTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaaa                                              >  1:4575029/1‑93 (MQ=255)

                                           gTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaaa                                              >  2:7872893/1‑93 (MQ=255)

                                            tAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAg                                             <  1:8773651/93‑1 (MQ=255)

                                            tAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAg                                             <  2:9524355/93‑1 (MQ=255)

                                            tAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAg                                             <  1:1874235/93‑1 (MQ=255)

                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAgg                                            >  1:4869464/1‑93 (MQ=255)

                                              aaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGc                                           >  1:8583139/1‑93 (MQ=255)

                                              aaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGc                                           <  2:3634398/93‑1 (MQ=255)

                                              aaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGc                                           <  2:2244579/93‑1 (MQ=255)

                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCGTCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGcc                                          >  2:5875747/1‑93 (MQ=255)

                                                aaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCg                                         >  2:1620700/1‑93 (MQ=255)

                                                aaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCg                                         >  1:6880628/1‑93 (MQ=255)

                                                aaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCg                                         >  1:149235/1‑93 (MQ=255)

                                                 agcttAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGa                                        >  2:4297211/4‑93 (MQ=255)

                                                 aGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGa                                        >  2:3540585/1‑93 (MQ=255)

                                                 aGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGa                                        >  1:10043121/1‑93 (MQ=255)

                                                  gATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGAc                                       >  2:6049957/1‑93 (MQ=255)

                                                     tAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAg                                    >  2:6904259/1‑93 (MQ=255)

                                                       gcgcGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgtt                                  <  1:2891304/93‑1 (MQ=255)

                                                          cgAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgttgtt                               <  2:8860894/93‑1 (MQ=255)

                                                           gAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTc                              >  2:3780282/1‑93 (MQ=255)

                                                           gAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTc                              <  2:4071818/93‑1 (MQ=255)

                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             >  1:9558919/1‑93 (MQ=255)

                                                             ggAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCg                            >  2:859196/1‑93 (MQ=255)

                                                             ggAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCg                            >  1:898051/1‑93 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           <  2:7370673/93‑1 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           >  2:5598277/1‑93 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           <  1:6811438/93‑1 (MQ=255)

                                                               aaGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTg                          >  1:5417794/1‑93 (MQ=255)

                                                                 gATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGa                        >  1:982835/1‑93 (MQ=255)

                                                                   tCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaaa                      <  1:8100427/93‑1 (MQ=255)

                                                                     aGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAgg                    >  2:9108920/1‑93 (MQ=255)

                                                                      gCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   <  1:1208263/93‑1 (MQ=255)

                                                                      gCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   <  2:5403049/93‑1 (MQ=255)

                                                                       cctCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAgggt                  <  1:2134792/91‑3 (MQ=255)

                                                                       cGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggtt                  <  2:3415476/93‑1 (MQ=255)

                                                                        gTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttg                 >  1:2770328/1‑93 (MQ=255)

                                                                         tCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttgg                <  2:1027708/93‑1 (MQ=255)

                                                                         tCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttgg                >  2:9648214/1‑93 (MQ=255)

                                                                         tCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttgg                >  2:8108994/1‑93 (MQ=255)

                                                                          cACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttggt               <  2:3284809/93‑1 (MQ=255)

                                                                          cACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttggt               <  1:7492923/93‑1 (MQ=255)

                                                                           cccGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttggtt              <  1:5167710/92‑1 (MQ=255)

                                                                           cccGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttggtt              <  2:3170765/92‑1 (MQ=255)

                                                                           aCCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttggtt              <  1:1330892/93‑1 (MQ=255)

                                                                            ccGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTa             <  1:3023823/93‑1 (MQ=255)

                                                                              gCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGc           >  2:7699203/1‑93 (MQ=255)

                                                                               cTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCa          <  2:673509/93‑1 (MQ=255)

                                                                                tGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         >  1:2233046/1‑93 (MQ=255)

                                                                                 gAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGCAAGGCCGACCAGCTGTTCCGTGGAAAGGTTGGTTAGCAtt        >  1:5554179/1‑93 (MQ=255)

                                                                                 gAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAtt        >  1:7103258/1‑93 (MQ=255)

                                                                                  aaGTCACTCCCCAGCCCATCACCGACATGGCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTa       >  2:8944545/1‑93 (MQ=255)

                                                                                  aaGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTa       >  2:2055053/1‑93 (MQ=255)

                                                                                  aaGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTa       <  2:9214049/93‑1 (MQ=255)

                                                                                   aGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTaa      <  1:3268200/93‑1 (MQ=255)

                                                                                   aGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTaa      >  1:5666501/1‑93 (MQ=255)

                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     >  2:8220467/1‑93 (MQ=255)

                                                                                       aCTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGaa  <  1:4025449/93‑1 (MQ=255)

AGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA  >  NC_000913/3857418‑3857597

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