BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,975,051	G→A	A38V (GCC→GTC)	motA ←	MotA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,975,051	0	G	A	99.3%	496.7 / ‑5.3	144	A38V (GCC→GTC)	motA	MotA
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (83/59); total (84/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 149 total aligned reads displayed.

GAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACC  >  NC_000913/1974961‑1975140

gAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACa                                                                                         <  2:2351051/93‑1 (MQ=255)

 aGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACaa                                                                                        <  2:517614/93‑1 (MQ=255)

    gacgaAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGTCaataa                                                                                     <  2:4607762/93‑1 (MQ=255)

    gacgaAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACaataa                                                                                     >  1:3865880/1‑93 (MQ=255)

        aaaCAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCAc                                                                                 >  1:697017/1‑93 (MQ=255)

        aaaCAGCAACGGCAGCACCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCAc                                                                                 >  2:3183512/1‑93 (MQ=255)

          aCAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGCCGATAAACGACCCAATCCCTGCACCGGCAATAATCACCa                                                                               >  1:4040689/1‑93 (MQ=255)

          aCAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCa                                                                               <  2:147267/93‑1 (MQ=255)

           cAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAg                                                                              >  2:4370229/1‑93 (MQ=255)

           cAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAg                                                                              <  1:1479998/93‑1 (MQ=255)

           cAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAg                                                                              <  2:718690/93‑1 (MQ=255)

            aGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGt                                                                             <  2:53054/93‑1 (MQ=255)

             gCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGtt                                                                            >  2:1214792/1‑93 (MQ=255)

              cAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTc                                                                           <  1:1374642/93‑1 (MQ=255)

               aaCGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCa                                                                          >  2:3700778/1‑93 (MQ=255)

               aaCGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCa                                                                          >  2:3444817/1‑93 (MQ=255)

                  ggCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCg                                                                       <  1:1376937/93‑1 (MQ=255)

                  ggCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCg                                                                       <  1:1072934/93‑1 (MQ=255)

                   gCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCgg                                                                      >  1:939835/1‑93 (MQ=255)

                    cAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCggg                                                                     <  1:674969/93‑1 (MQ=255)

                     aGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGt                                                                    >  2:2351691/1‑93 (MQ=255)

                     aGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGt                                                                    >  2:3948631/1‑93 (MQ=255)

                        gcCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGtt                                                                   <  2:4342344/91‑1 (MQ=255)

                        gcCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGtt                                                                   >  1:4342344/1‑91 (MQ=255)

                        gcCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGa                                                                 <  1:2928659/93‑1 (MQ=255)

                        gcCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGa                                                                 >  1:2977409/1‑93 (MQ=255)

                         ccTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGAt                                                                <  2:3659517/93‑1 (MQ=255)

                         ccTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGAt                                                                >  1:2145004/1‑93 (MQ=255)

                           ttCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATag                                                              >  2:242410/1‑93 (MQ=255)

                            tCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATaga                                                             >  2:3911950/1‑93 (MQ=255)

                              aGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGt                                                           <  2:3320819/93‑1 (MQ=255)

                               gCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTg                                                          >  2:3967504/1‑93 (MQ=255)

                               gCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTg                                                          >  1:2078221/1‑93 (MQ=255)

                               gCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTg                                                          <  1:1586001/93‑1 (MQ=255)

                               gCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTg                                                          >  1:3361729/1‑93 (MQ=255)

                               gCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCACAGAGTTGATAGAGTg                                                          >  1:3361728/1‑93 (MQ=255)

                                 gTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCt                                                        >  1:2530156/1‑93 (MQ=255)

                                  ggCCTTTGATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTc                                                       >  2:2883050/2‑93 (MQ=255)

                                   gCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTcc                                                      >  2:3235814/1‑93 (MQ=255)

                                    ccTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCATGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCa                                                     >  1:3701190/1‑93 (MQ=255)

                                     cTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCaa                                                    <  1:1636066/93‑1 (MQ=255)

                                     cTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCaa                                                    >  1:4749646/1‑93 (MQ=255)

                                         aaTCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCt                                                <  1:2163614/93‑1 (MQ=255)

                                         aaTCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCt                                                >  2:4598828/1‑93 (MQ=255)

                                         aaTCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCt                                                >  1:4257331/1‑93 (MQ=255)

                                         aaTCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCt                                                >  2:3238639/1‑93 (MQ=255)

                                          aTCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCtt                                               <  2:902144/93‑1 (MQ=255)

                                           tCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTc                                              <  1:4118365/93‑1 (MQ=255)

                                             gCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCa                                            <  2:4464986/93‑1 (MQ=255)

                                             gCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCa                                            >  1:3153193/1‑93 (MQ=255)

                                               tttGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCAcc                                          >  2:1013237/1‑93 (MQ=255)

                                               tttGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCAcc                                          >  2:4160120/1‑93 (MQ=255)

                                               tttGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCAcc                                          >  1:587092/1‑93 (MQ=255)

                                                ttGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCg                                         <  2:196649/93‑1 (MQ=255)

                                                  tccATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGt                                       <  1:3588736/92‑1 (MQ=255)

                                                   ccATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtc                                      >  2:1938574/1‑93 (MQ=255)

                                                     attattGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtcat                                    >  1:3779646/1‑93 (MQ=255)

                                                      ttattGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtcatc                                   >  2:2373128/1‑93 (MQ=255)

                                                      ttattGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtcatc                                   >  2:3593831/1‑93 (MQ=255)

                                                       tattGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtcatca                                  >  2:210501/1‑93 (MQ=255)

                                                       tattGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGCGCTCCAAGGCTTCCACCGGtcatca                                  >  2:2600124/1‑93 (MQ=255)

                                                         ttGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCaaa                                <  1:241708/93‑1 (MQ=255)

                                                             cgacgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAc                            <  2:2892114/93‑1 (MQ=255)

                                                             cgacgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAc                            >  2:3154234/1‑93 (MQ=255)

                                                              gacgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAcc                           <  2:1823815/93‑1 (MQ=255)

                                                              gacgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAcc                           <  2:1656671/93‑1 (MQ=255)

                                                              gacgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAcc                           <  1:4652836/93‑1 (MQ=255)

                                                              gacgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAcc                           <  1:1469949/93‑1 (MQ=255)

                                                                cgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAccgc                         <  1:3830347/93‑1 (MQ=255)

                                                                 gaTAAAGGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAccgcc                        >  2:4668968/1‑93 (MQ=255)

                                                                 gaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAccgcc                        >  2:4100763/1‑93 (MQ=255)

                                                                 gaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAccgcc                        <  1:236235/93‑1 (MQ=255)

                                                                  aTAACCGACCCACTCCGTGCACCGACAATAATCACCAGTGCAGCGGGTTGATAGAGTGCTCCAAGGCTGCCACCGGTCATCAAATAAccgccg                       <  1:4549203/93‑1 (MQ=255)

                                                                  aTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAccgccg                       >  1:2503350/1‑93 (MQ=255)

                                                                  aTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAccgccg                       >  1:3599310/1‑93 (MQ=255)

                                                                  aTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAccgccg                       >  2:2107211/1‑93 (MQ=255)

                                                                      aCGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGagaa                   >  1:4172127/1‑93 (MQ=255)

                                                                      aCGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGaaaa                   >  1:2554606/1‑93 (MQ=255)

                                                                       cGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAAc                  <  1:2600136/93‑1 (MQ=255)

                                                                           ccAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTa              >  2:1832813/1‑93 (MQ=255)

                                                                            cAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAc             >  1:3348799/1‑93 (MQ=255)

                                                                            cAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAc             >  2:163875/1‑93 (MQ=255)

                                                                            cAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAc             >  1:1834393/1‑93 (MQ=255)

                                                                            cAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAc             <  1:581200/93‑1 (MQ=255)

                                                                            cAATCCCTCCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAc             <  2:4285393/93‑1 (MQ=255)

                                                                             aaTCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAcc            <  1:1005496/93‑1 (MQ=255)

                                                                             aaTCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAcc            >  2:615989/1‑93 (MQ=255)

                                                                              aTCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCg           >  1:1157342/1‑93 (MQ=255)

                                                                               tCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCga          <  1:4290321/93‑1 (MQ=255)

                                                                                cccTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACTgag         >  2:4184742/1‑93 (MQ=255)

                                                                                 cccacACCGACAATAACCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCTCCGGTAATCAAATACCCGCCGAAAACTGGACCgaga        <  2:1414344/89‑1 (MQ=255)

                                                                                  cTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaa       >  1:2050741/1‑93 (MQ=255)

                                                                                  cTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaa       <  1:4489766/93‑1 (MQ=255)

                                                                                   tGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaac      >  2:920849/1‑93 (MQ=255)

                                                                                      aCCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaacaac   >  1:932770/1‑93 (MQ=255)

                                                                                      aCCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaacaac   >  1:391574/1‑93 (MQ=255)

                                                                                      aCCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGCCATCCAATAACCACCGAAAACTGTACAGAGaacaac   >  2:1540469/1‑93 (MQ=255)

                                                                                      aCCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTACACCGGTCCTCAAATAACCGCCGAAAACTGTACCGCGaacaac   >  1:391689/1‑93 (MQ=255)

                                                                                       ccGCCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAAcc  <  2:2572190/93‑1 (MQ=255)

                                                                                       ccGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAAcc  <  1:2077997/93‑1 (MQ=255)

GAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACC  >  NC_000913/1974961‑1975140

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