BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,533,160	A→G	L116P (CTG→CCG)	envZ ←	protein histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,533,160	0	A	G	97.4%	472.5 / 0.5	156	L116P (CTG→CCG)	envZ	protein histidine kinase
Reads supporting (aligned to +/- strand): ref base A (0/4); new base G (76/76); total (76/80)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 164 total aligned reads displayed.

GTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAG  >  NC_000913/3533071‑3533249

gTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGtt                                                                                        <  2:2712635/93‑1 (MQ=255)

 tACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtg                                                                                       >  2:4381031/1‑93 (MQ=255)

  tcccAGATTTTTGGGGACAGCCAGGGTTTTAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgc                                                                                      <  1:1932674/92‑1 (MQ=255)

  aCCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgc                                                                                      >  2:2701974/1‑93 (MQ=255)

  aCCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTtgc                                                                                      <  1:3988936/93‑1 (MQ=255)

  aCCCAGACATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgc                                                                                      <  2:279804/93‑1 (MQ=255)

   cccAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgct                                                                                     <  1:2792405/93‑1 (MQ=255)

   cccAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgct                                                                                     >  2:4571019/1‑93 (MQ=255)

   cccAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgct                                                                                     >  1:4235108/1‑93 (MQ=255)

    ccAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgctg                                                                                    <  1:285620/93‑1 (MQ=255)

      aGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTgcg                                                                                  >  2:1042725/1‑93 (MQ=255)

       gATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTgcgc                                                                                 >  2:955088/1‑93 (MQ=255)

       gATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTCGCTgcgc                                                                                 >  1:646239/1‑93 (MQ=255)

        atatTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGCGCCGCCCGGTTGCTGCGcc                                                                                >  2:1144706/1‑93 (MQ=255)

          atTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCAt                                                                              <  2:2943405/93‑1 (MQ=255)

          atTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCAt                                                                              >  1:1189186/1‑93 (MQ=255)

           ttGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATc                                                                             <  1:523598/93‑1 (MQ=255)

           ttGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATc                                                                             <  1:4093412/93‑1 (MQ=255)

             gggCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTg                                                                           <  1:714701/93‑1 (MQ=255)

                 gACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATgg                                                                       >  1:3814259/1‑93 (MQ=255)

                  aCAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGc                                                                      <  1:4218977/93‑1 (MQ=255)

                    aGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCtt                                                                    >  1:660981/1‑93 (MQ=255)

                      ccAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTaa                                                                  >  1:2270566/1‑93 (MQ=255)

                        aGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGa                                                                <  1:1586281/93‑1 (MQ=255)

                         ggTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGCTGGCTTAAGaa                                                               >  1:4176280/1‑93 (MQ=255)

                         ggTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGaa                                                               <  1:4422108/93‑1 (MQ=255)

                         ggTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGaa                                                               <  1:3351183/93‑1 (MQ=255)

                            tttCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTc                                                            <  2:1675526/93‑1 (MQ=255)

                               cAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATa                                                         <  2:1049460/93‑1 (MQ=255)

                                aGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAg                                                        >  1:3422951/1‑93 (MQ=255)

                                 gCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgt                                                       >  2:1167814/1‑93 (MQ=255)

                                  ccAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgtg                                                      <  1:987491/93‑1 (MQ=255)

                                   cAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgtgt                                                     >  2:2922915/1‑93 (MQ=255)

                                   cAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgtgt                                                     >  1:3084730/1‑93 (MQ=255)

                                   cAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgtgt                                                     <  2:783846/93‑1 (MQ=255)

                                    aGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGtt                                                    <  1:2035828/93‑1 (MQ=255)

                                    aGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGtt                                                    <  1:1681714/93‑1 (MQ=255)

                                    aGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGtt                                                    >  2:23177/1‑93 (MQ=255)

                                    aGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGtt                                                    <  2:2360809/93‑1 (MQ=255)

                                    aGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGtt                                                    <  2:583347/93‑1 (MQ=255)

                                    aGACGAAAGGCGAACATATGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCATAAGAATTCATAGTAtt                                                    <  2:4211641/93‑1 (MQ=255)

                                     gacgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTg                                                   >  2:1574142/1‑93 (MQ=255)

                                     gacgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTg                                                   >  1:572771/1‑93 (MQ=255)

                                     gacgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTg                                                   >  1:2938374/1‑93 (MQ=255)

                                       cgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTgcg                                                 >  2:1364030/1‑93 (MQ=255)

                                        gacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTgcgc                                                <  1:1013150/93‑1 (MQ=255)

                                         acAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGcc                                               >  2:4697049/1‑93 (MQ=255)

                                           aGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCa                                             >  1:2654302/1‑93 (MQ=255)

                                           aGGCGAACTGTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCa                                             <  2:337638/93‑1 (MQ=255)

                                             gCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAAc                                           <  2:290996/93‑1 (MQ=255)

                                               gAACTTTTGTTGACCTCAACGCGCTCTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGc                                         <  2:809422/93‑1 (MQ=255)

                                               gAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGc                                         <  2:3287256/93‑1 (MQ=255)

                                                aaCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCa                                        <  1:1068196/93‑1 (MQ=255)

                                                aaCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCa                                        >  1:3166578/1‑93 (MQ=255)

                                                 aCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAg                                       >  1:2574468/1‑93 (MQ=255)

                                                 aCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTCAAGAATTCATAGTGTTGCGCCCAACGCAg                                       >  1:1734406/1‑93 (MQ=255)

                                                      tgttgACCTCAACGCGCTCGTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTg                                  <  1:1516749/93‑1 (MQ=255)

                                                       gttgGCCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGc                                 <  2:1292707/93‑1 (MQ=255)

                                                       gttgACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGc                                 <  1:2747548/93‑1 (MQ=255)

                                                       gttgACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGc                                 <  2:4736267/93‑1 (MQ=255)

                                                          gACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCctc                              >  1:1240420/1‑93 (MQ=255)

                                                          gACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCctc                              >  1:3578679/1‑93 (MQ=255)

                                                           aCCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCctct                             >  1:716077/1‑93 (MQ=255)

                                                           aCCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCAACTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCctct                             <  2:1551374/93‑1 (MQ=255)

                                                            ccTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCtt                            <  1:15165/93‑1 (MQ=255)

                                                             cTCAACGCGCAGTTCCGTCGTGCCGCCCGGTTGCTGCGCAATCTGATGTCTTAAGAATTCATAGTGTTGCGCCCAACTCAGCCCTGCCTCTTc                           <  2:1790966/93‑1 (MQ=255)

                                                              tCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCg                          >  2:558370/1‑93 (MQ=255)

                                                              tCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCg                          >  1:4670534/1‑93 (MQ=255)

                                                              tCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCg                          >  1:373088/1‑93 (MQ=255)

                                                                aaCGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGc                        >  2:847316/1‑93 (MQ=255)

                                                                 aCGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCa                       <  1:645420/93‑1 (MQ=255)

                                                                  cgcgCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAg                      >  1:2197271/1‑93 (MQ=255)

                                                                   gcgcACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGc                     >  1:1666717/1‑93 (MQ=255)

                                                                    cgcACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGcc                    <  1:296470/93‑1 (MQ=255)

                                                                    cgcACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGcc                    >  1:4097145/1‑93 (MQ=255)

                                                                    cgcACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGcc                    >  2:3552417/1‑93 (MQ=255)

                                                                    cgcACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGcc                    >  2:3552418/1‑93 (MQ=255)

                                                                     gcACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCt                   >  2:3778604/1‑93 (MQ=255)

                                                                      cACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTc                  >  1:898920/1‑93 (MQ=255)

                                                                      cACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTc                  <  1:2039979/93‑1 (MQ=255)

                                                                      cACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTc                  >  2:1980334/1‑93 (MQ=255)

                                                                      cACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTc                  <  2:174281/93‑1 (MQ=255)

                                                                       aCTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCg                 >  1:4182823/1‑93 (MQ=255)

                                                                        cTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTTAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGt                >  2:479235/1‑93 (MQ=255)

                                                                        cTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGt                <  1:2549596/93‑1 (MQ=255)

                                                                        cTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGt                >  2:189054/1‑93 (MQ=255)

                                                                         ttCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGtt               >  1:318962/1‑93 (MQ=255)

                                                                           ccGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTgg             >  2:2959562/1‑93 (MQ=255)

                                                                           ccGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTgg             <  1:3212954/93‑1 (MQ=255)

                                                                            cGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGa            <  2:3118051/93‑1 (MQ=255)

                                                                             gTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAg           <  1:3222354/93‑1 (MQ=255)

                                                                                gggCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTag        <  1:3587537/93‑1 (MQ=255)

                                                                                gggCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTag        <  1:4167224/93‑1 (MQ=255)

                                                                                gggCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTag        >  1:894389/1‑93 (MQ=255)

                                                                                 ggCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTaga       >  2:4042354/1‑93 (MQ=255)

                                                                                 ggCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTaga       <  2:2302939/93‑1 (MQ=255)

                                                                                  gccgccCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagag      <  2:3272511/93‑1 (MQ=255)

                                                                                   ccgccCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagaga     <  1:2420452/93‑1 (MQ=255)

                                                                                    cgccCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagagag    <  2:4389160/93‑1 (MQ=255)

                                                                                    cgccCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagagag    >  2:1387998/1‑93 (MQ=255)

GTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAG  >  NC_000913/3533071‑3533249

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