BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	439,695	A→G	Y211H (TAT→CAT)	ispA ←	geranyltranstransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	439,695	0	A	G	98.9%	550.0 / ‑3.3	179	Y211H (TAT→CAT)	ispA	geranyltranstransferase
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (79/97); total (80/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 185 total aligned reads displayed.

CACCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTT  >  NC_000913/439605‑439784

cACCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCt                                                                                         <  1:402583/93‑1 (MQ=255)

cACCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCt                                                                                         >  1:4108759/1‑93 (MQ=255)

 aCCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCtt                                                                                        >  1:4288761/1‑93 (MQ=255)

  cccTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTg                                                                                       <  1:315761/93‑1 (MQ=255)

    cTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTc                                                                                     <  1:1565090/93‑1 (MQ=255)

    cTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTc                                                                                     >  1:4228118/1‑93 (MQ=255)

        cGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGTCAAGGCCGATGCTCTCTGCATGCTTGTCGAGt                                                                                 <  2:2419387/93‑1 (MQ=255)

         gTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTa                                                                                >  1:1913624/1‑93 (MQ=255)

         gTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGCCCGATGCTCTCTGCATGCTTGTCGAGTa                                                                                >  1:583918/1‑93 (MQ=255)

          ttttCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTAc                                                                               >  1:4234363/1‑93 (MQ=255)

           tttCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTAcc                                                                              <  2:692807/93‑1 (MQ=255)

            ttCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCg                                                                             >  1:2279104/1‑93 (MQ=255)

               ccAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCa                                                                          <  2:769474/93‑1 (MQ=255)

                cAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCag                                                                         >  1:4537348/1‑93 (MQ=255)

                     ttGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCac                                                                    <  2:1015449/93‑1 (MQ=255)

                     ttGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCac                                                                    >  1:1395835/1‑93 (MQ=255)

                      tGCAGGATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTCCCGGCAGAGCacg                                                                   <  1:778756/93‑1 (MQ=255)

                       gCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCacga                                                                  <  1:841086/93‑1 (MQ=255)

                        cAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCacgac                                                                 >  2:1052103/1‑93 (MQ=255)

                        cAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCacgac                                                                 <  1:1764238/93‑1 (MQ=255)

                          gTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGCAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTAGCGGCAGGGCcagacgc                                                               >  1:918382/1‑92 (MQ=255)

                          gTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGt                                                               >  1:2362505/1‑93 (MQ=255)

                          gTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGt                                                               <  2:2509327/93‑1 (MQ=255)

                          gTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGt                                                               <  2:2660961/93‑1 (MQ=255)

                           tATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTc                                                              <  2:818005/93‑1 (MQ=255)

                           tATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTc                                                              >  1:123525/1‑93 (MQ=255)

                            aTCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTcc                                                             >  1:3463015/1‑93 (MQ=255)

                             tctcCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCt                                                            >  1:2611459/1‑93 (MQ=255)

                              ctcCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCtt                                                           <  2:2392877/93‑1 (MQ=255)

                              ctcCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCtt                                                           >  2:2427627/1‑93 (MQ=255)

                              ctcCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCtt                                                           >  2:3187112/1‑93 (MQ=255)

                              ctcCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCtt                                                           >  2:2913228/1‑93 (MQ=255)

                               tcCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCttt                                                          >  1:546827/1‑93 (MQ=255)

                                 ccaccaCATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTCCCGGCAGAGCACGACGTCCTTTAt                                                        <  1:3202959/93‑1 (MQ=255)

                                 ccaccaCATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTAt                                                        >  2:3335427/1‑93 (MQ=255)

                                  caccacATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTAtc                                                       <  1:2275121/93‑1 (MQ=255)

                                  caccacATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTAtc                                                       <  1:2318896/93‑1 (MQ=255)

                                  caccacATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTAtc                                                       <  1:2882919/93‑1 (MQ=255)

                                    ccacATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTAtctc                                                     <  2:3704723/93‑1 (MQ=255)

                                     cacaTCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTcc                                                    >  1:760860/1‑93 (MQ=255)

                                     cacaTCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTcc                                                    <  2:3369397/93‑1 (MQ=255)

                                     cacaTCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTcc                                                    >  2:466058/1‑93 (MQ=255)

                                     cacaTCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTcc                                                    <  1:1623680/93‑1 (MQ=255)

                                      acaTCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCg                                                   >  2:547554/1‑93 (MQ=255)

                                      acaTCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCg                                                   >  1:3814166/1‑93 (MQ=255)

                                      acaTCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCg                                                   <  2:4096182/93‑1 (MQ=255)

                                      acaTCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCg                                                   >  1:2019135/1‑93 (MQ=255)

                                         tCCAGGATGTTATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGGCCTTTATCTCCGgct                                                >  2:2583534/1‑92 (MQ=255)

                                         tCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGgcg                                                <  1:1100834/93‑1 (MQ=255)

                                           cAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCt                                              >  1:2257231/1‑93 (MQ=255)

                                           cAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCt                                              >  2:1784662/1‑93 (MQ=255)

                                            aGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCtt                                             >  2:2757153/1‑93 (MQ=255)

                                              gATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGGCGTCCTTTATCTCCGGCGCTTaa                                           >  1:2479224/1‑93 (MQ=255)

                                              gATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTaa                                           <  1:2181067/93‑1 (MQ=255)

                                               aTGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAAt                                          <  2:1089717/93‑1 (MQ=255)

                                                 gTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGc                                        >  2:3081835/1‑93 (MQ=255)

                                                 gTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGc                                        <  2:3402529/93‑1 (MQ=255)

                                                  tCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCa                                       >  1:368638/1‑93 (MQ=255)

                                                    aTCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCAcc                                     <  2:3122042/93‑1 (MQ=255)

                                                    aTCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCAcc                                     >  2:3797180/1‑93 (MQ=255)

                                                    aTCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCAcc                                     <  1:3004279/93‑1 (MQ=255)

                                                     tCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCa                                    >  2:406415/1‑93 (MQ=255)

                                                     tCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCa                                    >  1:3313678/1‑93 (MQ=255)

                                                     tCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTAGAGTACCGGCAGAGCACGCCGTCCTTTATCTCCGGCGCTTACTGCACCa                                    >  1:3228662/1‑93 (MQ=255)

                                                      ccTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCaa                                   <  1:470018/93‑1 (MQ=255)

                                                       cTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAg                                  >  1:373418/1‑93 (MQ=255)

                                                        tGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAgg                                 >  2:4488056/1‑93 (MQ=255)

                                                        tGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAgg                                 >  2:4174997/1‑93 (MQ=255)

                                                        tGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAgg                                 <  2:307156/93‑1 (MQ=255)

                                                         gAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGc                                <  2:4338921/93‑1 (MQ=255)

                                                         gAACCTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGc                                >  2:1519577/1‑93 (MQ=255)

                                                         gAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGc                                <  1:834076/93‑1 (MQ=255)

                                                            ccTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGaa                             <  1:849861/93‑1 (MQ=255)

                                                            ccTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGaa                             <  2:3141373/93‑1 (MQ=255)

                                                             cTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAAc                            <  2:2340280/93‑1 (MQ=255)

                                                             cTGGAAGGCAAGGCCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAAc                            >  1:2554071/1‑93 (MQ=255)

                                                                  aggcaaggcCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAAcggcgg                       <  1:864969/93‑1 (MQ=255)

                                                                       aggcCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGa                  >  1:4391285/1‑93 (MQ=255)

                                                                       aggcCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGa                  >  2:1248106/1‑93 (MQ=255)

                                                                        ggcCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGaa                 <  1:4439557/93‑1 (MQ=255)

                                                                        ggcCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGaa                 >  1:4667523/1‑93 (MQ=255)

                                                                         gcCGATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAAt                <  1:4519010/93‑1 (MQ=255)

                                                                            gATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGACTCaa             >  1:1509824/1‑93 (MQ=255)

                                                                            gATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCaa             <  2:2963967/93‑1 (MQ=255)

                                                                            gATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCaa             >  2:51631/1‑93 (MQ=255)

                                                                            gATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCaa             <  2:1904182/93‑1 (MQ=255)

                                                                            gATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCaa             <  2:3596278/93‑1 (MQ=255)

                                                                            gATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCaa             <  2:3038816/93‑1 (MQ=255)

                                                                            gATGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCCCTTAATGCACCAAGGCGAACGGCGGCGCGAATCaa             <  2:4550215/93‑1 (MQ=255)

                                                                             aTGCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAAt            <  2:1116156/93‑1 (MQ=255)

                                                                               gCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCTTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATgc          <  2:1596366/93‑1 (MQ=255)

                                                                               gCTCTCTGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATgc          <  1:1955023/93‑1 (MQ=255)

                                                                                ctctctGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATgcg         <  1:1044331/93‑1 (MQ=255)

                                                                                  ctctGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGcc       <  2:12149/93‑1 (MQ=255)

                                                                                  ctctGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGcc       <  2:3287686/93‑1 (MQ=255)

                                                                                  ctctGCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGcc       <  1:3511302/93‑1 (MQ=255)

                                                                                      gCATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGtt   >  2:3981210/1‑93 (MQ=255)

                                                                                       cATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGttt  <  2:4278336/93‑1 (MQ=255)

                                                                                       cATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGttt  <  2:2106263/93‑1 (MQ=255)

                                                                                       cATGCTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGttt  <  2:609168/93‑1 (MQ=255)

CACCCTGGCGTTTTCCCAACGTTGCAGTATCTCCCACCACATCCAGGATGTCATCCTGAACCTGGAAGGCAAGGCCGATGCTCTCTGCATACTTGTCGAGTACCGGCAGAGCACGACGTCCTTTATCTCCGGCGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTT  >  NC_000913/439605‑439784

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