BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,400,051	T→G	E346A (GAA→GCA)	mreB ←	cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,400,051	0	T	G	96.5%	501.4 / 5.8	146	E346A (GAA→GCA)	mreB	cell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand): ref base T (2/3); new base G (68/71); total (72/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 151 total aligned reads displayed.

GCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACT  >  NC_000913/3399955‑3400146

gCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCg                                                                                               >  2:1245625/1‑99 (MQ=255)

 ccACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAACCAGACACTTCCCCTGCCTGCATTCGATTACTCTgcgc                                                                                              >  2:3566048/1‑99 (MQ=255)

 ccACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTgcgc                                                                                              >  1:326763/1‑99 (MQ=255)

 ccACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACACACACTTCCCCTGCCTGCATCCGATTACTCTgcgc                                                                                              >  2:2465350/1‑99 (MQ=255)

 ccACGGCTAAAAATTGGCTTCATAAGTTATACGTATTCTCGTATCAGACCAGTCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTgcgc                                                                                              >  1:4354099/1‑99 (MQ=255)

  cACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCCGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCt                                                                                             <  2:3362400/99‑1 (MQ=255)

     ggCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGaa                                                                                          <  1:3520488/99‑1 (MQ=255)

     ggCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGaa                                                                                          <  1:3118329/99‑1 (MQ=255)

     ggCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGaa                                                                                          <  2:2649318/99‑1 (MQ=255)

     ggCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAAACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGaa                                                                                          >  1:825246/1‑99 (MQ=255)

        tAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAg                                                                                       >  1:2728537/1‑99 (MQ=255)

        tAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATAACTCTGCGCTGAACAg                                                                                       >  1:4010960/1‑99 (MQ=255)

         aaaaaTGGGCTTCATAAGTTATGCGTATTCGCGTATCAGACCAGTCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAgg                                                                                      <  1:2162623/99‑1 (MQ=255)

          aaaaTTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGt                                                                                     <  1:2508826/99‑1 (MQ=255)

            aaTTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTcg                                                                                   <  1:2068839/99‑1 (MQ=255)

            aaTTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTcg                                                                                   <  1:3366696/99‑1 (MQ=255)

            aaTTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTcg                                                                                   <  1:2381562/99‑1 (MQ=255)

             aTTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTcgc                                                                                  <  1:148370/99‑1 (MQ=255)

               tGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTcgccg                                                                                <  2:914536/99‑1 (MQ=255)

                 gCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCgccgcc                                                                              >  2:4130420/1‑99 (MQ=255)

                  cTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCTCCTGCCTGCATCCGATTAATCTGCGCTGAACAGGTCGccgacg                                                                             >  2:161991/1‑99 (MQ=255)

                  cTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGccgccg                                                                             <  1:1535187/99‑1 (MQ=255)

                   ttCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGt                                                                            >  2:1557978/1‑99 (MQ=255)

                     cATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGc                                                                          <  1:2208325/99‑1 (MQ=255)

                     cATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGc                                                                          >  1:1993652/1‑99 (MQ=255)

                     cATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGc                                                                          <  1:1695554/99‑1 (MQ=255)

                     cATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGc                                                                          <  2:3474399/99‑1 (MQ=255)

                      aTAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACGCTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCa                                                                         >  1:3888113/1‑99 (MQ=255)

                      aTAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCa                                                                         >  2:4152542/1‑99 (MQ=255)

                       tAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCAt                                                                        >  1:2594497/1‑99 (MQ=255)

                          gTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTc                                                                     <  1:4124847/99‑1 (MQ=255)

                            tATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGa                                                                   <  2:2172097/99‑1 (MQ=255)

                             aTGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGAt                                                                  >  1:2005704/1‑99 (MQ=255)

                              tGCGTATTCGCGTATCAGACCAGGCAGGGTAAACAGACGCTGCCCCTGCCTGGATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATc                                                                 <  2:2253825/99‑1 (MQ=255)

                               gCGTATTCGCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCa                                                                >  1:4239061/1‑99 (MQ=255)

                                cGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCAt                                                               >  2:1719948/1‑99 (MQ=255)

                                cGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCAt                                                               >  2:2036676/1‑99 (MQ=255)

                                 gTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCAtt                                                              >  1:3741246/1‑99 (MQ=255)

                                 gTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCAtt                                                              <  2:3870725/99‑1 (MQ=255)

                                   aTTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTc                                                            >  2:467438/1‑99 (MQ=255)

                                   aTTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTc                                                            >  2:3554392/1‑99 (MQ=255)

                                      ctcGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAg                                                         <  1:2859443/99‑1 (MQ=255)

                                       tcGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAgc                                                        >  2:1954359/1‑99 (MQ=255)

                                       tcGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAgc                                                        >  1:3785588/1‑99 (MQ=255)

                                         gTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAgcgc                                                      >  1:2223141/1‑99 (MQ=255)

                                          tATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCt                                                     <  1:3842220/99‑1 (MQ=255)

                                           aTCAGACCAGGCAGGGTAAACAGACACTTCCCCTACGTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCtt                                                    >  2:4168555/1‑99 (MQ=255)

                                            tCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCCATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCttt                                                   <  1:915863/99‑1 (MQ=255)

                                             cAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTg                                                  <  1:109322/99‑1 (MQ=255)

                                               gACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgcc                                                >  2:2579363/1‑99 (MQ=255)

                                               gACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgcc                                                >  2:1891131/1‑99 (MQ=255)

                                                aCCAGGCAGGGTAAACAGACACTTCCGCTGCCTGCATCCGATTACTCAGCGCTCAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgcag                                               >  2:4257181/1‑97 (MQ=255)

                                                aCCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccg                                               >  1:2503883/1‑99 (MQ=255)

                                                aCCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCGGCGCTTAACAGGTCGCCGCCTTGCATTTCGATCATTTCCATCGCTTTgccg                                               >  2:2789186/1‑99 (MQ=255)

                                                  caggcaggGTAAACAGACACTTCCCCTGCCTGCCTCCGATTACTCTCCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccgcc                                             >  2:2408842/1‑99 (MQ=255)

                                                  caggcaggGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccgcc                                             >  2:626570/1‑99 (MQ=255)

                                                         ggTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACcgcgc                                      <  2:1859476/99‑1 (MQ=255)

                                                         ggTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAAAAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACcgccc                                      >  2:4247917/1‑97 (MQ=255)

                                                          gTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCTATCATTTCCAGCGCTTTGCCGCCACcgcgcg                                     >  1:1466038/1‑99 (MQ=255)

                                                          gTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACcgcgcg                                     <  1:1537481/99‑1 (MQ=255)

                                                           tAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCgcgcgc                                    <  2:255762/99‑1 (MQ=255)

                                                           tAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCgcgcgc                                    >  1:608890/1‑99 (MQ=255)

                                                                aGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCcacac                               >  2:1438982/1‑99 (MQ=255)

                                                                aGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCcacac                               >  2:3940113/1‑99 (MQ=255)

                                                                 gACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCcacaca                              >  2:2374469/1‑99 (MQ=255)

                                                                  acacTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAg                             <  1:4032681/99‑1 (MQ=255)

                                                                  acacTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAg                             <  1:2793659/99‑1 (MQ=255)

                                                                    acTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATGTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGt                           >  2:3600876/1‑99 (MQ=255)

                                                                      ttCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCa                         >  1:2891345/1‑99 (MQ=255)

                                                                      ttCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCa                         >  2:2483986/1‑99 (MQ=255)

                                                                      ttCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATCTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTcc                         >  2:3445215/1‑98 (MQ=255)

                                                                        ccccTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGc                       >  2:3697784/1‑99 (MQ=255)

                                                                        ccccTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGc                       <  1:2866616/99‑1 (MQ=255)

                                                                         cccTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCg                      >  1:2072668/1‑99 (MQ=255)

                                                                          cctgcctgCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCgg                     <  2:3596629/99‑1 (MQ=255)

                                                                          cctgcctgCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCgg                     <  1:2780968/99‑1 (MQ=255)

                                                                             gcctgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTc                  <  2:3062859/99‑1 (MQ=255)

                                                                             gcctgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTc                  >  1:3945151/1‑99 (MQ=255)

                                                                              cctgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCt                 >  1:1065866/1‑99 (MQ=255)

                                                                               ctgcATACGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCtt                <  2:1009885/99‑1 (MQ=255)

                                                                                tgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTc               <  1:4064585/99‑1 (MQ=255)

                                                                                 gcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCa              <  1:2314159/99‑1 (MQ=255)

                                                                                  cATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAg             <  1:299979/99‑1 (MQ=255)

                                                                                  cATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAg             >  1:2205632/1‑99 (MQ=255)

                                                                                   aTCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGc            <  1:3652071/99‑1 (MQ=255)

                                                                                   aTCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGc            <  1:4196445/99‑1 (MQ=255)

                                                                                     ccGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGcaa          <  2:1552196/99‑1 (MQ=255)

                                                                                      cGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGcaac         <  1:3501027/99‑1 (MQ=255)

                                                                                        aTTACTCTGCGCTGAACAGGTCGCCGCCGTGCGTGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGcaacaa       <  2:3717134/99‑1 (MQ=255)

                                                                                         ttgctctGCGCTCAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCGTCAGCaacaac      <  2:3358751/96‑1 (MQ=255)

                                                                                         ttACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCaacaac      <  1:1941661/99‑1 (MQ=255)

                                                                                          tACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACg     >  1:4204451/1‑99 (MQ=255)

                                                                                            ctctGCGCTGAACAGTTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGAc   >  2:4127806/1‑99 (MQ=255)

                                                                                            ctctGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGAc   >  1:1315263/1‑99 (MQ=255)

                                                                                            ctctGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGAc   >  1:1211464/1‑99 (MQ=255)

                                                                                            ctctGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGAc   <  1:2947761/99‑1 (MQ=255)

                                                                                            ctctGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGAc   <  2:3085959/99‑1 (MQ=255)

                                                                                             tctGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACt  >  1:991011/1‑99 (MQ=255)

                                                                                             tctGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACt  >  1:2891889/1‑99 (MQ=255)

                                                                                             tctGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACt  >  1:3219051/1‑99 (MQ=255)

GCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACT  >  NC_000913/3399955‑3400146

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