BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,210,781	T→C	M68V (ATG→GTG)	b2125 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,210,781	0	T	C	99.3%	445.3 / ‑4.5	138	M68V (ATG→GTG)	b2125	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (75/62); total (76/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

Only 100 of 145 total aligned reads displayed.

TCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCATCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATCCC  >  NC_000913/2210693‑2210868

tCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACcccc                                                                                     >  1:2995301/1‑93 (MQ=255)

acAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACcccc                                                                                     <  1:1447199/92‑1 (MQ=255)

 cAAAGTCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCcca                                                                                    <  2:108212/93‑1 (MQ=255)

   aaGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacc                                                                                  >  2:941052/1‑93 (MQ=255)

   aaGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacc                                                                                  >  2:221680/1‑93 (MQ=255)

   aaGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacc                                                                                  >  2:297432/1‑93 (MQ=255)

    aGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacca                                                                                 >  2:2032488/1‑93 (MQ=255)

     ggCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCAt                                                                                >  1:4119407/1‑93 (MQ=255)

      gCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCAtt                                                                               >  1:2872367/1‑93 (MQ=255)

        aTGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTc                                                                             >  2:1043865/1‑93 (MQ=255)

        aTGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTc                                                                             >  2:3596724/1‑93 (MQ=255)

         tGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTcc                                                                            >  1:4200848/1‑93 (MQ=255)

         tGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTcc                                                                            <  2:2607111/93‑1 (MQ=255)

          gTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCa                                                                           <  1:857487/93‑1 (MQ=255)

            tcttcAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGa                                                                         >  2:1930785/1‑93 (MQ=255)

            tcttcAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTACAGa                                                                         >  2:1930915/1‑93 (MQ=255)

              ttcAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGAcc                                                                       >  2:3059368/1‑93 (MQ=255)

                cAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCAc                                                                     >  2:2322992/1‑93 (MQ=255)

                 aaaGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACt                                                                    >  2:2566618/1‑93 (MQ=255)

                 aaaGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACt                                                                    >  2:2404005/1‑93 (MQ=255)

                  aaGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTg                                                                   <  1:2002074/93‑1 (MQ=255)

                   aGGCTTTAATTGCTTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGa                                                                  <  2:3551307/93‑1 (MQ=255)

                   aGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGa                                                                  <  2:3125937/93‑1 (MQ=255)

                   aGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGa                                                                  <  1:3018812/93‑1 (MQ=255)

                   aGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGa                                                                  <  1:1904404/93‑1 (MQ=255)

                    ggCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGAt                                                                 <  2:3405577/93‑1 (MQ=255)

                     gCTTTAATTGCGTATTCGTCAAACGCAGCGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATg                                                                >  2:3652110/1‑93 (MQ=255)

                      cTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgc                                                               <  1:3842076/93‑1 (MQ=255)

                      cTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgc                                                               >  1:3902225/1‑93 (MQ=255)

                       tttAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgcg                                                              >  2:2103814/1‑93 (MQ=255)

                        ttAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgcgc                                                             >  2:233067/1‑93 (MQ=255)

                        ttAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgcgc                                                             <  1:2423226/93‑1 (MQ=255)

                         tAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGcgcg                                                            >  1:3233031/1‑93 (MQ=255)

                          aaTTGCGGATTCGTTAAACCCAGGGAGGAAAACAATATACGGGCGGTGTTCCGGGTCCAGCACCCCCACCATTTCCAGACCACTGATGCGCgg                                                           >  2:2195122/1‑93 (MQ=255)

                           aTTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGc                                                          >  1:321784/1‑93 (MQ=255)

                             tGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCAt                                                        >  2:2023715/1‑93 (MQ=255)

                             tGCGTATTCGAAAAAAGCAGTGAAAAAAAAAATATGGGGGCGATGTCCCGGGTCAAGCCCCCCCACCATTTCCAGACCACTGATGCGCGGCAt                                                        <  2:4653124/93‑1 (MQ=255)

                              gCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATc                                                       <  2:2954877/93‑1 (MQ=255)

                              gCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATc                                                       <  1:1097885/93‑1 (MQ=255)

                                gTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTg                                                     <  2:463657/93‑1 (MQ=255)

                                gTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTg                                                     <  2:1062084/93‑1 (MQ=255)

                                gTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTg                                                     >  1:1806280/1‑93 (MQ=255)

                                 tATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTgg                                                    >  2:946027/1‑93 (MQ=255)

                                 tATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTgg                                                    <  1:1090712/93‑1 (MQ=255)

                                 tATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTgg                                                    <  1:2798063/93‑1 (MQ=255)

                                     cGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGatat                                                >  1:1027629/1‑93 (MQ=255)

                                      gTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATc                                               <  1:2094656/93‑1 (MQ=255)

                                      gTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATc                                               >  1:275708/1‑93 (MQ=255)

                                       tCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCg                                              >  1:4628004/1‑93 (MQ=255)

                                       tCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCg                                              >  2:4131365/1‑93 (MQ=255)

                                       tCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCg                                              >  2:2852582/1‑93 (MQ=255)

                                         aaaCGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCgag                                            >  2:2478169/1‑93 (MQ=255)

                                          aaCGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCgaga                                           <  2:4521503/93‑1 (MQ=255)

                                          aaCGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCgaga                                           <  1:3467264/93‑1 (MQ=255)

                                          aaCGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGCCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCgaga                                           >  2:4200200/1‑93 (MQ=255)

                                          aaCGCAGTGAGAAAAACAATATACGGGAGATGTTCAGGGTCAAGAACACCCACCATTTCCAGACCACTGATGCGCGCCATCTGGATATCgagt                                           >  2:4001472/1‑92 (MQ=255)

                                           aCGCAGTGAGAAAAACAAGATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGaa                                          <  1:1157986/93‑1 (MQ=255)

                                             gCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAAc                                        <  1:2232961/93‑1 (MQ=255)

                                              cAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACa                                       <  2:4417999/93‑1 (MQ=255)

                                              cAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACa                                       <  1:176827/93‑1 (MQ=255)

                                              cAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACa                                       >  2:3241466/1‑93 (MQ=255)

                                                gTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGc                                     <  1:3097762/93‑1 (MQ=255)

                                                  gagaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGcac                                   >  2:4730952/1‑93 (MQ=255)

                                                     aaaaaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATc                                <  1:116464/93‑1 (MQ=255)

                                                        aaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCgg                             <  1:356329/93‑1 (MQ=255)

                                                        aaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCgg                             >  1:634478/1‑93 (MQ=255)

                                                          cAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGgc                           <  1:4180427/93‑1 (MQ=255)

                                                            atatACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGgcgc                         >  2:630612/1‑93 (MQ=255)

                                                            atatACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGgcgc                         <  1:4084309/93‑1 (MQ=255)

                                                            atatACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGgcgc                         >  1:3136264/1‑93 (MQ=255)

                                                             tataCGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCa                        >  2:2004579/1‑93 (MQ=255)

                                                              ataCGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAg                       >  1:1332677/1‑93 (MQ=255)

                                                                aCGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCGGGATATCGAGAAACAGCACATCCGGGCGCAGtt                     <  2:3517797/93‑1 (MQ=255)

                                                                  gggCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCCGTTTa                   >  1:4174075/1‑93 (MQ=255)

                                                                  gggCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTa                   >  1:3388727/1‑93 (MQ=255)

                                                                   ggCGTTGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCTGGGCGTAGTTTag                  >  2:4044146/1‑92 (MQ=255)

                                                                     cGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGc                <  1:7726/93‑1 (MQ=255)

                                                                      gATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCa               <  2:1747481/93‑1 (MQ=255)

                                                                      gATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCa               <  2:559000/93‑1 (MQ=255)

                                                                       aTGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCAc              <  2:1133578/93‑1 (MQ=255)

                                                                       aTGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGAAc              >  1:3120664/1‑93 (MQ=255)

                                                                       aTGTTCCGGGTCAAGCACCCCCACCATTGCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCAc              <  1:458671/93‑1 (MQ=255)

                                                                        tGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCAcc             >  1:2828525/1‑93 (MQ=255)

                                                                        tGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCAcc             <  1:2551790/93‑1 (MQ=255)

                                                                        tGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCAcc             <  2:3855957/93‑1 (MQ=255)

                                                                         gTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACcg            >  2:200745/1‑93 (MQ=255)

                                                                          ttCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACcgc           <  2:2849759/93‑1 (MQ=255)

                                                                          ttCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACcgc           >  2:2732126/1‑93 (MQ=255)

                                                                          ttCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACcgc           <  1:2607587/93‑1 (MQ=255)

                                                                          ttCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACcgc           >  1:1601012/1‑93 (MQ=255)

                                                                            ccGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCgcgc         >  1:1716183/1‑93 (MQ=255)

                                                                            ccGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCgcgc         <  2:1794178/93‑1 (MQ=255)

                                                                            ccGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCgcgc         <  1:1496767/93‑1 (MQ=255)

                                                                             cGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGcc        <  2:583731/93‑1 (MQ=255)

                                                                             cGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGcc        <  1:569981/93‑1 (MQ=255)

                                                                              gggTCAAGCCCCCCCACCATTTCCAGACCACTCATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCCCCGCGCCg       <  2:1243017/93‑1 (MQ=255)

                                                                              gggTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCg       >  1:3937954/1‑93 (MQ=255)

                                                                                gTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGAt     <  2:3940084/93‑1 (MQ=255)

                                                                                 tCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATc    <  1:129781/93‑1 (MQ=255)

                                                                                   aaGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATccc  >  1:2865136/1‑93 (MQ=255)

TCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCATCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATCCC  >  NC_000913/2210693‑2210868

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