BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,953,910	A→G	N320S (AAC→AGC)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,953,910	0	A	G	98.3%	540.3 / ‑3.2	175	N320S (AAC→AGC)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (95/77); total (97/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 183 total aligned reads displayed.

CGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAACTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAG  >  NC_000913/3953820‑3953999

cgctggcgctggCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCt                                                                                         <  1:921282/93‑1 (MQ=255)

cgctggcgctggCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCt                                                                                         <  2:2439918/93‑1 (MQ=255)

  ctggcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTc                                                                                       >  2:3872057/1‑93 (MQ=255)

   tggcgctggcgGGAATGCACAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTcc                                                                                      <  2:4223144/93‑1 (MQ=255)

    ggcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCa                                                                                     >  2:2938714/1‑93 (MQ=255)

    ggcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAACTTCCa                                                                                     >  2:2341525/1‑93 (MQ=255)

       gctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACgg                                                                                  <  2:3738914/93‑1 (MQ=255)

       gctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCAAACGTGAGCTTCCACgg                                                                                  >  1:133937/1‑93 (MQ=255)

       gctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCCGTGCCCACGTGAGCTTCCAccg                                                                                  >  2:904834/1‑91 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                >  1:1212494/1‑93 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                <  2:863146/93‑1 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                >  2:1693897/1‑93 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                <  1:1031584/93‑1 (MQ=255)

         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                <  1:586395/93‑1 (MQ=255)

           gcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTg                                                                              >  2:4268179/1‑93 (MQ=255)

           gcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTg                                                                              >  1:4403810/1‑93 (MQ=255)

            cgGGAATTAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTTCACGGCCTgc                                                                             >  2:4564615/1‑93 (MQ=255)

            cgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTgc                                                                             <  1:2860201/93‑1 (MQ=255)

                aaTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTa                                                                         >  1:1626928/1‑93 (MQ=255)

                aaTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTa                                                                         >  1:464426/1‑93 (MQ=255)

                 aTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTAc                                                                        >  2:4447755/1‑93 (MQ=255)

                  tGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACg                                                                       >  1:4624431/1‑93 (MQ=255)

                  tGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACg                                                                       >  2:2186566/1‑93 (MQ=255)

                    aaaaaaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGtc                                                                     >  2:2007651/1‑93 (MQ=255)

                       aaaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCa                                                                  >  1:3257373/1‑93 (MQ=255)

                       aaaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGGGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCa                                                                  >  2:4543454/1‑93 (MQ=255)

                        aaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTATTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAg                                                                 >  1:1980493/1‑93 (MQ=255)

                          tatatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGaa                                                               <  1:2975625/93‑1 (MQ=255)

                           atatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAAc                                                              >  2:1205642/1‑93 (MQ=255)

                           atatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAAc                                                              <  1:1042158/93‑1 (MQ=255)

                            tatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACg                                                             >  1:3158164/1‑93 (MQ=255)

                             atCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCCGAACGc                                                            <  2:4608121/93‑1 (MQ=255)

                                gCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGc                                                         <  2:3713247/93‑1 (MQ=255)

                                gCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGc                                                         >  1:2830673/1‑93 (MQ=255)

                                gCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGc                                                         >  1:1758533/1‑93 (MQ=255)

                                gCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGc                                                         >  2:4683424/1‑93 (MQ=255)

                                  ccTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGa                                                       >  1:1289180/1‑93 (MQ=255)

                                  ccTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGa                                                       >  2:2639818/1‑93 (MQ=255)

                                  ccTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGa                                                       >  2:1307334/1‑93 (MQ=255)

                                   cTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGaa                                                      >  1:4589274/1‑93 (MQ=255)

                                   cTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGaa                                                      <  1:2587863/93‑1 (MQ=255)

                                     gCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACt                                                    >  1:4179113/1‑93 (MQ=255)

                                      cacaACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTg                                                   >  2:2351880/1‑93 (MQ=255)

                                       acaacaTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTgg                                                  >  1:455366/1‑93 (MQ=255)

                                        caacaTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTggg                                                 <  2:1460100/93‑1 (MQ=255)

                                         aacaTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTTCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGc                                                >  1:414910/1‑93 (MQ=255)

                                            aTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGaa                                             >  2:3769392/1‑93 (MQ=255)

                                            aTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGaa                                             <  1:3303311/93‑1 (MQ=255)

                                             ttCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAAc                                            <  2:2767801/93‑1 (MQ=255)

                                             ttCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAAc                                            >  1:4149599/1‑93 (MQ=255)

                                                gcggcgAACGGCTTGCGCATATTCTTTCCTGTGCCAACGTGAGCTTCCACGGCCTCCGCTACGTCTCAGAACCCCGCGAAATGGGCGcacacc                                         >  2:374317/1‑91 (MQ=255)

                                                gcggcgAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGc                                         <  2:795286/93‑1 (MQ=255)

                                                gcggcgAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGc                                         <  1:341753/93‑1 (MQ=255)

                                                     gAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGaa                                    >  1:1779652/1‑93 (MQ=255)

                                                     gAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGaa                                    >  1:4288552/1‑93 (MQ=255)

                                                     gAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGaa                                    <  2:4198514/93‑1 (MQ=255)

                                                     gAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCCGCGAACTGGGCGAACAGCGTGaa                                    >  2:3058224/1‑93 (MQ=255)

                                                      aaCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCTCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAg                                   <  2:3616206/93‑1 (MQ=255)

                                                      aaCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAg                                   <  1:2370804/93‑1 (MQ=255)

                                                      aaCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAg                                   <  1:1866899/93‑1 (MQ=255)

                                                      aaCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAg                                   >  2:2463731/1‑93 (MQ=255)

                                                      aaCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAg                                   >  1:1149443/1‑93 (MQ=255)

                                                      aaCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAg                                   <  1:1104440/93‑1 (MQ=255)

                                                       aCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGc                                  >  1:703564/1‑93 (MQ=255)

                                                       aCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAACGTGAAgg                                  >  1:1548118/1‑92 (MQ=255)

                                                        cGGCTGGCTCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCGGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCg                                 <  2:2812596/93‑1 (MQ=255)

                                                         ggCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGCAGCgt                                >  2:3934779/1‑93 (MQ=255)

                                                          gCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgtt                               <  1:2429710/93‑1 (MQ=255)

                                                          gCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgtt                               <  1:1415709/93‑1 (MQ=255)

                                                          gCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgtt                               <  2:1417829/93‑1 (MQ=255)

                                                           cTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttg                              >  2:731893/1‑93 (MQ=255)

                                                           cTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttg                              <  1:658447/93‑1 (MQ=255)

                                                            tGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgt                             <  2:610989/93‑1 (MQ=255)

                                                                gcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTggg                         >  1:734326/1‑92 (MQ=255)

                                                                gcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGc                         <  1:2731780/93‑1 (MQ=255)

                                                                 cATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCg                        >  1:3666713/1‑93 (MQ=255)

                                                                  atatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCgg                       >  1:2719899/1‑93 (MQ=255)

                                                                  atatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCgg                       >  1:524701/1‑93 (MQ=255)

                                                                   tatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCggg                      >  1:1709278/1‑92 (MQ=255)

                                                                   tatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGt                      >  1:2559036/1‑93 (MQ=255)

                                                                   tatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTCGGCGAACAGCGTGAAGCGTTGTTGGCGGt                      >  1:3223390/1‑93 (MQ=255)

                                                                    atTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTg                     <  2:1865229/93‑1 (MQ=255)

                                                                      tCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGAc                   >  2:1502682/1‑93 (MQ=255)

                                                                      tCTTTCCGGTGCCAACGTGAACTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGAc                   >  1:3090103/1‑93 (MQ=255)

                                                                        tttCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCa                 >  2:2547476/1‑93 (MQ=255)

                                                                        tttCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCa                 <  2:4317351/93‑1 (MQ=255)

                                                                         ttCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCAt                <  2:2017274/93‑1 (MQ=255)

                                                                           ccGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTc              <  2:4152740/93‑1 (MQ=255)

                                                                            cGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTcc             <  2:4368483/93‑1 (MQ=255)

                                                                            cGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTcc             >  1:3787493/1‑93 (MQ=255)

                                                                              gTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTTACCATTCCgg           >  1:2032742/1‑93 (MQ=255)

                                                                              gTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCgg           >  2:593568/1‑93 (MQ=255)

                                                                                gCCAACGTGAACTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGgaa         <  1:2487299/93‑1 (MQ=255)

                                                                                   aaCGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAAGAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGgaagaa      >  2:2088956/1‑93 (MQ=255)

                                                                                    aCGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaa     <  1:3557525/93‑1 (MQ=255)

                                                                                    aCGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaa     >  2:3831833/1‑93 (MQ=255)

                                                                                     cGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaaa    <  2:2606922/93‑1 (MQ=255)

                                                                                      gTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaaaa   >  2:1798053/1‑93 (MQ=255)

                                                                                       tGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAg  >  2:4126796/1‑93 (MQ=255)

                                                                                       tGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAg  <  2:3358238/93‑1 (MQ=255)

CGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAACTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAG  >  NC_000913/3953820‑3953999

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