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.2%	466.4 / ‑5.3	132	A38V (GCC→GTC)	motA	MotA
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (71/60); total (71/61)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 136 total aligned reads displayed.

GAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACC  >  NC_000913/1974961‑1975140

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

     acgaAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACaataat                                                                                    <  2:3196828/93‑1 (MQ=255)

     acgaAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACaataat                                                                                    >  1:3967350/1‑93 (MQ=255)

      cgaAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATc                                                                                   >  2:3993522/1‑93 (MQ=255)

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

         aaCAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCAcc                                                                                <  2:2608867/93‑1 (MQ=255)

          aCAGCAACGGCAGCGCCTTCAGCTTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCa                                                                               <  2:4177745/93‑1 (MQ=255)

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

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

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

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

            aGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGt                                                                             >  2:1511917/1‑93 (MQ=255)

             gCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGtt                                                                            <  2:4184996/93‑1 (MQ=255)

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

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

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

              cAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTc                                                                           >  2:4553214/1‑93 (MQ=255)

               aaCGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCCGTTCa                                                                          >  1:3674328/1‑93 (MQ=255)

               aaCGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCa                                                                          <  2:3005389/93‑1 (MQ=255)

                aCGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAg                                                                         <  1:2348314/93‑1 (MQ=255)

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

                    cAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCggg                                                                     >  1:923419/1‑93 (MQ=255)

                     aGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGt                                                                    <  1:677731/93‑1 (MQ=255)

                     aGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGt                                                                    <  1:404889/93‑1 (MQ=255)

                     aGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGt                                                                    <  1:1793091/93‑1 (MQ=255)

                      gcgcCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGtt                                                                   >  2:3572493/1‑93 (MQ=255)

                       cgcCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCGAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGttt                                                                  >  1:3488452/1‑92 (MQ=255)

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

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

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

                         ccTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGAACCAATCCCTGCACCGACAATAATCACTAGTTCAGCGTGCTGAt                                                                >  2:507236/1‑93 (MQ=255)

                          cTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATa                                                               >  1:3256905/1‑93 (MQ=255)

                          cTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATa                                                               <  2:3776421/93‑1 (MQ=255)

                           ttCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATag                                                              <  1:3350613/93‑1 (MQ=255)

                           gtCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATag                                                              <  1:1516614/92‑1 (MQ=255)

                            tCAGCGTGCCTTTAATCGCTTTTCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATaga                                                             <  1:2735093/93‑1 (MQ=255)

                            tCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATaga                                                             <  1:182540/93‑1 (MQ=255)

                            tCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATaga                                                             <  1:1416435/93‑1 (MQ=255)

                             cAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATagag                                                            <  1:855595/93‑1 (MQ=255)

                                cGTGTCTTTAATCGCTTTGCCATTATTCCCGACGATAACCGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGc                                                         <  2:883634/93‑1 (MQ=255)

                                  tGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTc                                                       >  2:2406728/1‑93 (MQ=255)

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

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

                                    ccTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCa                                                     >  1:3574539/1‑93 (MQ=255)

                                     cTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCGCCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCaa                                                    >  2:4619034/1‑93 (MQ=255)

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

                                      tttAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAg                                                   >  2:1661304/1‑93 (MQ=255)

                                       ttAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAgg                                                  >  2:2749138/1‑93 (MQ=255)

                                            cGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTcc                                             <  1:73576/93‑1 (MQ=255)

                                              cTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCAc                                           <  1:1768205/93‑1 (MQ=255)

                                              cTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCAc                                           <  2:1479979/93‑1 (MQ=255)

                                              cTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCAc                                           >  2:2541757/1‑93 (MQ=255)

                                                 tGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCgg                                        >  2:1943183/1‑93 (MQ=255)

                                                   ccATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtc                                      <  2:2611345/93‑1 (MQ=255)

                                                    cATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtca                                     >  1:1129441/1‑93 (MQ=255)

                                                    cATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtca                                     <  1:2130663/93‑1 (MQ=255)

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

                                                      ttattGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtcatc                                   <  2:4308065/93‑1 (MQ=255)

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

                                                         ttGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCaaa                                >  2:813723/1‑93 (MQ=255)

                                                           gCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATa                              >  2:2552301/1‑93 (MQ=255)

                                                           gCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATa                              <  2:1625634/93‑1 (MQ=255)

                                                           gCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATa                              <  2:1577823/93‑1 (MQ=255)

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

                                                              gacgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGGTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAcc                           >  2:4610587/1‑93 (MQ=255)

                                                              cacgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAcc                           <  2:2130166/92‑1 (MQ=255)

                                                                cgaTAAACGACCCAATCCCTGCACCGACAATAATCGCCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAccgc                         <  1:3876178/93‑1 (MQ=255)

                                                                cgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGCTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAccgc                         >  2:1325355/1‑93 (MQ=255)

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

                                                                   tAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGa                      >  2:4251086/1‑93 (MQ=255)

                                                                   tAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGa                      >  2:1451322/1‑93 (MQ=255)

                                                                   tAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGa                      >  1:2976921/1‑93 (MQ=255)

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

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

                                                                      aCGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGaaaa                   <  2:3290372/93‑1 (MQ=255)

                                                                       cGACCCAATCCCTGCACCTACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAAc                  <  1:1418391/93‑1 (MQ=255)

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

                                                                       cGACACAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAAc                  <  1:929055/93‑1 (MQ=255)

                                                                        gACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACt                 <  1:1215735/93‑1 (MQ=255)

                                                                        gACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACt                 >  2:1372069/1‑93 (MQ=255)

                                                                          cccAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGt               >  1:1401630/1‑93 (MQ=255)

                                                                           ccAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTa              <  1:4487365/93‑1 (MQ=255)

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

                                                                             caccccTGCACCGCCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAcc            <  1:1137634/90‑1 (MQ=255)

                                                                             aaTCCCTGCACCAACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAcc            >  1:2532300/1‑93 (MQ=255)

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

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

                                                                               tCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCga          >  1:4102873/1‑93 (MQ=255)

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

                                                                               tCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCga          >  2:2270704/1‑93 (MQ=255)

                                                                                cccTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCgag         >  2:2416338/1‑93 (MQ=255)

                                                                                 ccTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCgaga        <  2:4370833/93‑1 (MQ=255)

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

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

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

                                                                                   tGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaac      <  2:3415896/93‑1 (MQ=255)

                                                                                     cACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaacaa    <  2:4061055/93‑1 (MQ=255)

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

                                                                                       ccGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAAcc  >  1:330160/1‑93 (MQ=255)

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

GAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACC  >  NC_000913/1974961‑1975140

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