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.6%	832.3 / ‑5.5	233	A38V (GCC→GTC)	motA	MotA
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (119/113); total (120/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

Only 100 of 238 total aligned reads displayed.

GAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACC  >  NC_000913/1974961‑1975140

gAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACa                                                                                         >  1:555595/1‑93 (MQ=255)

 aGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCCTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACaa                                                                                        >  2:109907/1‑93 (MQ=255)

    gacgaAACAGCAACGGCAGCTCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACaataa                                                                                     <  2:1764473/93‑1 (MQ=255)

    gacgaAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACaataa                                                                                     <  1:4246393/93‑1 (MQ=255)

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

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

       gaAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCa                                                                                  >  2:2738515/1‑93 (MQ=255)

       gaAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCa                                                                                  >  1:7857759/1‑93 (MQ=255)

        aaaCAGCAGCGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCAc                                                                                 <  2:5577242/93‑1 (MQ=255)

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

        aaaCAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCAc                                                                                 <  1:7083065/93‑1 (MQ=255)

          aCAGCAGCGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCa                                                                               <  2:4550299/93‑1 (MQ=255)

          aCAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCa                                                                               >  2:2231379/1‑93 (MQ=255)

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

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

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

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

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

                aCGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAg                                                                         >  1:1964597/1‑93 (MQ=255)

                aCGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAg                                                                         >  2:154335/1‑93 (MQ=255)

                  ggCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCg                                                                       >  1:2177031/1‑93 (MQ=255)

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

                     aGCGCCTTCAGCGTGCCTTTAATCGCTTTTCCATTATTGCCGACGATCAACGACCCAATCCTTGCACCGCCAATAATCACCAGTTCAGCGCGt                                                                    >  1:1012992/1‑93 (MQ=255)

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

                       cgcCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTg                                                                  <  1:6950480/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

                            tCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGCTCAGCGGGTTGATaga                                                             >  1:1292237/1‑93 (MQ=255)

                              aGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGt                                                           >  1:6809456/1‑93 (MQ=255)

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

                               gCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTAATAGAGTg                                                          <  1:1548732/93‑1 (MQ=255)

                                cGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGGTGATAGAGAGc                                                         >  1:6369138/1‑93 (MQ=255)

                                 gTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCt                                                        <  2:6735671/93‑1 (MQ=255)

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

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

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

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

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

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

                                      tttAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCAACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAg                                                   >  2:7234022/1‑93 (MQ=255)

                                       ttAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAgg                                                  <  1:1853469/93‑1 (MQ=255)

                                        tAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGc                                                 <  2:7194076/93‑1 (MQ=255)

                                        tAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGc                                                 <  2:6794564/93‑1 (MQ=255)

                                        tAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGc                                                 <  2:6075065/93‑1 (MQ=255)

                                         aaTCTCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCt                                                <  2:5904928/93‑1 (MQ=255)

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

                                           tCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTc                                              >  2:5617006/1‑93 (MQ=255)

                                           tCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTc                                              >  2:1376685/1‑93 (MQ=255)

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

                                            cGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTcc                                             >  1:2414075/1‑93 (MQ=255)

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

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

                                               tttGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCAcc                                          <  1:6581221/93‑1 (MQ=255)

                                                  gCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGt                                       >  1:5307024/1‑93 (MQ=255)

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

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

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

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

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

                                                      ttagtGCCGACCATAAACGGCCCAAGCCCTTCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtcatc                                   <  2:5906869/93‑1 (MQ=255)

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

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

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

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

                                                             cgacgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCCCCAAGGCTTCCACCGGTCATCAAATAAc                            <  2:4176040/93‑1 (MQ=255)

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

                                                               acgaTAAACGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAAccg                          >  1:7728925/1‑93 (MQ=255)

                                                                     aaCGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGaaa                    >  1:7304526/1‑93 (MQ=255)

                                                                     aaCGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGaaa                    >  1:4112412/1‑93 (MQ=255)

                                                                     aaCGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGaaa                    <  2:5653362/93‑1 (MQ=255)

                                                                     aaCGACCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGaaa                    <  1:3337510/93‑1 (MQ=255)

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

                                                                         aCCCAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTg                >  1:5131025/1‑93 (MQ=255)

                                                                          cccAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTga               >  1:5754359/1‑92 (MQ=255)

                                                                          cccAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGt               <  2:3105589/93‑1 (MQ=255)

                                                                          cccAATCCCTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGt               <  1:3569280/93‑1 (MQ=255)

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

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

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

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

                                                                            cAATCCCTACACCGACAATTATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGAAc             <  2:7025584/93‑1 (MQ=255)

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

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

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

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

                                                                                 ccTGCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCgaga        >  1:2864807/1‑93 (MQ=255)

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

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

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

                                                                                    gCACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaaca     <  1:1149772/93‑1 (MQ=255)

                                                                                     cACCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAAGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaacaa    >  1:4385571/1‑93 (MQ=255)

                                                                                      aCCGACAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaacaac   <  1:4550057/93‑1 (MQ=255)

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

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

GAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACC  >  NC_000913/1974961‑1975140

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