BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,395,834	C→G	A282G (GCC→GGC)	mutL →	MutL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,395,834	0	C	G	93.4%	526.9 / 24.5	181	A282G (GCC→GGC)	mutL	MutL
Reads supporting (aligned to +/- strand): ref base C (7/5); new base G (85/84); total (92/89)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

Only 100 of 196 total aligned reads displayed.

GAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTT  >  NC_000913/4395743‑4395923

gAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGcc                                                                                          >  1:3256230/1‑93 (MQ=255)

    ttCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATc                                                                                      >  1:3908792/1‑93 (MQ=255)

    ttCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATc                                                                                      <  2:17802/93‑1 (MQ=255)

       aGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGc                                                                                   >  1:2555099/1‑93 (MQ=255)

        gTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCa                                                                                  >  2:3276445/1‑93 (MQ=255)

        gTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCAAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCa                                                                                  >  1:1890621/1‑93 (MQ=255)

          aTTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAAc                                                                                <  1:3420598/93‑1 (MQ=255)

           ttGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAAcc                                                                               <  2:2362545/93‑1 (MQ=255)

           ttGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAAcc                                                                               >  2:3574030/1‑93 (MQ=255)

            tGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCg                                                                              <  2:291497/93‑1 (MQ=255)

              cTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGc                                                                            >  1:4436901/1‑93 (MQ=255)

               tACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCa                                                                           >  2:3674741/1‑93 (MQ=255)

               tACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCa                                                                           >  1:3954738/1‑93 (MQ=255)

               tACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCa                                                                           >  2:2030693/1‑93 (MQ=255)

               tACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCa                                                                           >  2:151944/1‑93 (MQ=255)

                 cGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCAtt                                                                         <  1:4170722/93‑1 (MQ=255)

                   tGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGTAGACAAACTAAGGGAAGATAAACAAGCGGGATTtg                                                                       >  2:4101927/1‑93 (MQ=255)

                   tGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTtg                                                                       <  1:1690292/93‑1 (MQ=255)

                   tGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTtg                                                                       <  1:1509208/93‑1 (MQ=255)

                    gAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTtgt                                                                      <  1:2940893/93‑1 (MQ=255)

                    gAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTtgt                                                                      <  2:4275573/93‑1 (MQ=255)

                    gAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTtgt                                                                      <  2:1183087/93‑1 (MQ=255)

                     aaCGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTtgtg                                                                     <  2:2708232/93‑1 (MQ=255)

                      aCGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTgtgt                                                                    <  2:1239075/93‑1 (MQ=255)

                      aCGGTCCCATGATGCGCGACCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTgtgt                                                                    <  2:2997737/93‑1 (MQ=255)

                       cGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGtgtt                                                                   >  2:500162/1‑93 (MQ=255)

                        ggTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGtgttg                                                                  <  1:2124052/93‑1 (MQ=255)

                        ggTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGtgttg                                                                  >  2:2113805/1‑93 (MQ=255)

                         gTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGtgttgt                                                                 >  2:2449081/1‑93 (MQ=255)

                          tCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTa                                                                <  1:846324/93‑1 (MQ=255)

                          tCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTa                                                                <  1:3824381/93‑1 (MQ=255)

                           cGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTAt                                                               >  1:2698464/1‑93 (MQ=255)

                            gCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATc                                                              >  1:2915332/1‑93 (MQ=255)

                            ccATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATc                                                              <  2:699078/92‑1 (MQ=255)

                             cATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCt                                                             <  1:2443318/93‑1 (MQ=255)

                             cATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCt                                                             >  1:2018179/1‑93 (MQ=255)

                             cATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCt                                                             <  2:3956224/93‑1 (MQ=255)

                              atgatgCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTg                                                            >  1:2291587/1‑93 (MQ=255)

                              atgatgCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTg                                                            >  1:3601576/1‑93 (MQ=255)

                              atgatgCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTg                                                            >  1:3549672/1‑93 (MQ=255)

                              atgatgCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTg                                                            >  2:2801838/1‑93 (MQ=255)

                               tgatgCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTgg                                                           >  1:2670504/1‑93 (MQ=255)

                                gatgCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGga                                                          <  2:746713/93‑1 (MQ=255)

                                 atgCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGgag                                                         <  1:127042/93‑1 (MQ=255)

                                    cgcgATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATc                                                      <  2:2318339/93‑1 (MQ=255)

                                    cgcgATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATc                                                      <  2:4239595/93‑1 (MQ=255)

                                      cgATCGCCTGATCAATCACGCGATCCTCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGa                                                    <  1:1962395/93‑1 (MQ=255)

                                       gATCGCCTGATCAATCAGGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGAc                                                   <  1:4226025/93‑1 (MQ=255)

                                       gATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGAc                                                   <  2:4625494/93‑1 (MQ=255)

                                       gATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGAc                                                   >  2:2528213/1‑93 (MQ=255)

                                        aTCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGAcc                                                  >  2:24331/1‑93 (MQ=255)

                                        aTCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGAcc                                                  <  1:1621589/93‑1 (MQ=255)

                                        aTCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGAcc                                                  <  2:1061658/93‑1 (MQ=255)

                                           gCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCcac                                               >  2:2626621/1‑93 (MQ=255)

                                              tGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATc                                            >  2:2749681/1‑93 (MQ=255)

                                              tGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATc                                            >  1:1785701/1‑93 (MQ=255)

                                              tGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATc                                            <  2:3450407/93‑1 (MQ=255)

                                               gATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCa                                           >  2:3158307/1‑93 (MQ=255)

                                               gATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCa                                           >  1:321669/1‑93 (MQ=255)

                                                  caatcaCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGt                                        <  1:3944979/93‑1 (MQ=255)

                                                  caaccaCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGt                                        <  1:3987020/93‑1 (MQ=255)

                                                   aatcaCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTg                                       <  2:4440793/93‑1 (MQ=255)

                                                    atcaCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTgg                                      >  2:3236313/1‑93 (MQ=255)

                                                     tcaCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGa                                     >  1:4544444/1‑93 (MQ=255)

                                                       aCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACg                                   >  1:2170226/1‑93 (MQ=255)

                                                        cgcgATCCGCCAGGCCTGCGAAGACAAACTGGGAGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGt                                  <  1:647629/93‑1 (MQ=255)

                                                         gcgATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTc                                 <  1:11246/93‑1 (MQ=255)

                                                         gcgATCCGCCAGGCCCGGGAAGACAAACTGGCGGGCGATCAGCAACCGGCATTTGTGGTGTATCTGGAGATCGACACACATCACGTGGACgcc                                 >  2:563910/1‑91 (MQ=255)

                                                          ggATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAGCCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCa                                <  1:2636457/92‑1 (MQ=255)

                                                           gATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCaa                               >  1:3716275/1‑93 (MQ=255)

                                                           gATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCaa                               <  1:2346336/93‑1 (MQ=255)

                                                           gATCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCaa                               >  2:483103/1‑93 (MQ=255)

                                                            aTCCGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAAc                              >  1:620627/1‑93 (MQ=255)

                                                              ccGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGt                            >  1:3063343/1‑93 (MQ=255)

                                                              ccGCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGt                            <  1:4143914/93‑1 (MQ=255)

                                                                gCCAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGc                          >  2:1950912/1‑93 (MQ=255)

                                                                  cAGGCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCAc                        >  2:3680191/1‑93 (MQ=255)

                                                                    ggCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCAccc                      <  2:1646355/93‑1 (MQ=255)

                                                                    ggCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCAccc                      >  1:1818698/1‑93 (MQ=255)

                                                                    ggCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCAccc                      >  1:3632813/1‑93 (MQ=255)

                                                                     gCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCAcccc                     <  1:2313992/93‑1 (MQ=255)

                                                                     gCCTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCAcccc                     <  1:1735237/93‑1 (MQ=255)

                                                                       cTGCGAAGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGc                   >  2:155038/1‑93 (MQ=255)

                                                                            aaGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAac              <  2:404895/93‑1 (MQ=255)

                                                                             aGACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAaca             <  1:2479047/93‑1 (MQ=255)

                                                                              gACAAACTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAacac            >  2:640985/1‑93 (MQ=255)

                                                                                 aaaCTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGac         >  2:3525723/1‑92 (MQ=255)

                                                                                 aaaCTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGaa         <  2:1321932/93‑1 (MQ=255)

                                                                                 aaaCTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGaa         <  1:1389513/93‑1 (MQ=255)

                                                                                 aaaCTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGaa         <  1:3972653/93‑1 (MQ=255)

                                                                                  aaCTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAg        <  2:88890/93‑1 (MQ=255)

                                                                                  aaCTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAg        <  1:286579/93‑1 (MQ=255)

                                                                                  aaCTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAg        <  1:4334712/93‑1 (MQ=255)

                                                                                   aCTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAgg       >  1:3448052/1‑92 (MQ=255)

                                                                                   aCTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGt       <  1:2771263/93‑1 (MQ=255)

                                                                                   aCTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGt       <  2:2183191/93‑1 (MQ=255)

                                                                                    cTGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTg      <  2:2864789/93‑1 (MQ=255)

                                                                                     tGGGGGGCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCCAACACGAAGTGc     >  2:4340850/1‑93 (MQ=255)

                                                                                     tGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGc     <  1:1369630/93‑1 (MQ=255)

                                                                                        ggggCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGtt  <  1:525438/93‑1 (MQ=255)

GAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTT  >  NC_000913/4395743‑4395923

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